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I.  PREFACE 

The  nearshore  currents  generated  by  obliquely  incident  breaking  waves  within  the 
surf  zone  can  exceed  1  ms"1  during  high  energy  conditions,  and  are  some  of  the  strongest 
currents  with  the  greatest  shear  in  all  the  oceans.  These  strong  currents  play  a  major  role  in 
oil  spills  and  dispersion  of  other  pollutants,  and  also  transport  significant  amounts  of 
sediment,  having  a  dramatic  effect  on  beach  erosion.  Additionally,  the  modeling  of 
nearshore  circulation  assumes  a  special  significance  for  military  operations,  where  the 
success  or  failure  of  an  amphibious  assault,  or  a  special  forces  operation  is  directly  related 
to  our  ability  to  forecast  the  conditions  within  this  unique  operational  theater. 

The  physical  understanding  and  mathematical  modeling  of  hydrodynamic  processes 
in  the  nearshore  zone  have  substantially  improved  during  the  last  two  decades.  Yet,  most 
existing  models  are  formulated  based  on  a  monochromatic  wave  description  and  have  been 
tested  and  calibrated  with  laboratory  data,  due  to  difficulties  associated  with  measurements 
within  this  harsh  environment.  In  this  study,  monochromatic  formulations  are  extended  to 
random  wave  fields  by  using  a  probabilistic  approach  (Thornton  and  Guza,  1983).  The 
numerical  models  developed  here  are  compared  with  data  acquired  during  the  DUCK94 
(October  1994)  experiment  conducted  on  a  barred  ocean  beach  at  the  U.S.  Army  Corps  of 
Engineers  Field  Research  Facility  near  Duck,  North  Carolina.  The  objective  of  this 
dissertation  is  to  improve  our  current  understanding  of  the  vertical  and  cross-shore 
distributions  of  both  longshore  and  cross-shore  currents,  providing  the  first  steps  towards  a 
comprehensive  three-dimensional  model. 

In  Chapter  II,  observed  vertical  structure  of  the  mean  longshore  current  over  a  barred 


beach  is  examined  to  test  the  hypothesis  that  the  turbulent  bottom  boundary  layer  of  the  mean 
longshore  currents  is  logarithmic  and  to  investigate  the  influence  of  wave-breaking  generated 
turbulence  on  the  mean  longshore  current  profile.  The  relationship  between  bottom  shear 
stress  and  bottom  roughness,  including  the  influence  of  ripples  and  mega-ripples  is  also 
investigated.  This  Chapter  consists  of  a  journal  article  accepted  by  the  Journal  of 
Geophysical  Research. 

The  vertical  structure  of  the  mean  cross-shore  current  (undertow)  over  a  barred  beach 
is  investigated  in  Chapter  HI.  A  surface  mass  flux  model,  including  contributions  from  white 
water  surface  rollers  on  breaking  wave  crests  (Lippmann  and  Thornton,  1997),  is  used  to 
predict  the  cross-shore  distribution  of  depth-averaged  undertow.  The  impact  of  using  linear 
and  non-linear  wave  theories  to  estimate  surface  mass  flux  is  investigated.  The  vertical 
structure  of  the  undertow  is  modeled  using  an  eddy  viscosity  closure  to  estimate  the  turbulent 
shear  stress.  The  effects  of  different  depth-dependent  formulations  for  the  eddy  viscosity  as 
well  different  boundary  conditions  on  the  vertical  structure  of  the  undertow  are  evaluated. 

In  Chapter  IV,  the  models  of  the  vertical  profiles  of  longshore  and  cross-shore  mean 
currents  developed  in  the  preceding  Chapters  are  combined  to  formulate  a  quasi  three- 
dimensional  model  to  describe  the  cross-shore  distribution  of  the  longshore  current.  This 
model  includes  turbulent  mixing  due  to  the  cross-shore  advection  of  mean  momentum  of  the 
longshore  current  by  the  mean  cross-shore  current  (Putrevu  and  Svendsen,  1994),  and 
contributions  from  breaking  wave  rollers.  The  effect  of  introducing  the  observed  cross-shore 
variation  of  bottom  shear  stress  (Chapter  I)  on  the  longshore  current  prediction  is 
investigated. 


Each  of  the  three  dissertation  Chapters  represents  a  separate  paper  for  publication. 
As  such,  they  are  self-contained  with  their  own  introduction,  theory  description,  data  results, 
conclusions,  and  references. 


II.  VERTICAL  PROFILES  OF  LONGSHORE  CURRENTS  AND  RELATED  BED 
SHEAR  STRESS  AND  BOTTOM  ROUGHNESS 


(This  chapter  consists  of  a  journal  article  accepted  by 
the  Journal  of  Geophysical  Research  in  May,  1997) 


VERTICAL  PROFILES  OF  LONGSHORE  CURRENTS  AND  RELATED  BED 

SHEAR  STRESS  AND  BOTTOM  ROUGHNESS 

A.F.  Garcez  Faria,  E.B.  Thornton,  T.P.  Stanton,  C.V.  Soares1  and  T.C.  Lippmann2 

Naval  Postgraduate  School,  Monterey,  California  93940 

ABSTRACT 
The  vertical  structure  of  the  mean  wave-driven  longshore  current  over  a  barred  beach 
is  examined  on  three  strong  current  days  during  the  DUCK94  experiment,  and  it  is  found  that 
the  bottom  boundary  layer  is  well  described  by  a  logarithmic  profile  (mean  correlation 
coefficient  for  all  22  profiles,  0.98).  The  logarithmic  profile  fits  better  in  the  trough  where 
turbulent  bottom  boundary  layer  processes  predominate,  than  over  the  bar  where  breaking- 
wave  induced  turbulence  generated  at  the  surface  modifies  the  profile.  The  surface  layer  in 
the  presence  of  waves  is  well  described  by  adjusting  the  logarithmic  profile  for  the 
intermittent  presence  of  water  and  adding  the  alongshore  component  of  the  mass  transport 
velocity  (slope  of  the  least-squares  linear  regression  between  model  predictions  and 
observations,  1 .005  and  rms  error  of  7  percent). 

Bed  shear  stresses  calculated  from  logarithmic  velocity  profiles  are  equated  to  a 
quadratic  bottom  shear  stress  formulation.  The  associated  bed  shear  stress  coefficients  vary 
by  more  than  an  order  of  magnitude  across  the  surf  zone  (0.0006-0.012).  Bottom  roughness 
was  measured  throughout  the  nearshore  using  a  sonic  altimeter  mounted  on  a  moving 
platform.  The  bed  shear  stress  coefficients  are  positively  correlated  with  bottom  roughness 


'Presently  at  Marinha  -  Instituto  Hidrografico,  1296  Lisboa  Codex,  Portugal 
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(linear  correlation  coefficient,  0.6).  A  higher  linear  correlation  coefficient  (0.8)  is  obtained 
by  subtracting  skin  friction  from  the  total  bed  shear  stress. 

INTRODUCTION 

Knowledge  of  the  bottom  boundary  and  surface  layers  is  fundamental  to 
understanding  nearshore  hydrodynamics  and  sediment  processes.  For  steady  flow,  such  as 
in  a  river,  the  bottom  boundary  layer  is  well  described  by  a  logarithmic  profile.  Only  limited 
wave-driven,  longshore  current  vertical  profiles  have  been  measured.  Visser  (1986)  measured 
wave-driven  longshore  currents  in  a  laboratory  experiment  using  micro-propeller  and  laser 
Doppler  velocimeters  and  found  profiles  approached  a  logarithmic  form.  In  a  similar 
laboratory  wave-driven  longshore  current  experiment  using  laser  velocimeters,  Simons  et 
al.  (1992)  verified  that  the  vertical  profiles  tended  to  be  logarithmic. 

The  velocity  profile  of  a  steady  current  is  modified  by  the  presence  of  waves.  The 
superposition  of  waves  on  the  mean  current  produces  enhanced  bottom  friction  (e.g.,  Grant 
and  Madsen,  1979;  Christoffersen  and  Jonsson,  1985;  Myrhaug  and  Slaattelid,  1989;  and 
Sleath,  1990 ).  As  a  result,  the  vertical  gradient  of  the  mean  current  near  the  bed  is  increased 
and  a  more  uniform  profile  can  be  expected  throughout  most  of  the  water  column.  The 
influence  of  waves  is  inversely  proportional  to  water  depth  and  wave  frequency,  with 
decreasing  importance  of  nonlinear  interactions  between  waves  and  currents  within  the 
wave-bottom  boundary  layer  with  increasing  depth  and  frequency. 

Turbulence  induced  by  breaking  waves  modifies  the  vertical  profile  of  longshore 
currents.  The  downward  momentum  mixing  produced  from  wave  breaking-injected 
turbulence  results  in  a  more  uniform  velocity  profile.  Fredsoe  and  Deigaard  (1992),  and 


Church  and  Thornton  (1993)  argue  this  would  result  in  an  increased  bed  shear  stress. 

In  the  following,  mean  longshore  current  profiles  obtained  over  a  barred  beach  are 
examined  with  the  objectives  of:  1 )  testing  the  hypothesis  that  the  turbulent  bottom  boundary 
layer  of  the  mean  longshore  currents  is  logarithmic,  2)  investigating  the  influence  of  the 
surface  layer  on  the  mean  longshore  current  profile,  and  3)  examining  the  relationship 
between  bottom  shear  stress  and  bottom  roughness,  including  the  influence  of  ripples  and 
mega-ripples. 

THEORY 

The  vertical  profile  of  longshore  currents  is  significantly  affected  by  the  bottom 
boundary  and  surface  layers.  The  bottom  boundary  layer  assumes  more  importance  since  it 
determines  the  general  logarithmic  profile  shape  over  most  of  the  water  column.  However, 
processes  in  the  surface  layer  can  modify  the  profile  in  the  presence  of  waves  and  wind. 
Bottom  Boundary  Layer 

Neglecting  molecular  viscous  stresses,  the  alongshore  momentum  equation  (  y- 
direction)  is  written 


dpv        dpvu        dpv 2        dpwv        _  dp 
dt  dx  dy  dz  dy 


The  velocities  are  expanded  into  mean,  turbulent,  and  wave-induced  components, 
w;.=  U .+  w+  ui and  w  =  w+  w,  where  (i  =  1,2)  refers  to  horizontal  coordinates  (x,y)  and 
the  mean  vertical  velocity  is  assumed  equal  to  zero.  After  time  averaging,  (1)  can  be 
simplified  with  the  aid  of  the  following  assumptions:  1)  straight  and  parallel  contours,  e.g., 
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3  -\ 

—  (    )  =0  (overbar  indicates  time  averaging);  2)  steady  state  conditions,  i.e.,  — (    )  =0; 

dy  dt 

3)    wave-induced    and   turbulent    velocity   components    are    statistically    independent 
(uncorrelated),  i.e.,  wv  =  0;  4)  horizontal  turbulent  momentum  flux  is  small  compared  to 


wave  induced  momentum  flux  (Stive  and  Wind,  1982),  i.e.,— - —  «     ^vu  ,  and  can  be 

dx  dx 

neglected.  Contributions  due  to  momentum  mixing  caused  by  interactions  between  cross- 


shore  and  longshore  currents  (— )are  usually  not  negligible  (Svendsen  and  Putrevu, 

dx 

1994;  and  Garcez  Faria  et  ah,  1996).  Nevertheless,  they  are  neglected  here  for  simplicity 
sake.  Applying  these  assumptions,  the  alongshore  momentum  equations  can  be  written 


dpwv        dpvu        dpwv 

— —  =  -±—  +  — —  (2) 

dZ  dx  dZ 


which  says  that  the  sum  of  cross-shore  changes  in  the  wave-induced  alongshore  momentum 
flux  and  vertical  changes  in  the  wave-induced  Reynolds  stress  are  balanced  by  vertical 
changes  in  alongshore  turbulent  shear  stress. 

Applying  a  first  order  turbulence  closure,  the  alongshore  turbulent  shear  stress  can 
be  defined 


dV 
dz 


\(z)-  =   -pwv  =  p/;  —  (3) 


where  jj    is  the  turbulent  eddy  viscosity,  which  is  assumed  to  be  uniform  with  depth. 
Substituting  (3)  into  (2) 

dty(z)    _   dpvii        dpwv  . 

dz  dx  dz 


The  cross-shore  gradient  of  vii  is  constant  in  the  alongshore  direction  for  straight  and 
parallel  contours.  Within  the  surf  zone,  the  shallow  water  approximation  holds,  thus  vu  is 
assumed  independent  of  depth.  The  wave-induced  Reynolds  stress  term  (vvv)  can  arise  from 
sloping  bottom  effects  as  well  as  from  wave  amplitude  gradient  effects  and  it  would  have  a 
non-zero  contribution  even  for  linear  wave  theory  (Deigaard  and  Fredsoe,  1989). 
Nevertheless,  Rivero  and  Arcilla  (1995),  showed  that  vvv  is  a  linear  function  of  depth,  and 
hence  its  vertical  gradient  is  a  constant.  Consequently,  the  right  hand  side  of  (4)  is 
independent  of  depth  and  the  shear  stress  profile  is  determined  by  integrating  (4)  over  depth, 
to  give 


xy(z)  =  [xy{0)-Ty(-h)]  j.  (5) 


Linearly  varying  shear  stress  occurs  in  flows  driven  by  uniform  hydrostatic  pressure 
gradients,  such  as  in  steady  open  channel  flows,  which  are  well  described  by  a  logarithmic 
velocity  profile.  Therefore  it  is  hypothesized  that  a  steady,  uniform,  turbulent  boundary  layer 
flow  over  a  rough  surface  in  the  alongshore  direction  can  be  described  by  a  logarithmic 
profile: 


v  7    +    h 

V(z)  =  —  In  (i 1)  (6) 

K  Z„ 


where  z  is  positive  upwards  from  the  surface,  h  is  the  mean  water  depth,  k  is  the  Von 
Karman  constant  (0.4),  v,  is  the  alongshore  shear  stress  velocity  and  z0  is  the  physical 
roughness  height,  determined  by  bottom  topography  and  sediment  grain  size.  When  waves 
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are  present,  nonlinear  interactions  between  waves  and  currents  within  the  bottom  boundary 
layer  increase  the  bottom  shear  stress.  Following  Grant  and  Madsen  (1979),  this  additional 
stress  can  be  modeled  by  an  apparent  roughness  height  za ,  that  is  analogous  to,  but  larger 
than  z  . 


The  mean  bottom  shear  stress  (  xv(-/z))  is  related  to  the  longshore  shear  stress 
velocity  ( v  J  through 


Ty(z=-h)  =  pv/  (7) 

In  addition,  a  bed  shear  stress  coefficient,  C,,  can  be  calculated  assuming  a  quadratic  bed 
shear  stress  relationship 


xy(-h)  =  pCf  (w2+  v2)1/2v  (8) 


and  combining  with  (7),  gives 


2 


Cf  =  =====  (9) 

(u2  +  v2)1/2v 


which  includes  contributions  from  both  steady  (U,V)  and  non-steady  (  u,v,u,v)  velocity 

components. 

Surface  Layer 

The  surface  layer  is  governed  by  waves  and  wind.  Wave  effects  are  manifested  in 
three  ways:  1)  an  undulating  boundary  is  imposed  on  a  mean  current,  which  is  being 
measured  in  an  Eulerian  frame;  2)  the  mass  transport  velocity  in  the  crest-trough  region  due 
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to  obliquely  incident  waves  contributes  to  the  longshore  flow;  3)  modifications  of  the  profile 
associated  with  wave-breaking  generated  turbulence  and  the  eventual  contribution  of  wave 
rollers  to  the  mass  transport. 

The  shape  of  the  mean  current  profile  in  this  layer  is  determined  mostly  by  the  fact 
that  the  current  meter  is  intermittently  in  and  out  of  the  water  in  the  wave  crest-trough  region, 
e.g.,  for  a  linear  wave  the  current  meter  is  out  of  the  water  half  the  time  at  the  mean  sea  level 
(MSL  hereafter),  and  the  time-averaged  current  is  only  50  percent  of  the  expected  value  from 
the  logarithmic  profile.  To  account  for  this,  the  surface  elevation  probability  distribution 
function  {pdf  hereafter)  is  applied  to  the  expected  mean  current  profile  in  the  absence  of 
waves.  The  percent  of  time  the  current  meter  is  in  the  water  is  given  by  1-P(r|),  where  P(t|) 
is  the  cumulative  surface  elevation  pdf.  In  an  Eulerian  frame  of  reference,  the  modified  mean 
current  within  the  surface  layer  is  given  by 

V{z)  =  U-P(r\)]V'(z)  (10) 

where  V'(z)  represents  the  logarithmic  profile  in  the  absence  of  waves  (6).  For  moderate 
wave  conditions  in  deep  water,  the  surface  elevation  pdf  is  well  described  by  the  Gaussian 
pdf.  As  will  be  seen,  the  measured  pdf  in  the  surf  zone  are  slightly  positively  skewed  from 
the  Gaussian  distribution. 

The  mass  transport  velocity  associated  with  obliquely  incident  waves  can  contribute 
to  the  longshore  flow  in  the  upper  layer.  Assuming  irrotational  flow,  the  mass  flux  (transport) 
in  the  direction  of  wave  transport,  M,  can  be  evaluated  by  considering  separately  two  regions 
in  an  Eulerian  frame  of  reference  (Philips,  1977).  In  the  first  region,  from  the  bottom  to  the 
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MSL,  the  contribution  is  zero  for  irrotational  flow.  Within  the  second  region,  from  the  MSL 
to  the  water  surface  (rj),  u  is  not  defined  for  linear  wave  theory.  A  Taylor  series  expansion 
about  z=0  is  used  to  extend  defined  values  of  w(0)  at  the  surface,  giving  a  second  order 
approximation  to  M.  The  mass  transport  for  a  single  wave  is  interpreted  to  this  order,  in  an 
Eulerian  reference  frame,  as  due  to  a  uniform  velocity  confined  to  the  crest-trough  region. 
Applying  linear  wave  theory,  a  mass  transport  velocity  is  defined  by 

ti  M 

U  -  ^M)  (11) 

where  A  is  the  wave  amplitude. 

For  random  waves,  the  wave  amplitudes  can  be  described  by  the  Rayleigh  distribution 
to  a  first  approximation,  even  within  the  surf  zone  (Thornton  and  Guza,  1983).  The  only 
waves  that  contribute  at  any  elevation  z,  will  have  an  amplitude  A  >  Izl.  Assuming 
directionally  narrow  banded  waves,  the  ensemble  averaged  mass  transport  velocity  profile 
in  the  direction  of  wave  travel  is  obtained  by  applying  the  wave  amplitude  probability  density 
function 

<U(z)>  =   f"    U(A)p(A)dA  (12) 

where  <  >  represents  ensemble  averaging.  The  Rayleigh  probability  density  function  is  given 
by 

24    2 

8A        "(^_) 
p(A)  =  -^-  e      ™  where  0<A<h  (13) 
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and  H      is  the  rms  wave  height.  Substituting  (13)  back  into  (12),  applying  linear  wave 
theory,  and  performing  the  integration,  gives 


-(4^)2 


<U(Z)>  =  -    Hrms    [—  e    H™     +  &  ^c(^-)]  /or  lzl<A      (14) 


8  tanh  kh  H  2  H 

rms  rms 


where  erfc(x)  is  the  complementary  error  function,  u)  is  the  wave  radial  frequency  and  k  is 
the  radial  wavenumber.  The  alongshore  component  of  the  mass  transport  velocity  is  defined 
by 

<V(z)>  =  <U(z)>sind  (15) 

where  0  is  the  mean  incident  wave  angle  with  respect  to  the  beach  normal,  for  the  assumed 
directionally  narrow  banded  waves. 

The  mean  longshore  current  within  the  surface  layer,  as  observed  in  an  Eulerian 
reference  frame,  is  modeled  by  adding  the  corrected  logarithmic  profile  for  the  cumulative 
surface  elevation/?^  (10)  with  the  alongshore  component  of  the  mass  transport  velocity  (15) 

Vsfc(z)  =  [\-P(r\)]VJ(z)  *  <U(z)>smQ  (16) 

DUCK  94  EXPERIMENT 

The  measurements  described  here  are  part  of  the  comprehensive  nearshore  DUCK94 
experiment  conducted  during  October  1994  at  the  U.S.  Army  Corps  of  Engineers  Field 
Research  Facility  (FRF),  Duck,  North  Carolina.  The  FRF  is  located  on  the  Outer  Banks,  a 
barrier  island  formation  with  no  major  coastal  structures  to  obstruct  nearshore  flows.  The 
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beach  is  a  two-bar  system  with  a  dynamic  inner  bar  (30-120  m  offshore)  and  a  secondary  bar 
with  lower  amplitude  (300-400  m  offshore).  The  mean  foreshore  slope  of  the  beach  is 
approximately  0.08  (1:12)  and  the  slope  offshore  of  the  bars  is  approximately  0.006  ( 1 :  170) 
(Lippmann  etal,  1993).  The  mean  tidal  range  is  1.0  m.  Sediments  within  the  surf  zone  are 
well  sorted  with  a  mean  grain  size  of  0.2  mm.  Sediments  on  the  foreshore  are  poorly  sorted 
with  larger  mean  grain  size  (>0.4  mm). 

The  data  analysis  presented  is  for  the  October  phase  of  the  DUCK94  experiment. 
The  weather  during  October  was  climatologically  characterized  by  three  distinct  phases: 
weak  currents  and  winds  from  north  (4-9  Oct.),  relatively  strong  currents  from  north  (0.6- 
1 .0  m/s)  caused  by  a  storm  with  predominant  winds  and  waves  from  north  (10-17  Oct.),  and 
variable  currents  and  winds  from  north/south  (17-21  Oct.).  For  the  first  phase  (4-9  Oct.),  the 
currents  were  weak  and  barely  above  sensor  accuracy.  During  the  last  phase  (17-21  Oct.), 
currents  and  winds  were  highly  variable  and  a  hole  developed  in  the  bar  associated  with  a  rip 
current  system  (observed  with  dye),  such  that  the  basic  assumptions  of  steady  state  and 
straight  and  parallel  contours  are  violated.  Within  the  second  phase  (10-17  Oct.), 
observations  were  limited  due  to  the  Coastal  Research  Amphibious  Buggy  (CRAB)  not  being 
able  to  operate  during  the  height  of  the  storm  on  the  15th  and  being  limited  to  the  area  inside 
the  bar  on  the  13th,  14th,  16th,  and  17th  due  to  large  waves.  Therefore,  the  data  selected  for 
analysis  are  from  10-12  Oct.  (Fig.  1),  during  the  strong  longshore  currents  period  when 
observations  spanned  the  entire  surf  zone  and  conditions  approximate  the  assumptions  of 
steady  state  and  straight  and  parallel  contours.  Within  this  period  a  logarithmic  profile  is  well 
defined. 
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A  specially  designed  sled  was  used  as  a  platform  to  mount  instruments  (Fig.  2).  The 
sled  is  constructed  of  a  3  x  4  m  six-inch  aluminum-pipe  frame  with  two  5  m  length,  20  cm 
diameter  pipe  runners.  This  low-profile  structure  was  stabilized  by  180  Kg  of  lead  weight 
plus  approximately  450  Kg  of  sand  inside  the  runners.  In  addition,  there  were  four  fins  (45 
cm  wide)  extending  60  cm  into  the  sand  to  insure  that  the  sled  did  not  move  while  on  station. 

Currents  were  measured  using  a  vertical  stack  of  8  Marsh-McBimey  two-component 
electromagnetic  current  meters  (ems,  hereafter)  with  2.5  cm  diameter  spherical  probes 
mounted  along  a  2.5  m  mast  (Fig.  2).  The  ems  elevations  above  the  bed  were  23, 42,  68,  101 , 
147,  179,  224  and  257  cm.  The  ems  were  displaced  at  least  1  m  from  the  sled,  and  the  sled 
was  oriented  such  that  the  vertical  stack  of  ems  was  on  the  "up-current"  side  of  the  sled  to 
avoid  flow  contamination  by  the  sled  structure  during  observations.  The  ems  were  pre-  and 
post-calibrated  in  a  tow  tank  at  the  Naval  Postgraduate  School  with  an  agreement  of  1 .9 
percent  in  gain.  An  in  situ  determination  for  the  offset  is  used,  which  was  obtained  by 
reversing  the  orientation  of  the  ems  on  a  very  slow  longshore  current  day  (8  Oct.)  by  turning 
the  sled  around  and  returning  it  to  the  same  location  (within  1  m)  within  one  hour.  The  in  situ 
determined  offsets  were  within  1  cm/s. 

The  sled  orientation  was  determined  using  a  digital  compass  mounted  on  the  sled 
with  accuracy  0(1  degree).  Measured  two-component  velocities  were  reduced  to  a  shore 
normal  right-handed  coordinate  system  (positive  offshore  and  to  the  south)  by  using  compass 
data  and  adding  at  each  sled  position  any  deviation  of  the  contour  line,  as  measured  by  the 
CRAB,  from  a  shore  parallel  direction.  Velocity  errors  associated  with  the  rotation  of  the 
coordinate  system  were  smaller  than  the  determined  offsets,  and  therefore  were  neglected. 
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Waves  and  mean  water  level  were  measured  using  an  array  of  five  pressure  sensors 
configured  in  a  3  m  square  with  sensors  at  each  corner  and  one  in  the  center. 

The  data  were  digitally  encoded  on  the  sled  to  14  bit  precision  at  36  samples/second 
and  transmitted  to  shore  via  a  fibre-optic  cable  where  signals  were  monitored  and  recorded. 
Short  cables  from  the  sensors  to  the  data  acquisition  system  on  the  sled  (<  7  m)  resulted  in 
low  noise  ems  and  pressure  sensor  signals.  An  armored  cable,  married  to  the  sled  chain 
tether,  provided  power  and  controller  signals  for  the  instruments  via  two  conductors  and 
returned  the  digitized  signals  and  video  via  a  fibre-optic  line. 

The  sled  was  towed  to  the  farthest  offshore  location  for  the  first  run  (approximately 
160  m  from  the  shoreline)  by  the  1 1  m  high,  motorized,  three-wheel  CRAB.  A  four-wheel 
drive  forklift  pulled  the  sled  shoreward  10  to  30  meters  for  subsequent  runs  (each  run  was 
nominally  one  hour).  Five  to  eight  runs  were  made  across  a  transect  each  day.  The  data  were 
acquired  during  daylight  to  early  night,  which  happened  to  span  the  high  tide  during  this 
period. 

The  morphology  of  the  bottom  (bathymetry)  was  measured  at  various  scales  from  the 
CRAB.  Large-scale  variations  of  bathymetry  were  obtained  by  using  an  autotracking  laser 
ranging  system  to  measure  the  CRAB  position  approximately  every  meter  with  a  vertical 
accuracy  of  less  than  3  cm  rms.  Small-scale  vertical  bottom  variations  relative  to  the  CRAB, 
including  ripples  and  megaripples,  were  measured  with  a  1  MHZ  sonic  altimeter  mounted 
on  the  CRAB,  70  cm  from  the  bed.  The  altimeter  has  a  3.4  degree  beam  width  which 
translates  into  an  approximate  4  cm  footprint  and  a  nominal  sampling  rate  of  25  Hz,  which 
resulted  in  a  sample  spacing  of  2-4  cm  (dependent  on  CRAB  speed)  with  mm  vertical 
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resolution  and  accuracy  less  than  2  cm  (Gallagher  et  al.,  1996).  The  decrease  in  accuracy 
relative  to  resolution  is  due  to  the  changing  reflective  surface  owing  to  the  bed  dialating  or 
sediment  transported  along  the  bed  as  waves  pass  overhead.  The  CRAB  survey  and  altimeter 
measurements  were  combined  to  obtain  a  high  resolution  description  of  the  bottom 
(Thornton  et  al.,  1997).  Contour  plots  of  the  bathymetry  for  the  days  selected  for  analysis 
show  the  alongshore  contours  in  the  vicinity  of  the  sled  measurements  to  be  essentially 
straight  and  parallel  (Fig.  3;  1 1  Oct.  is  typical).  Bathymetry  for  the  three  selected  days 
(Figs.  4-6,  upper  panel)  show  a  pronounced  bar  progressively  moving  offshore  and 
significant  small-scale  morphology  in  the  trough.  Areal  variations  were  determined  using  a 
500  KHz  side-scan  sonar  also  mounted  on  the  CRAB. 

Meteorological  information  of  wind,  air  temperature,  atmospheric  pressure  and  sea 
surface  temperature  were  recorded  simultaneously  at  the  seaward  end  of  the  600  m  long  FRF 
pier  and  atop  the  FRF  building  in  front  of  the  pier. 

DATA  RESULTS 

The  data  are  qualitatively  sorted  by  location  into  the  two  regions  of  over  the  bar  and 
in  the  trough.  This  sorting  allows  a  better  identification  of  the  possible  correlations  among 
variables,  as  wave-breaking,  which  is  a  major  controlling  factor  within  the  surf  zone, 
significantly  changes  for  these  regions. 

Since  the  logarithmic  velocity  profile  hypothesis  is  a  bottom  boundary  layer  concept, 
information  from  ems  near  the  surface  influenced  by  the  effects  of  waves  and  wind  as  well 
as  from  coming  in  and  out  of  the  water  are  not  included.  A  criterion  is  established  such  that 
only  data  from  ems  below  (MSL-  H     )  are  considered  to  define  the  logarithmic  profile  of 
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the  bottom  boundary  layer.  This  criterion  assures  that  the  ems  used  in  the  analysis  came  out 
of  the  water  less  than  0.25  percent  of  the  time  based  on  a  Gaussian  distribution,  which  is 
conservative  for  the  measured  positively  skewed  distributions. 

The  rms  wave  height  is  approximated  by  Hrms  =  ySo2  where  a2  is  the  variance 
calculated  from  the  surface  elevation  time  series.  Surface  elevation  was  calculated  by  Fourier 
transforming  a  one-hour  pressure  record,  applying  a  linear  wave  theory  transfer  function  to 
the  complex  Fourier  amplitudes  in  the  frequency  domain,  and  inverse  transforming  to  obtain 
the  surface  elevation  time  series  (Thornton  and  Guza,  1982). 

All  22  vertical  profiles  of  longshore  currents  obtained  during  these  three  days  are 
analyzed.  The  profiles  are  based  on  the  measurements  by  four  to  seven  ems  over  the  vertical. 
The  em  closest  to  the  sea  bed  was  not  used  because  of  malfunction.  Mean  alongshore 
velocities  are  1  -hour  averaged  data  with  the  exception  of  three  runs,  run  7  on  1 1  Oct.  and 
runs  6  and  7  on  12  Oct.,  which  are  40-minutes  averaged  data.  In  examining  the  bottom 
boundary  and  surface  layers,  the  ems  data  are  treated  separately. 

For  the  bottom  boundary  layer,  logarithmic  profiles  are  fit  to  the  data  based  on  a 
linear-regression  least-squares  method.  The  value  of  za  is  calculated  from  the  z  intercept 
of  the  linear-regression  on  a  semi-log  plot  of  z  versus  V(z) ,  and  the  shear  stress  velocity  v, 
is  calculated  from  the  slope. 

Cf  values  are  calculated  using  measured  velocities  (u,v)  in  (9)  time-averaged  over 
the  record  length,  with  vt  determined  by  least-squares  fit.  Cf  values  are  dependent  on  the 
elevation  of  the  measured  u,v  values.  Typically,  Cf  is  calculated  using  measured  u,v  values 
referenced  to  one  meter  above  the  bed.  A  sensitivity  analysis  on  the  variation  of  Cf  to  the 
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elevation  of  the  selected  gage  was  performed,  showing  that  the  variation  of  C,  was 
dominated  by  the  time-averaged  alongshore  velocity  (V).  This  term  was  not  only  the  largest 
in  the  denominator  of  (9),  but  also  varied  the  most  with  depth.  Therefore,  to  minimize  the 

depth-dependence  of  C,    the  depth-averaged  velocity  calculated  from  the  logarithmic 

v»  h       za 

velocity  profile  (  Vd  =  —  [In ( — )  +  —  -  1  ])  is  used  to  specify  the  dominant  time-averaged 
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alongshore  velocity.  The  em  located  at  the  elevation  of  1  m  above  the  bed  is  used  to  measure 


the  smaller  terms  (total  cross-shore  velocity  u  and  the  non-steady  component  of  the 
alongshore  velocity  v  +  v  =  v  -  V),  which  have  weak  vertical  variation  outside  the  bottom 
boundary  and  surface  layers.  The  total  alongshore  velocity  is  obtained  by  adding  Vd  to  the 
measured  non-steady  component  v  +  v.  The  use  of  this  method  resulted  in  a  mean  variation 
of  C,  with  depth  of  only  7  percent  with  a  standard  deviation  of  6  percent  and  a  maximum 
variation  of  24  percent,  provided  the  gage  used  to  measure  the  cross-shore  velocity  and  non- 
steady  component  of  the  alongshore  velocity  was  not  in  the  surface  layer  (z  >  MSL  -  H^J. 
Calculated  Cf ,  za  and  v„  values  are  listed  in  Table  1. 

Error  estimates  of  za  and  v,  based  on  the  linear  regression  correlation  coefficients 
(C)  are  calculated  using  (Gross  and  Nowell,  1983  and  Cacchione  et  al,  1987) 

V.    ±   V.    V2,l-a/2)    (^-)m  (17) 

n-2 
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where  n  is  the  number  of  ems  used  for  the  regression,  and  t._2  j         is  the  Student's  t 
distribution  for  ( 1  -a) confidence  interval  with  (n  -2) degrees  of  freedom. 

The  uncertainties  of  Cf  values  can  be  determined  from  the  error  estimates  of  v,  and 

the     time     and     depth-averaged     velocity     calculated     from     the     measurements 

o     
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Vm  =  —         [(v  +u  )    v]     dz   (mean  velocity,  hereafter) 
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where  v,     is   the  relative  error  of  the  shear  stress  velocity  at  a  given  confidence  level 

err 

calculated  from  (17)  and  Vm     is  the  relative  error  of  the  mean  velocity  determined  by 
Vm     =  — ,  assuming  that  the  mean  velocity  has  a  constant  absolute  error  equal  to  the 

m 

in  situ  determined  offset  (1  cm/s).  Error  estimates  for  Cf,z,  and  vt  are  listed  in  Table  2, 
for  a  95  percent  confidence  interval. 

The  measured  cumulative  surface  elevation  pdf  is  used  to  describe  the  surface  layer, 
and  compares  reasonably  well  with  the  Gaussian  cumulative  pdf  (Fig.  7).  The  Gaussian 
distribution  has  zero  skewness  and  kurtosis  equal  to  three.  The  measured  skewness  values 
ranged  from  0.26-0.57  and  kurtosis  values  ranged  from  2.7-3.5,  indicative  of  weakly 
nonlinear  waves. 

In  the  surface  layer,  the  upper  ems  sometimes  came  in  and  out  of  water,  which 
causes  noise  in  the  ems  outputs.  To  eliminate  this  noise,  the  current  velocities  were  set  to 
zero  when  the  ems  were  within  5  cm  of  the  surface  as  determined  from  the  surface  elevation 
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time  series. 

The  combined  surface  and  bottom  boundary  layer  solution  is  the  logarithmic  profile 
given  by  (6)  from  za  up  to  (MSL  -  Brms),  and  the  modified  logarithmic  profile  by  the 
observed  cumulative  surface  elevation  pdf  plus  the  alongshore  component  of  the  mass 
transport  velocity  given  by  (16)  above  this  level.  As  a  typical  example,  the  profile  for  the 
seventh  run  of  10  Oct.  is  enlarged  in  Fig.  8,  where  the  solid  line  is  the  logarithmic  profile  and 
the  dashed  line  the  modified  profile  within  the  surface  layer. 

The  observed  and  model  predicted  velocity  profiles  at  successive  offshore  positions 
(runs)  that  the  sled  occupied  during  a  transect  are  shown  in  the  upper  panels  of  Figs.  4-6.  The 
largest  deviation  occurs  for  the  fourth  run  of  1 1  Oct.,  which  generates  an  outlier  in  the 
parameter  estimates  if  the  data  from  the  em  at  an  elevation  of  42  cm  above  the  sea  bed  (Fig. 
5)  is  included  in  the  linear  regression.  A  comparative  analysis  for  all  ems  for  this  run  showed 
the  energy  density  spectra  to  be  nearly  uniform  (average  variance  of  0.072  m  2and  standard 
deviation  of  0.013  m2)  with  exception  of  this  em  (variance  of  0.021  ra2).  Thus,  data  from 
this  em  is  disregarded  and  a  new  regression  is  calculated  with  the  remaining  six  ems.  This 
procedure  eliminates  the  outlier. 

DISCUSSION 

Three  days  of  the  DUCK94  experiment  are  examined  when  strong  longshore  currents 
occurred.  Mean  longshore  current  profiles  obtained  using  four  to  seven  ems  spaced  from  42 
to  257  cm  above  the  sea  bottom  are  used  to  test  the  validity  of  the  logarithmic  profile 
hypothesis  within  the  energetic  surf  zone  region. 
Bottom  Boundary  Layer 
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A  high  correlation  coefficient  for  the  linear  regression,  is  commonly  accepted  as  an 
indicator  of  the  validity  of  the  logarithmic  approach  (Grant  et  ai,  1984;  Gross  et  al.,  1994; 
Li,  1994).  The  linear  correlation  coefficients  for  all  profiles  ranged  from  0.95  to  0.99 
(Table  2),  with  an  average  value  of  0.98,  and  the  largest  deviations  occurring  over  the  bar, 
where  wave  breaking  was  strong.  Anomalous  high  correlation  coefficients  can  be  obtained 
for  profiles  calculated  using  a  small  number  of  ems  for  the  linear  regression,  such  as 
occurred  for  stations  very  close  to  the  shoreline  due  to  shallow  water.  Thus,  a  high 
correlation  coefficient,  although  necessary,  is  not  sufficient  to  validate  the  logarithmic 
profile  approach. 

Other  measures  of  how  well  the  logarithmic  model  describes  the  mean  alongshore 
current  profile  within  the  bottom  boundary  layer  are  the  uncertainties  of  the  calculated  shear 
stress  velocity  v„  (17)  and  apparent  roughness  height  za  (18).  These  uncertainties  reflect 
both  the  correlation  coefficients  and  the  number  of  ems  used  in  the  regression.  An  example 
is  comparing  the  fourth  run  of  1 1  Oct.  with  the  seventh  run  of  10  Oct.  (Table  2) .  The  latter 
has  a  higher  correlation  coefficient  (0.991),  but  only  4  ems  were  used  in  the  linear  regression; 
thus  uncertainties  in  the  calculated  v^  and  z  values  are  larger  for  this  run,  even  though  the 
former  has  a  lower  correlation  coefficient  (0.987).  Therefore,  the  uncertainties  are  not  biased 
by  an  anomalous  high  correlation  coefficient  due  to  a  small  number  of  ems  used  in  the 
regression.  The  degenerate  case  is  when  only  two  ems  are  used  for  the  linear  regression 
(eighth  run  of  1 1  Oct.),  resulting  in  a  correlation  coefficient  of  one,  which  does  not  allow 
the  calculation  of  the  uncertainties  for  vt  and  za.  Consequently,  data  from  this  run  are 
disregarded. 
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There  is  some  uncertainty  of  the  exact  distance  of  the  ems  from  the  bed,  because  the 
sled  runners  sank  into  the  sand  an  unknown  amount  depending  on  the  bearing  capacity  of  the 
bed  (theory  and  divers  suggest  3-8  cm)  and  owing  to  the  undulating  bottom,  particularly 
over  mega-ripples.  To  test  the  sensitivity  of  the  logarithmic  profile  due  to  the  uncertainty  in 
elevation,  the  elevation  of  the  ems  was  shifted  ±5  cm  in  steps  of  1  cm  and  the  linear 
regression  least-squares-fit  recalculated  until  the  best  correlation  between  the  data  and  the 
logarithmic  model  predictions  was  obtained  for  each  run  (Grant  et  al,  1984). 

As  the  bathymetry  is  not  uniform  over  the  transect  covered  by  the  sled,  no  unique 
zero-shift  value  should  be  expected  to  optimize  the  correlation  coefficient  for  all  runs. 
Another  effect  on  the  error  calculation  to  be  considered  with  an  ad  hoc  decrease  (increase) 
in  the  elevation  of  the  ems  is  that  an  em  previously  located  above  (below)  the  model  cut-off 
level  (MSL  -  H^)  can  be  relocated  to  a  new  elevation  below  (above)  this  level.  For 
example,  an  em  previously  neglected  by  the  model  is  now  included  in  the  linear  regression. 
The  zero-shifting  modifies  the  error  estimates  of  v,  ,  za,  and  C,  (17,  18,  and  19)  by 
changing  the  correlation  coefficients  and  the  number  of  ems  used  for  the  regression. 

With  the  uncertainties  involved,  the  zero-shifting  method  was  only  applied  whenever 
a  significant  reduction  (greater  than  10  percent)  was  obtained  in  the  shear  stress  velocity 
error  estimate.  Applying  this  criteria,  only  three  runs  were  shifted  (by  the  same  -  5  cm): 
seventh  runs  of  10  and  11  Oct.,  and  third  run  of  12  Oct.  The  sensitivity  of  the  ems  to 
elevation  from  the  bed  as  inferred  from  the  zero-shifting  method  for  these  three  runs  is 
shown  in  Table  3.  The  general  trend  being  that  both  the  shear  stress  velocity  and  apparent 
roughness  height  errors  decreased  with  increased  negative  shift,  although  the  improvement 
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is  only  significant  for  the  shear  stress  velocity.  The  values  listed  in  Tables  1  and  2  for  these 
runs  were  calculated  after  applying  the  zero-shift. 

The  observed  and  model  predicted  velocity  profiles  at  successive  offshore  positions 
(runs)  that  the  sled  occupied  during  a  transect  are  shown  in  the  upper  panels  of  Figs.  4-6.  The 
data  agree  well  with  the  model  indicating  that  the  contribution  from  momentum  mixing 


arising  from  interactions  between  cross-shore  and  longshore  currents  (  — - )  neglected  for 

dx 

simplicity  sake,  does  not  significantly  modify  the  vertical  profile  of  longshore  currents, 


suggesting  that  — is  independent  of  depth. 

dx 

The  largest  discrepancies  between  measured  and  modeled  profiles  occur  over  the  bar, 
where  wave  breaking  is  most  intense  (x-distance  between  220  and  240  m).  The  surface  layer 
during  wave  breaking  is  typified  by  an  intense  production  of  turbulence  that  eventually  is 
dissipated  in  the  shear  layer  at  the  lower  boundary  of  the  surface  roller  for  spilling  type 
breakers  and  more  intense  injection  of  turbulence  for  plunging  type  breakers.  The  increase 
of  turbulent  mixing  due  to  wave  breaking  produces  a  more  uniform  vertical  profile  of  the 
mean  longshore  current  within  the  bottom  boundary  layer  for  a  given  bottom  shear  stress, 
compared  with  profiles  in  the  absence  of  breaking.  Therefore,  larger  discrepancies  between 
observations  and  logarithmic  profile  predictions  (lower  correlation  coefficients)  would  be 
expected  for  increased  turbulent  mixing  caused  by  wave  breaking. 

To  test  this  hypothesis,  percentages  of  waves  breaking  were  determined  from  video 
recordings.  The  number  of  waves  breaking  is  determined  using  the  methods  of  Lippmann  and 
Holman  (1991),  while  the  total  number  of  waves  is  found  applying  the  zero-up-crossing 
method  to  the  surface  elevation  time  series  (see  Lippmann  and  Thornton,  1997,  for  details). 
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The  general  decrease  of  the  correlation  coefficient  between  observations  and  logarithmic 
profile  predictions  with  increasing    percentage  of  wave  breaking  (Fig.  9)  supports  this 
conclusion. 
Surface  Layer 

The  mean  current  profile  within  the  surface  layer  is  modified  in  the  presence  of  winds 
and  waves.  The  winds  act  directly  to  generate  longshore  currents  via  the  alongshore  surface 
wind  stress  component,  and  indirectly  by  generating  obliquely  incident  waves  which  then 
force  the  longshore  currents  via  changes  in  the  radiation  stress.  Whitford  and  Thornton 
(1993)  measured  the  various  terms  in  the  alongshore  momentum  balance  including  wind 
stress  and  wave  forces  during  the  SUPERDUCK  experiment  (15-18  October  1986)  at  the 
same  location  as  the  DUCK94  experiment.  They  found  an  average  wind  force  to  wave  force 
ratio  for  this  four-day  period  of  0.1 1  (range  0.02  -  0.33)  for  a  mean  alongshore  component 
of  wind  speed  of  3.6  m/s  (range  1 .9  -  5.3  m/s).  During  the  three-day  period  studied  here  the 
mean  alongshore  wind  speed  was  9. 1  m/s  (range  7.7  - 12.6  m/s)  resulting  in  an  average  wind 
force  to  wave  force  ratio  of  0.21  (0.10  -  0.50).  Despite  the  wave  force  dominance 
characterizing  wave-driven  longshore  currents,  the  wind  force  contribution  is  not  negligible. 

Winds  and  currents  are  approximately  from  the  same  direction  (Fig.  1 ),  thus  the  wind 
force  effect  is  to  increase  the  alongshore  currents.  As  the  logarithmic  profiles  are  fit  to  the 
data  based  on  a  linear-regression  least-squares  method,  the  direct  wind  force  effect  is  already 
included  in  the  logarithmic  profile.  The  indirect  wind  force  effect  of  wave  generation 
resulting  in  increased  wave  heights  (higher  H^  values)  is  also  included  in  the  alongshore 
component  of  mass  transport  velocity.  The  theoretical  exponential  decay  of  wind-induced 
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current  speed  with  depth  was  not  observed  within  the  surface  layer  due  to  the  dominant  effect 
of  the  ems  intermittently  being  out  of  the  water,  and  therefore  is  not  included  in  the  model. 

The  modification  of  the  mean  current  profile  within  the  surface  layer  is  modeled  here 
by  correcting  the  logarithmic  profile  predictions  for  measurements  in  an  Eulerian  frame  with 
an  undulating  boundary  and  adding  the  alongshore  component  of  the  mass  transport 
velocity  (16).  As  few  ems  were  located  above  the  (MSL  -  H^j)  level  for  each  run,  a  local 
correlation  coefficient  as  an  indicator  of  the  validity  of  this  approach  is  meaningless. 
Therefore,  the  validity  of  the  model  is  evaluated  by  comparing  model  predictions  with 
observations  for  the  ensemble  of  ems  located  within  the  surface  layer. 

The  linearity  of  (16)  allows  an  evaluation  of  the  contribution  of  each  term 
(intermittent-wetting  and  alongshore  component  of  mass  transport  velocity)  separately.  For 
the  intermittent-wetting  term  only,  a  plot  of  predicted  versus  observed  velocities  (Fig.  10) 
shows  good  agreement,  but  with  the  modified  logarithmic  profile  by  the  cumulative  surface 
elevation  pdf  slightly  under-predicting  the  observed  velocities.  The  rms  error  and  slope  of 
the  least-squares  linear  regression  between  predictions  and  observations  are  respectively  1 3 
percent  and  0.96,  giving  a  slope-error  of  -4  percent.  Next,  including  the  alongshore 
component  of  the  mass  transport  velocity  (Fig.  1 1 )  results  in  reduced  rms  and  slope  errors 
of  7  percent  and  0.5  percent.  Despite  the  small  magnitude  of  the  mass  transport  term  (on 
average  only  1 1  percent  of  the  intermittent-wetting  term),  its  addition  corrects  the  small 
under-prediction  of  using  only  (10)  and  improves  the  overall  agreement  with  the  data  by 
reducing  both  the  rms  and  slope  errors. 

Although  the  results  obtained  by  this  first  attempt  (to  the  knowledge  of  the  authors) 
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to  model  mean  longshore  currents  within  the  surface  layer  are  encouraging,  improvements 
can  be  achieved  by  including  contributions  from  at  least  two  physical  processes.  The  first, 
arises  from  the  additional  mass  transport  associated  with  the  presence  of  wave-breaking 
generated  rollers,  and  will  increase  the  model-predicted  velocity  within  the  surface  layer.  The 
second  is  associated  with  the  directional  spreading  of  random  waves  and  it  will  have  an 
opposite  effect,  reducing  the  model-predicted  velocity.  Preliminary  calculations  indicate  that 
these  terms  are  of  the  same  order  of  the  mass  transport  term,  and  thus  an  order  of  magnitude 
smaller  than  the  intermittent-wetting.  The  good  agreement  with  the  data  obtained  by  applying 
this  simple  model  is  attributed  to  the  dominance  of  the  intermittent-wetting  term  and 
compensating  effects  of  neglected  contributions  from  the  presence  of  wave  rollers  and 
directional  spreading. 
Bed  Shear  Stress  Coefficient 

The  bed  shear  stress  coefficient  (C,)  from  the  quadratic  friction  model  is  an 
important  parameter  in  both  nearshore  hydrodynamics  and  sediment  transport.  The  simplest 
formulation  for  longshore  currents  assumes  steady  state  wave  conditions  and  straight  and 
parallel  bottom  contours,  and  results  in  an  alongshore  balance  between  cross-shore  changes 
in  wave-induced  momentum  (radiation  stress,  S^)  with  the  bottom  shear  stress,  to  give 
(Thornton,  1970) 

1       dS 

=  ^Ff"a7  (20) 

pCJw|    ax 

where  \u\  is  the  magnitude  of  the  total  velocity  vector.  The  sediment  transport  formulation 
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of  Bailard  and  Inman  ( 1 98 1 ),  based  on  the  work  by  Bagnold  ( 1 966),  states  that  the  immersed- 
weight  sediment  transport  rate  is  proportional  to  the  local  rate  of  energy  dissipation 

it  -9Cf\ut\HKbt  +Ks)  (21) 

where  Kbt  and  Kst  are  dimensionless  time-varying  vectors  associated  with  the  bedload  and 
suspended  load  transport  rates.  It  is  noted  that  both  the  longshore  currents  and  sediment 
transport  rates  are  direct  functions  of  Cf.  Nevertheless,  the  physical  processes  governing 
these  phenomena  are  inherently  different.  The  longshore  currents  are  mainly  controlled  by 
form  drag  of  the  bed  forms  and  non-linear  interactions  between  waves  and  mean  current. 
Sediment  transport  is  mainly  related  to  skin  friction  due  to  sediment  grains  (Smith,  1977; 
Dyer,  1980),  although  the  presence  of  waves  increases  the  sediment  entrainment  rate  (Grant 
and  Madsen,  1979),  thus  increasing  the  transport  rate. 

The  bottom  shear  stress  coefficient  ( C, )  varied  by  an  order  of  magnitude  across  the 
surf  zone,  with  the  values  offshore  and  over  the  bar  in  the  order  of  10"3,  while  the  values  in 
the  trough  were  in  the  order  of  10"2.  An  attempt  was  made  to  find  empirical  relationships 

between  Cf  and  measured  physical  parameters  commonly  used  throughout  the  literature 

I  "J 
such  as  (ratio  of  near-bottom  wave  velocity  magnitude  from  linear  wave  theory  to 

mean  current  speed),  the  rms  bottom  roughness  (r)  and  percent  of  wave  breaking. 

\u,\ 
Surprisingly,  no  statistically  significant  correlation  was  found  between  Cf  and  . 

Cf  was  found  to  be  negatively  correlated  with  percentage  of  wave  breaking,  with  a 

linear  correlation  coefficient  of  -0.71,  which  is  statistically  significant  at  the  99  percent 

confidence  level  (Fig.  12).  This    is  contrary  to  the  theoretical  analyses  of  Fredsoe  and 
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Deigaard  (1992)  and  Church  and  Thornton  (1993).  An  effect  associated  with  wave  breaking 
in  the  surf  zone  is  the  generation  of  surface  rollers  that  increases  the  mass  transport  within 
the  surface  layer,  thus  increasing  the  cross-shore  return  flow,  or  undertow,  and  consequently 


the  average  total  velocity  (w2+  i^)1/2v.  This  effect  reduces  the  C,  values  calculated  using  (9) 
for  increasing  wave  breaking.  In  addition,  decreased  Cf  with  increased  wave  breaking  may 
be  physically  related  to  the  planing-off  of  wave-ripples  due  to  increased  near-bottom 
velocities  associated  with  wave  breaking  and  consequent  reduction  of  form  drag,  as  well  as 
to  stratification  caused  by  an  intense  vertical  gradient  of  suspended  sediment  near  the  sea 
bed.  This  stratification  can  diminish  the  turbulent  flow  intensity  in  the  bottom  boundary 
layer,  leading  to  reduced  bottom  stress  (Smith  and  McLean,  1977a;  and  Taylor  and  Dyer, 
1977). 

Cf  would  be  expected  to  be  related  to  bottom  roughness  due  to  enhanced  form  drag. 
The  bottom  roughness  is  examined  by  calculating  wavenumber  spectra  of  the  bed.  To 
calculate  wavenumber  spectra,  the  unevenly  spaced  data  from  the  combined  CRAB  surveys 
and  altimeter  measurements  are  linearly  interpolated  to  evenly  spaced  2  cm  increments  of  the 
cross-shore  distance.  The  small-scale  morphology  in  general  shows  large  cross-shore 
variation;  as  a  consequence,  the  condition  of  spatial  homogeneity  (stationarity)  required  for 
calculating  averaged  spectra  is  not  met.  Therefore,  continuous  bottom  roughness 
wavenumber  spectra  are  calculated  for  20  m  cross-shore  segments  at  increments  of  1  m 
across  the  surf  zone. 

Lowest  wave  numbers  are  filtered  by  subtracting  a  third-order  polynomial  best-fit 
curve  from  each  20  m  section.  A  10  percent  cosine-taper  data  window  is  applied  to  decrease 
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spectral  leakage.  The  spectra  are  summed  over  three  wavelength  bands  (0.2  -  1.0  m,  1.0  - 
3.33  m,  and  3.33  -  20  m)  plus  the  total  band  (0.2  -  20  m),  resulting  in  160,  28,  10,  and  198 
degrees  of  freedom  for  each  band  respectively.  The  wavelength  bands  chosen  are  based  on 
examination  of  individual  spectra.  The  spectra  were  generally  broad,  indicating  that  several 
ripple  wavelengths  coexisted  as  a  result  of  newly  formed  ripples  combined  possibly  with 
residual  ripples  from  the  past  to  form  a  complex  series  of  ripple  patterns,  plus  the  effects  of 
alignment  of  ripples  relative  to  the  cross-shore  measurement  axis.  The  rms  height  of  each 
band  is  calculated  as  the  square  root  of  the  variance  within  each  band.  Note  that  variances 
sum,  not  rms  heights,  such  that  the  rms  height  of  the  sum  of  the  three  bands  is  calculated 
from  the  square  root  of  the  sum  of  their  variances.  The  rms  bottom  roughness  (Figs.  4-6, 
lower  panel)  is  dominated  by  longer  wavelength  features  primarily  associated  with  mega- 
ripples,  with  a  small  contribution  from  the  superimposed  shorter  wavelength  ripples  0(1  m). 
The  general  trend  being  that  bottom  roughness  was  smoothest  offshore  and  over  the  bar 
where  wave  ripples  were  planed-off  due  to  higher  near-bottom  velocities,  with  increased 
roughness  within  the  trough  associated  with  mega-ripples  (see  Thornton,  et  ai,  1997,  for 
details).  The  rms  bottom  roughness  (r)  is  calculated  from  the  total  band  (0.2  -  20  m), 
and  thus  includes  contributions  from  both  large  and  small-scale  morphology.  Calculated 
values  are  listed  in  Table  1 . 

Roughness  measurements  were  made  in  the  cross-shore  direction  only.  To  relate  Cf, 
calculated  from  the  alongshore  bed  shear  stress,  to  roughness,  it  is  assumed  that  the  bed 
forms  are  quasi-isotropic,  which  is  not  unreasonable  for  mega-ripples  in  the  trough. 
However,  this  assumption  is  violated  for  long-crested  wave  ripples,  such  as  occurred  for  the 
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seventh  run  of  1 1  Oct.  when  large  wave  ripples  (observed  in  the  side-scan  sonar)  were 
oriented  parallel  to  the  flow  (alongshore  direction)  which  results  in  over-estimating  the  rms 
roughness  associated  with  longshore  currents  at  this  cross-shore  position.  Altimeter  data  for 
12  Oct.  is  not  available  from  210-260  m  due  to  the  high  false  return  levels  induced  by  strong 
scattering  of  the  acoustic  pulses  from  air-bubbles  entrained  by  waves  breaking  over  the  bar. 
Therefore,  the  CRAB  survey  data  is  used  as  the  profile  in  this  section,  and  hence  the  bottom 
is  not  as  highly  resolved  in  this  area,  compromising  roughness  estimates  for  the  third,  fourth, 
and  fifth  runs  of  12  Oct.  Thus,  roughness  data  from  these  four  runs,  as  well  data  from  the 
eighth  run  of  1 1  Oct.  (only  two  ems  below  the  MSL  -  HTms  level)  are  disregarded. 

If  data  from  runs  with  logarithmic  profile  correlation  coefficients  less  than  0.98  with 
resulting  relative  errors  for  C,  greater  than  50  percent  (Table  2)  are  also  disregarded,  only 
9  profiles  could  be  used,  making  the  ensemble  too  small  to  infer  any  statistically  reliable 
relationship  between  Cf  and  bottom  roughness.  Therefore,  it  was  decided  to  use  data  from 
the  remaining  17  profiles  with  the  associated  95  percent  confidence  error  estimates  to 
examine  the  relationship  between  Cf  and  r.  The  bed  shear  stress  coefficients  are  positively 
correlated  with  bottom  roughness  normalized  by  mean  water  depth  (Fig.  13),  with  a  linear 
correlation  coefficient  of  0.63,  which  is  statistically  significant  at  the  99  percent  confidence 
level. 

The  scatter  of  data  observed  in  Fig.  13  may  be  due  to  correlating  the  non-synoptic 
velocity  measurements  with  bottom  roughness  measured  once  in  the  morning,  prior  to  the 
positioning  of  the  sled  for  the  first  station.  Wave  forcing  quantified  by  deep  water  wave 
height  ( H  ),  period  of  peak  wave  frequency  ( T),  and  wave  direction  ( 6 )  changed  during  the 
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period  of  observation  (Fig.  1)  due  to  both  wind  and  tidal  variations.  Changes  in  wave  forcing 
have  a  direct  effect  in  the  measured  velocities  and  an  indirect  effect  in  bottom  roughness  due 
to  modification  of  wave  ripples  associated  with  variations  in  wave  height. 

For  steady  currents  (open  channel  flows)  C,  depends  only  on  bottom  roughness  when 
the  bed  is  hydraulically  rough,  and  several  empirical  relationships  (power-laws)  are  available 
throughout  the  literature.  One  of  the  most  widely  used  is  the  Manning-Strickler  equation 
(Sleath,  1984,  equation  5.86) 

f  k   - 

C   =  i™  =  0.015  (-^)3  (22) 

f        8  h 

where  fDW  is  the  Darcy  Weisbach  friction  coefficient  and  Jcs  is  the  bed  roughness  length  scale 
or  equivalent  Nikuradse  roughness  of  the  bed. 

In  our  definition  of  C,  (9),  contributions  from  both  steady  and  non-steady  velocity 
components  were  taken  into  account,  whereas  (22)  applies  to  steady  flow  only.  Nevertheless, 
during  the  period  being  analyzed  the  observed  strong  longshore  current  is  not  only  the 

dominant  term,  but  is  also  well  represented  by  a  logarithmic  profile  (6).  The  flow  can  be 

k  v, 
classified  as  hydraulically  rough  as  calculated  Reynolds  numbers  (— — ,  in  which  u  = 

molecular  kinematic  viscosity  coefficient)  were  greater  than  70  for  all  runs  (Schlichting, 

1979).  Thus,  if  measured  rms  bottom  roughness  is  representative  of  the  equivalent  Nikuradse 


roughness  of  the  bed  (k  ~r),  one  might  expect C,  to  be  related  to  (  — )   by  a  power  law 


similar  to  (22).  A  plot  of  predicted  C,  by  the  Manning-Strickler  equation  (22)  is  included  in 
Fig.  13,  which  shows  an  order  of  magnitude  agreement  between  the  data  and  predictions. 
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These  results  suggest  that  the  calculated  rms  bottom  roughness  (r)  is  representative  of  the 
bed  roughness  length  scale  or  equivalent  Nikuradse  roughness  of  the  bed  (k  ). 

The  structure  of  the  flow  inside  the  wave  boundary  layer  (WBL  hereafter)  is  out  of 
the  scope  of  this  paper  as  no  measurements  were  made  within  this  region.  Outside  the  WBL, 
za  is  the  proper  length  scale  and  not  z0  (Grant  and  Madsen,  1986).  Therefore,  the  equivalent 
Nikuradse  roughness  of  the  bed  (ks)  is  replaced  by  an  apparent  bed  roughness  length  scale 

( k  ).  The  flow  is  hydraulically  rough  during  the  period  being  analyzed,  and  thus  k   can  be 

k 
assumed  to  be  given  by  —  =  za  (Nikuradse,  1933).  Assuming  that  ka  is  representative  of 

the  bed  roughness  length  scale  when  waves  and  currents  are  present,  an  empirical  relation 

k 

(power  law)  between  C,  and  —  is  obtained  by  applying  a  linear  regression  least-squares  fit 

1  h 

to  the  data  to  give 

k    -L 
r   =  0.01 1(— )275  (23) 

f  h 


which  is  plotted  in  Fig.  14.  The  rms  error  between  measured  and  predicted  C,  by  (23)  is  18 
percent,  which  is  within  the  95  percent  confidence  band  of  the  measurements.  This  is  not 

surprising  as  the  coefficients  of  (23)  were  determined  by  fitting  the  data  in  a  least-squares 

k 

sense.  Nevertheless,  the  correlation  between  measured   C,  and   —  is  0.80  which  is 

f  h 

statistically  significant  at  the  99.5  percent  confidence  level  indicating  that  a  power  law 
relation  between  these  parameters  is  not  fortuitous. 

An  independent  check  of  (23)  is  obtained  by  including  the  results  from  velocity 
profiles  measured    by  Grant  et  al.  (1984)  on  the  northern  California  continental  shelf  in 
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90  m  depth  during  the  Coastal  Ocean  Dynamics  Experiment  (CODE)  using  four  vertically 


stacked  acoustic-travel-time  current  meters  (BASS)  mounted  on  tripods.  A  plot  of 

k 

measured  Cf  versus  —  for  the  CODE  data  (listed  in  table  2  of  their  paper)  is  included  in  Fig. 

J  h 

14.  The  agreement  between  C,  predicted  by  (23)  and  measured  for  the  CODE  data  set  is  even 
better  (rms  error  of  3.5  percent)  than  for  the  DUCK94  data  set  which  was  used  to  determine 
the  coefficients  of  (23). 

These  results  indicate  that  a  single  roughness  length  scale  can  be  used  to  characterize 
combined  flows  over  a  movable  bed,  which  is  a  basic  assumption  for  all  existing  theoretical 
models.  This  assumption  was  shown  to  be  valid  for  waves,  currents,  and  combined  flows 
over  a  fixed  rippled  bed  by  Mathisen  and  Madsen  (1996). 

Another  reason  for  the  scatter  of  data  observed  in  Fig.  1 3  arises  from  C,  values  being 
calculated  from  the  total  bottom  shear  stress,  which  has  stress  contributions  from  skin 
friction  related  to  sediment  grains,  waves-current  nonlinear  interactions  within  the  bottom 
boundary  layer  as  well  as  form  drag  related  to  bed  forms.  Smith  and  McLean  (1977b) 
linearly  partitioned  the  total  bed  shear  stress  into  skin  friction  and  form  drag,  and  found  good 
agreement  with  data  from  the  Columbia  River.  Nelson  and  Smith  (1989),  Wiberg  and  Nelson 
(1992),  and  Li  (1994)  subsequently  applied  the  linear  partition  concept  to  several  flume 
experiments  with  good  results.  Extending  this  concept  to  the  surf  zone  environment  requires 
including  an  additional  component  due  to  non-linear  interactions  between  waves  and 
currents  within  the  bottom  boundary  layer  to  the  total  bed  shear  stress.  Assuming  that  the 
linear  stress  partition  is  valid  within  the  surf  zone,  the  skin  friction  contribution  can  be 
removed  from  the  total  bottom  stress  and  a  new  bed  shear  stress  coefficient  Cj-  is  defined 
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2  "> 

V         +    V 
c  /    =  J  (24) 

(w2  +  v2)1/2v 


where  v„   and  v,    are  the  form  drag  and  waves-current  interactions  shear  stress  velocities. 

d  wc 

The  relationship  between  C,   and  the  total  bed  shear  stress  coefficient  (C,)  can  be 
determined  from  (9)  and  (24) 


Cj  =  [l-(— -)]  C.  (25) 

V  J 


where  v,  is  the  skin  friction  related  shear  stress  velocity. 

As  skin  friction  was  not  measured  during  the  DUCK94  experiment,  an  attempt  is 
made  to  isolate  its  contribution  from  the  total  bottom  shear  stress  by  applying  a  stress 
partitioning  model.  The  probabilistic  approach  used  to  quantify  bottom  roughness  does  not 
allow  the  adjustment  of  the  two  empirical  coefficients  CD  and  al  necessary  to  apply  the 
Smith  and  McLean  (1977b)  model.  Therefore,  the  empirical  relationships  obtained  by  Li 
(1994)  by  applying  linear  regression  to  the  mean  flow  laboratory  data  of  Kapdasli  and  Dyer 
(1986),  Paola  (1983) ,  as  well  as  his  own  data  obtained  from  mobile  sand  ripples,  are  used 
to  estimate  skin  friction  shear  stress  velocity  from  the  total  shear  stress  velocity  obtained 
from  the  logarithmic  profile 


v,  V  V 

-±  =  0.125(  — )+  0.373  for    —  <  2.3  (26) 

v  R  R  K     } 
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V.  V  V 

— ^  =  0.107(  — )  +0.266  for    —  >  2.3  (27) 

v  R  R  K     ' 


where  R  is  ripple  height,  that  is  assumed  here  to  be  equal  to  the  measured  rms  bottom 

roughness  (R  =  r)  .    For  — -  >  6.86,  the  predicted  v^    by  (27)  is  greater  than  the  total 

R 

measured  shear  stress  velocity,  and  it  is  assumed  that  v.  is  equal  to  99  percent  of  the  total 

s 

shear  stress  whenever  this  occurred.  Calculated  values  of  vt  are  listed  in  Table  1 . 

s 

The  recalculated  bed  shear  stress  coefficients  Cf  by  (25)  show  a  higher  linear 
correlation  coefficient  (0.76)  with  bottom  roughness  normalized  by  mean  water  depth  (Fig. 
15).  Theoretically  this  is  expected  as  the  removal  of  the  skin  friction  component  from  the 
total  stress  should  enhance  the  form  drag  contribution,  and  consequently  increase  the 
correlation  between  bed  shear  stress  and  bottom  roughness.  The  improved  correlation 
obtained  by  applying  Li  (1994)  empirical  relationships  indicates  that  the  linear  stress 
partition  concept  introduced  by  Smith  and  McLean  (1977b)  can  be  extended  to  the  surf  zone 
environment,  although  a  more  detailed  data  set  is  necessary  to  validate  these  expressions. 
Apparent  Roughness  Height 

The  change  of  mean  current  velocity  profile  due  to  the  superposition  of  waves  has 
been  theoretically  investigated  over  the  last  two  decades  and  several  models  have  been 
presented  (Lundgren,  1972;  Smith,  1977;  Bakker  and  Van  Doom,  1978;  Grant  and  Madsen, 
1979;  Fredsoe,  1984;  Christoffersen  and  Jonsson,  1985;  Myrhaug  and  Slaattelid,  1989; 
Sleath,  1991).  A  general  approach  common  to  all  these  works  is  splitting  the  mean  velocity 

37 


profile  in  an  inner  region  close  to  the  bed  affected  by  turbulence  produced  in  the  wave 
boundary  layer,  and  an  outer  region  above  it,  where  the  current  is  not  directly  affected  by 
wave-current  non-linear  interactions  so  that  the  logarithmic  profile  (6)  may  be  assumed  to 
apply.  The  increased  bed  shear  stress  associated  with  wave-current  interactions  is  modeled 

by  replacing  zowith  za'm  (6). 

z 
The  apparent  roughness  increase  —  is  theoretically  expected  to  depend  on  the 

\u,\  Z°  A, 

relative  current  strength  ,  the  relative  roughness  —  (ratio  of  the  near-bottom  wave 

V  k 

s 

semi-orbital  excursion  from  linear  wave  theory  to  equivalent  Nikuradse  roughness),  and  the 

z 
angle  between  wave  and  current  direction  ( <t>).  Several  empirical  relations  for  —  based  on 

these  parameters  have  been  obtained  from  laboratory  experiments.  Nevertheless,  the 

applicability  of  these  empirical  relations  to  field  data  has  not  been  verified,  as  previous  field 

experiments,  although  containing  reliable  data  on  za ,  lack  corresponding  measurements  of 

bottom  roughness  (Grant  et  al.,  1984;  Cacchione  et  al.,  1987;  Lambrakos  et  al.,  1988; 

Slaattelidefa/.,  1990). 

z 
Three  widely  applied  empirical  relationships  for  —  (Coffey  and  Nielsen,  1986; 

Zo 

Sleath,  1991 ;  and  Van  Rijn,  1993)  are  tested  using  DUCK94  data,  as  measurements  of  both  za 

and  bottom  roughness  were  made  during  this  experiment.  For  hydraulically  rough  flows,  the 

k 
physical  roughness  height  can  be  estimated  from  z    =  —  (Nikuradse,  1933),  and  once  again 

the  measured  rms  bottom  roughness  is  used  to  represent  the  equivalent  Nikuradse  roughness 

of  the  bed  (jfc    =  r). 

z 
No  statistically  significant  correlation  was  found  between  observed  —  and  the 

I  "J     Ad 
parameters  ,  — ,  and  (j).  Consequently,  none  of  these  empirical  relations  are  able  to 

s 
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accurately  predict  the  apparent  roughness  increase  for  this  data.  These  results,  together  with 

z 
the  scatter  of  available  laboratory  data  on  —  (Nielsen,  1992,  Fig.  1.5.13)  and  difficulties 

associated  with  bottom  roughness  and  velocity  measurements  within  the  WBL  in  the  field, 

indicate  that  much  remains  to  be  done  before  a  reliable  empirical  or  theoretical  relation  for 

the  apparent  roughness  height  increase  experienced  by  the  mean  current  profile  in  the 

presence  of  waves  is  obtained. 

CONCLUSIONS 

The  vertical  structure  of  mean  longshore  currents  on  a  barred  beach  is  well  described 
by  a  logarithmic  profile  for  the  three  strong  longshore  current  days  examined.  This 
hypothesis  works  better  in  the  trough  where  turbulent  bottom  boundary  layer  processes  are 
more  dominant  than  over  the  bar,  where  breaking-wave  induced  turbulence  generated  at  the 
surface  modifies  the  profile. 

The  modification  of  the  mean  longshore  current  profile  within  the  surface  layer  in  the 
presence  of  winds  and  waves  is  modeled  by  correcting  the  logarithmic  profile  predictions  for 
measurements  in  an  Eulerian  frame  with  an  undulating  boundary  for  all  ems  located  above 
the  (MSL  -  Hrms)  level.  The  addition  of  the  alongshore  component  of  the  mass  transport 
velocity  corrects  the  under-prediction  of  the  modified  logarithmic  profile  and  improves  the 
overall  agreement  with  the  data.  Therefore,  it  can  be  concluded  that  this  simple  model 
provides  a  first  order  approximation  that  is  sufficiently  accurate  to  predict  mean  longshore 
currents  within  the  surface  layer. 

The  data  indicates  that  wave  breaking  inside  the  surf  zone  decreases  the  bottom  shear 
stress  coefficient  which  is  contrary  to  the  theoretical  analysis  of  Fredsoe  and  Deigaard 

39 


(1992),  and  Church  and  Thornton  (1993).  This  may  be  physically  related  to  wave-breaking 
generated  surface  rollers  that  increases  the  cross-shore  return  flow  (undertow)  and 
consequently  the  average  total  velocity,  decreasing  Cf values  calculated  using  (9);  as  well  as 
to  the  planing-off  of  wave-ripples  due  to  increased  near-bottom  velocities  and  consequent 
reduction  of  form  drag,  and  to  stratification  caused  by  an  intense  vertical  gradient  of 
suspended  sediment  near  the  sea  bed. 

The  bed  shear  stress  coefficient  ( Cf )  varied  by  an  order  of  magnitude  across  the  surf 
zone  (0.0006-0.012),  with  the  values  offshore  and  over  the  bar  O(10"3),  while  the  values  in 
the  trough  were  O(10"2).  Thus,  longshore  current  and  sediment  transport  models  that  assume 
Cf  to  be  constant  or  mildly  changing  should  be  revised. 

C,  was  found  to  be  directly  proportional  to  bottom  roughness,  and  hence,  bottom 
roughness  is  an  important  parameter  to  characterize  the  bottom  boundary  layer.  The 
empirical  relationships  obtained  by  Li  (1994)  were  used  to  remove  the  skin  friction 
contribution  from  the  total  bottom  shear  stress.  The  improved  correlation  between  bottom 
shear  stress  and  bottom  roughness  obtained,  although  not  conclusive  to  validate  these  simple 
expressions,  indicates  that  the  linear  stress  partition  concept  introduced  by  Smith  and 
McLean  (1977b)  can  be  extended  to  the  surf  zone  environment. 

An  empirical  relation  between  Cf  and  apparent  roughness  length  scale  of  the  bed 

k 
normalized  by  water  depth  ( — )  was  obtained  by  applying  a  linear  regression  least-squares 

h 

fit  to  the  DUCK94  data.  Good  agreement  was  also  found  between  this  relation  and  CODE 
data  (Grant  et  ah,  1984),  indicating  that  a  single  roughness  length  scale  can  be  used  to 
characterize  combined  wave-current  flows  over  a  movable  bed. 
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Surprisingly,  no  statistically  significant  correlation  was  found  between  observed 

z                                     \u.\     A. 
apparent  roughness  increase  ( — )  and  the  parameters ,  — ,  and  4> .  Consequently,  none 

Zo  V         K 

Z 

of  the  three  most  used  empirical  relations  for —  (Coffey  and  Nielsen,  1986;  Sleath,  1991; 

Zo 

and  Van  Rijn,  1993)  were  able  to  accurately  predict  the  apparent  roughness  increase.  These 
results  indicate  that  despite  improvements  in  our  knowledge  of  turbulent  boundary  layers 
during  the  last  two  decades,  we  are  still  not  able  to  estimate  bottom  roughness  for  field 
applications,  which  prevents  the  use  of  an  empirical  relation  such  as  (23)  or  any  other 
theoretical  expression  to  predict  friction  factor  for  combined  waves  and  current  flows. 
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Figure  1.  Climatology  for  the  three  days  considered  (10,  1 1  and  12  Oct.).  Currents  were 
measured  in  the  middle  of  the  trough.  H  is  the  significant  wave  height,  T  is  the  period  of 
peak  frequency,  and  Angle  is  the  mean  wave  angle  relative  to  beach  normal  (0). 
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Figure  2.  Photograph  of  the  Sled  being  pulled  off  the  beach  by  the  CRAB  (upper  right) 
during  DUCK94  experiment.  Vertical  array  of  ems  are  mounted  along  Sled  mast  (center 
left). 


49 


Instruments 
:  Guxo  &  Cigar 
-  Herbers  4  Guzo 
: Hoy  &  Bowen 
:Hom«5  4  Geltonboum 
: Beach.  Sternberg  4  Holmon 
:  White 

: Hones  4  Vincent 
:FRF 4  Howd 
:  Thornton  4  Stonton 
:  Lippmonn 

:  Livingston  4  Posework 
:  Holland  4  Holmon 


SO 


150 


250 

Cross-shore,  m 


350 


450 


Figure  3.  Bathymetry  and  instrument  locations  during  DUCK94  experiment  (1 1  Oct.).  Sled 
transect  line  was  940  m  alongshore. 
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Figure  4.  Measured  (  *  )  and  predicted  (  — )  vertical  profiles  of  mean  longshore  currents 
superposed  on  bottom  profile  with  tide  elevation  indicated  by  (  o  )  and  measured  C,  values 
(upper  panel).  Variation  of  band-limited  rms  bottom  roughness  with  cross-shore  distance  for 
10  Oct  (lower  panel). 
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Figure  5.  Same  as  in  Figure  4  but  for  1 1  Oct.  The  anomalous  measurement  point  indicated 
by  +  in  the  fourth  run  (upper  panel)  is  not  included  in  the  linear  regression. 
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Figure  6.  Same  as  in  Figure  4  but  for  12  Oct. 
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Figure  7.  Comparison  between  measured  and  Gaussian  cumulative  surface  elevation  pdf  for 
all  the  runs  during  the  three  days  considered  (10,  11  and  12  Oct.).  Mean  skewness  and 
kurtosis  for  the  three  days  are  indicated. 
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Figure  8.    Vertical  profile  combining  bottom  boundary  and  surface  layer  models  for  the 
seventh  run  of  10  Oct.  The  upper  three  ems  sometimes  were  out  of  the  water. 
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Figure  9.  Correlation  coefficient  between  observations  and  logarithmic  profile  predictions 
versus  percentage  of  wave-breaking  calculated  from  video  data.  The  line  represents  a  linear 
regression  with  a  regression  coefficient  of  -0.94. 
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Figure  10.  Predicted  (equation  10)  versus  observed  longshore  velocities  within  the  surface 
layer.  The  solid  line  represents  perfect  agreement  (slope=l)  and  the  dashed  line  represents 
a  linear  regression  with  a  slope  of  0.96. 
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Figure  11.  Same  as  in  Figure  10  but  including  the  alongshore  component  of  mass  transport 
velocity  (equation  16).  The  slope  of  the  linear  regression  is  1.005. 
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Figure  12.  Bed  shear  stress  coefficient  (  C, )  versus  percentage  of  wave-breaking  calculated 
from  video  data.  The  line  represents  a  linear  regression  with  a  regression  coefficient  of  -0.7 1 . 
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Figure  13.  Bed  shear  stress  coefficient  (Cf  )  versus  rras  bottom  roughness  normalized  by 

mean  water  depth  ( — ).  The  dashed  line  is  the  Manning-Strickler  equation  (22)  and  the  solid 

h 
line  is  obtained  by  adjusting  the  coefficients  of  (22)  by  applying  a  linear-regression  least- 
squares  method  to  the  data. 
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Figure  14.  Bed  shear  stress  coefficient  (  Cf )  versus  apparent  bed  roughness  length  scale 

k  } 

normalized  by  mean  water  depth  ( — ).  The  line  represents  predicted  C,  by  equation  23. 
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Figure  15.  Bed  shear  stress  coefficient  without  skin  friction  contribution  (Cf)  versus  rms 

bottom  roughness  normalized  by  mean  water  depth  (  — ).  The  line  represents  a  linear 

h 

regression  with  a  regression  coefficient  of  0.76. 
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Table  1.  Logarithmic  Profile  Fitting  Results. 
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Table  2.   95  percent  confidence  bands  on    vm  ,  za  ,  and    C,    .  Run  8  of  1 1  Oct  is  not 
included  as  only  two  ems  are  used  for  profile  fitting. 


Correlation 

95% 

Confidence 

Band 

Number 

Student's  t 

coefficient 

of  ems 

distribution 

between 

Shear 

Apparent 

Bed  shear 

Day 

Run 

used  for 

observations 

stress 

roughness 

stress 

the 

?(K-2,l-o/2) 

and  log-profile 

velocity 

height 

Coefficient 

regression 

predictions 

V, 

za 

Cf 

n 

C 

(%) 

(X/-T) 

J 

(%) 

10 

1 

7 

2.015 

0.983 

17 

2.3 

33 

10 

2 

5 

2.353 

0.953 

43 

7.2 

99 

10 

3 

5 

2.353 

0.952 

44 

7.4 

102 

10 

4 

6 

2.132 

0.992 

14 

1.9 

27 

10 

5 

6 

2.132 

0.992 

13 

1.9 

26 

10 

6 

5 

2.353 

0.999 

5 

1.3 

6 

10 

7 

4 

2.920 

0.991 

28 

3.5 

56 

11 

1 

6 

2.132 

0.969 

27 

3.6 

58 

11 

2 

5 

2.353 

0.963 

38 

5.8 

87 

11 

3 

5 

2.353 

0.975 

31 

4.1 

68 

11 

4 

6 

2.132 

0.987 

17 

2.4 

35 

11 

5 

6 

2.132 

0.991 

15 

2 

28 

11 

6 

5 

2.353 

0.998 

9 

1.5 

14 

11 

7 

5 

2.353 

0.996 

12 

1.8 

21 

12 

1 

6 

2.132 

0.966 

29 

3.9 

59 

12 

2 

4 

2.920 

0.968 

53 

10 

128 

12 

3 

5 

2.353 

0.989 

20 

2.5 

41 

12 

4 

6 

2.132 

0.998 

6 

1.3 

10 

12 

5 

6 

2.132 

0.997 

8 

1.5 

14 

12 

6 

5 

2.353 

0.995 

14 

1.9 

25 

12 

7 

5 

2.353 

0.992 

17 

2.2 

30 

64 


Table  3.  Sensitivity  of  the  ems  to  elevation  from  the  bed  as  inferred  from  the  zero-shifting 
method. 
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III.  MEAN  CROSS-SHORE  CURRENTS  OVER  A  BARRED  BEACH 


(This  chapter  consists  of  work  to  be  submitted 
to  the  Journal  of  Geophysical  Research) 
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CROSS-SHORE  MEAN  FLOW  OVER  A  BARRED  BEACH 

A.F.  Garcez  Faria,  E.B.  Thornton,  T.C.  Lippmann1,  T.P.  Stanton 

Naval  Postgraduate  School,  Monterey,  California  93940 

ABSTRACT 

The  spatial  distribution  of  the  mean  cross-shore  flow  (undertow)  over  a  barred  beach 
is  examined  with  field  data  obtained  on  three  energetic  wave  days  during  the  DUCK94 
experiment.  The  vertical  structure  of  the  undertow  is  modeled  using  a  turbulent  eddy 
viscosity  closure  and  includes  the  important  effects  of  wave  breaking  (described  using  the 
roller  concept)  and  convective  acceleration  of  the  current.  Other  than  a  more  realistic 
description  of  observed  turbulence  variations,  a  depth  dependent  eddy  viscosity  does  not 
improve  the  agreement  between  predicted  and  observed  undertow  profiles.  The  effects  of 
using  different  boundary  conditions  is  investigated  by  extending  the  monochromatic 
formulations  of  Stive  and  Wind  (1986)  and  Svendsen  et  al.  (1987)  to  random  waves  by 
ensemble  averaging  over  the  wave  height  distribution.  The  contribution  of  breaking  wave 
rollers  to  the  surface  mass  flux  can  be  of  the  same  order  of  the  contribution  associated  with 
the  organized  wave  motion  and  should  not  be  neglected.  The  largest  discrepancies  between 
model  predictions  and  observations  occur  over  the  sand  bar,  where  the  observed  strong 
undertow  jet  is  not  predicted  by  the  model.  It  is  postulated  that  the  downward  momentum 
mixing  induced  by  wave  breaking  produces  a  constricted  return  flow  layer  that  enhances 
undertow  velocities,  a  process  not  accounted  for  in  the  model.  Better  agreement  is  obtained 
for  this  region  with  an  empirical  relation  that  correlates  the  surface  mass  flux  with  the  mean 
rate  of  wave-energy  dissipation  per  unit  area. 
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INTRODUCTION 

The  influence  of  the  mean  cross-shore  return  flow,  or  undertow,  on  sediment 
transport  and  beach  profile  evolution  has  long  been  believed  to  be  important  (e.g.,  Bagnold, 
1940).  A  local  vertical  imbalance  between  the  wave  setup  pressure  gradient,  that  is  uniform 
with  depth,  and  the  depth-varying  wave  radiation  stress  is  conceptually  responsible  for 
driving  the  undertow  (Dyhr-Nielsen  and  S0rensen,  1970).  In  the  last  two  decades  several 
theoretical  models  for  the  vertical  structure  of  the  undertow  for  two-dimensional  beaches 
have  been  developed  (e.g.,  Dally,  1980;  Svendsen,  1984;  Dally  and  Dean,  1984;  Stive  and 
Wind,  1986;  Svendsen  and  Buhr  Hansen,  1988;  Deigaard  et  ai,  1991;  Stive  and  DeVriend, 
1994;  Haines  and  Sallenger,  1994).  All  models  use  an  eddy  viscosity  closure  scheme,  and 
solve  for  the  depth  dependent  undertow  by  integrating  the  cross-shore  momentum  equation 
twice  over  depth,  which  requires  two  boundary  conditions  to  evaluate  the  integration 
constants. 

There  is  a  general  consensus  throughout  the  literature  of  using  local  conservation  of 
mass  over  the  vertical  as  one  boundary  condition.  Commonly,  the  second  boundary  condition 
is  either  the  stress  at  the  trough  level  (Stive  and  Wind,  1986)  or  the  no-slip  condition  at  the 
bottom  combined  with  the  steady  streaming  generated  by  the  bottom  boundary  layer 
(Svendsen,  1984).  Despite  significant  physical  differences,  both  approaches  reduce  to  the 
same  form  between  the  trough  level  and  the  top  of  the  bottom  boundary  layer  within  the  surf 
zone,  as  contributions  from  steady  streaming  and  bed  shear  stress  are  outweighed  by  mean 
water  slope  and  wave  forcing  gradients. 

The  wave-induced  onshore  mass  flux  in  the  region  between  the  wave  crest  and  trough 
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is  critical  to  predicting  the  magnitude  of  the  undertow,  which  is  predicted  heuristically  by 
adding  the  contribution  from  breaking  wave  rollers  to  the  mass  transport  given  by  an  inviscid 
wave  theory  plus  (Svendsen,  1984).  Dally  and  Brown  (1995)  found  better  agreement  for 
laboratory  generated  regular  waves  breaking  over  a  planar  beach  using  stream  function  theory 
(Dean,  1974)  compared  with  linear  wave  theory  (LWT,  hereafter),  which  tended  to  over- 
predict  the  depth-averaged  undertow.  On  the  other  hand,  Masselink  and  Black  (1995),  using 
field  measurements  from  experiments  at  two  planar  beaches,  found  good  agreement  using 
shallow  water  LWT. 

Observations  of  onshore  mass  transport  by  waves  in  the  crest-trough  region  are  scarce 
owing  to  difficulties  associated  with  making  measurements  in  this  region.  As  a  result,  the 
performance  of  surface  mass  flux  models  is  usually  evaluated  by  assuming  that  mass  is 
conserved  over  the  vertical  and  comparing  model  predictions  with  measured  return  mass 
flux.  This  procedure  dictates  a  priori  determination  of  the  height  of  the  cross-shore  return 
flow  layer.  Direct  determination  of  the  height  of  this  layer  requires  measurements  of  the 
vertical  distribution  of  mean  cross-shore  flow  with  a  vertical  resolution  that  so  far  has  only 
been  achieved  in  laboratory  experiments  with  monochromatic  waves  using  Laser  Doppler 
Velocimeters  {e.g.,  Nadaoka  et  al.,  1989). 

Most  existing  undertow  models  show  good  agreement  with  laboratory  data  for 
monochromatic  waves  breaking  over  planar  beaches  when  the  depth-averaged  mean  return 
flow  is  adjusted  to  fit  the  data  (instead  of  using  predicted  mass  flux)  and  the  magnitudes  of 
the  two  largest  dynamical  forcing  terms  (wave  setup  and  radiation  stress  gradients)  are 
determined  from  data. 
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Validation  of  these  models  with  field  data  has  been  limited  by  the  lack  of  data.  Smith 
et  al.  (1992),  compared  a  model  based  on  Putrevu  and  Svendsen  (1991)  solution  for  the 
vertical  structure  of  the  undertow  to  data  from  the  1990  DELILAH  experiment  on  a  barred 
beach  near  Duck,  North  Carolina,  and  found  large  discrepancies  over  the  bar,  where  the 
model  under-predicted  observed  velocities.  This  strong  "undertow  jet"  over  the  bar  was  also 
present  during  an  earlier  field  experiment  (Sallenger  et  al.,  1983)  at  the  same  site. 

Within  the  surf  zone,  wave-breaking  generated  turbulence  dominates  bottom 
boundary  layer  processes,  with  the  turbulent  shear  stress  maximum  at  the  surface  decreasing 
to  a  minimum  at  or  near  the  bottom.  King  and  Kirby  (1994),  using  LDV  data  from  a  wave 
flume  experiment,  showed  that  the  primary  turbulence  generating  mechanism  in  the  surf  zone 
is  due  to  wave-breaking  at  the  wave-roller  interface  and  that  turbulence  intensities  decrease 
with  distance  from  the  surface.  Cox  and  Kobayashi  (1997),  using  laboratory  LDV 
measurements  of  regular  waves  spilling  on  a  rough,  plane  slope,  observed  that  the  shear 
stress  distribution  within  the  surf  zone  varies  linearly  with  depth  until  the  top  of  the  bottom 
boundary  layer,  and  that  the  eddy  viscosity  is  small  near  the  trough  level,  increases  to  a 
maximum  about  one-third  of  the  depth  below  the  trough  level,  and  then  decreases  toward  the 
bottom.  Svendsen  (1984),  investigated  the  effect  of  introducing  an  exponentially  varying 
eddy  viscosity  (u  )  to  the  predicted  vertical  profile  of  the  undertow.  In  his  formulation,  two 
free  parameters  are  necessary  to  define  the  magnitude  and  decay  rate  with  depth  of  u,,  and 
despite  the  more  realistic  depth-variation  of  uT,  only  marginal  improvements  are  obtained 
for  the  vertical  profile  of  the  undertow  (Svendsen  and  Buhr  Hansen,  1988). 

In  this  paper,  field  observations  of  vertical  profiles  of  the  mean  cross-shore  current 
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obtained  over  a  naturally  barred  beach  are  used  to  test  various  models.  Wave  height 
transformation  is  based  on  an  energy  flux  balance  including  the  effects  of  rollers  to  describe 
wave  breaking  and  a  probabilistic  description  of  wave  heights.  The  mass  flux  is  investigated 
using  both  linear  and  higher-order  wave  theories,  and  including  the  contribution  from  surface 
rollers,  to  estimate  the  cross-shore  distribution  of  the  undertow.  The  influence  of  depth- 
dependent  formulations  for  the  eddy  viscosity  and  different  choices  of  boundary  conditions 
on  the  vertical  structure  of  the  undertow  is  examined. 

THEORY 

In  the  following,  a  right-handed  coordinate  system  with  origin  at  the  shoreline  (x 
positive  offshore)  and  z  positive  upward  from  the  sea  surface  is  used.  Solutions  for  the 
surface  mass  flux,  setup  and  the  vertical  profile  of  the  undertow  are  described  assuming 
stationary  wave  conditions,  straight  and  parallel  depth  contours,  and  random  waves  that  are 
narrow-banded  in  both  frequency  and  direction. 
Wave  Transformation 

Wave  transformation  is  determined  from  a  probabilistic  breaking  wave  model  that 
includes  roller  energy  gradients  in  the  energy  flux  balance  (Lippmann  et  al.,  1996,  LBT 
hereafter).  This  model  assumes  that  wave  heights  both  inside  and  outside  the  surf  zone  can 
be  reasonably  described  by  the  Rayleigh  distribution  (Thornton  and  Guza,  1983)  and  has 
been  shown  to  give  accurate  results  for  random  waves  breaking  over  both  planar  and  barred 

beaches.  The  model  has  two  free  parameters:  o,  the  mean  angle  of  the  wave-roller  interface, 

H 

and  y .  a  measure  of  breaking  wave  intensity  y  =  — — ,  where//      is  the  rms  wave  height 

h  "" 

and  h  is  the  local  depth.  Although  the  model  is  insensitive  to  the  interfacial  angle  (kept 
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constant  at  a  =  10° for  all  runs)  based  on  the  results  obtained  by  LBT  at  the  same  location 
of  the  DUCK94  experiment,  it  is  necessary  to  adjust  y  to  fit  the  data. 

Wave  properties  are  ensemble  averaged  by  integrating  wave  heights  through  the 
assumed  Rayleigh  distribution,  p(H).  Following  Thornton  and  Guza  (1983),  the  fraction  of 
waves  that  are  breaking  is  found  by  integrating  through  the  breaking  wave  distribution 
pb(H)  =  W(H)  p(H),  where  the  weighting  function,  W(H),  is  given  by  Whitford  (1988). 
Surface  Mass  Flux 

The  conservation  of  mass  for  straight  and  parallel  contours  with  the  boundary 
condition  of  no  flow  through  the  beach  is  given  by 


|   p[U(z)  +  u(z)+u(z)]dz  =  0  (1) 

-h 

where  overbar  indicates  time  averaging,  x\  is  surface  elevation,  p  is  water  density,  and  the 
onshore  horizontal  velocity  has  been  partitioned  into  mean,  wave  and  turbulence 
contributions.  In  an  Eulerian  reference  frame  and  assuming  irrotational  flow  below  the 
trough  level,  there  is  a  net  shoreward  mean  mass  transport  by  waves  limited  to  an  upper 
region  between  the  crest  and  trough  that  is  given  by 


4w  =    /    pu(z)dz  (2) 

n, 

where  subscripts  (c,  t)  refer  to  (crest,  trough).  There  is  an  additional  contribution  to  the  mean 
mass  transport  (per  surface  area)  above  the  wave  trough  (surface  layer)  that  arises  from  the 
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presence  of  turbulent  wave  rollers  (Svendsen,  1984),  that  is  defined  by 

4r    =    PrA  f-  (3) 


where  pr  is  the  roller  density,  Lris  the  length  along  the  wave  covered  by  the  roller  cross- 
sectional  area  A,  and  c  is  the  velocity  that  the  roller  mass  is  advected  landward,  assumed  to 
be  given  by  the  wave  phase  speed.  The  roller  contribution  (3)  to  the  surface  mass  flux  is 
based  on  the  calculated  cross-sectional  area  of  the  roller.  Earlier  models  (Svendsen,  1984; 
Deigaard  et  al.,  1991)  assume  that  this  area  is  proportional  to  the  rms  wave  height.  As  a 
consequence,  the  largest  mean  return  flow  would  be  predicted  to  occur  at  the  breaker  point, 
which  is  contrary  to  laboratory  observations  (Nadaoka  and  Kondoh,  1982).  Dally  and  Brown 
(1995)  solve  for  the  area  of  the  roller  for  the  case  of  monochromatic  waves  breaking  over  a 
planar  beach  by  simultaneously  solving  the  depth-integrated  and  time-averaged  energy, 
continuity,  and  cross-shore  momentum  equations  and  found  good  agreement  with  existing 
laboratory  data.  However,  their  model  requires  observations  of  the  cross-shore  distribution 
of  either  the  mean  return  flow  or  setup  to  constrain  the  model,  and  therefore  cannot  be 
applied  to  the  present  data. 

Another  recent  roller  model  (Lippmann  and  Thornton,  1997;  LT97  hereafter)  presents 
an  independent  method  for  calculating  qr,  and  is  calibrated  with  video  observations  of  the 
cross-shore  variation  of  the  fraction  of  waves  that  are  breaking.  The  LT97  model  is  based  on 
the  energy  flux  balance  and  describes  energy  dissipation  following  Deigaard  (1993).  Two 
free  parameters  in  the  model,  B,  the  vertical  fraction  of  wave  height  covered  by  the  roller, 
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and  ijj ,  a  measure  of  the  average  wave  face  angle,  are  adjusted  to  give  a  rms  best  fit  to 
breaking  observations. 

The  onshore  mass  transport  in  the  upper  region  is  balanced  locally  by  a  mean  return 
flow  below  the  trough  (undertow) 


n, 


qw  +  qr  =  ~  (  pU(z)dz  =  -puRht  (4) 

-h 

where  UR  is  the  depth-averaged  return  flow  and  ht  is  the  depth  below  the  trough. 
Undertow 

The  time-averaged  cross-shore  momentum  equation,  neglecting  molecular  viscous 
stresses,  and  for  the  assumed  straight  and  parallel  contours  and  stationary  wave  conditions, 
is  given  by 


dpu~        dpuw       _  dp 
dx  dz  dx 


The  horizontal  and  vertical  velocities  (u,w)  are  expanded  into  mean,  turbulent,  and  wave- 
induced  components,  u  =  U  +  u  +  u  and  w  =  w  +  w ,  where  the  mean  vertical  velocity  is 
assumed  equal  to  zero.  The  time-averaged  pressure,  p  =  p  [g(r\  -  z)  -  vv  2  -  w  2] ,  is  obtained 
from  the  depth  integrated  vertical  momentum  equation,  after  neglecting  contributions  from 
the  cross-shore  gradient  of  the  vertically  integrated  wave  and  turbulent  shear  stress  (Stive  and 
Wind,  1982).  After  substituting  for  the  time-averaged  pressure,  (5)  can  be  further  simplified 
for  the  region  between  the  top  of  the  bottom  boundary  layer  and  trough  level  (middle  layer), 
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with  the  aid  of  the  following  assumptions:  1)  wave  and  turbulent  velocity  components  are 
statistically  independent;  2)  turbulence  is  near  isotropic  (Stive  and  Wind,  1982);  3)  the  wave 

shear  stress  is  given  by  — —  = -[u2  -w2]  (Ri  vero  and  Arcilla,  1 995);  and  4)  a  first 

dz  2  dx 

dU 

dz 
assumptions,  (5)  reduces  to 


order  eddy  viscosity  closure  for  the  turbulent  shear  stress  -  p  uw  =  p  u,-^— .  Applying  these 


d  ,       dU ..        1  dp  r—    — ,  dr\        dpU2       „,  .  ... 

—  [pu— -]  =  --^[m-w2]   +  pg— !■  +  -£—  =  F(x)  (6) 

dz        "  dz  2  ox  ox  ox 


where  u  is  the  time-invariant  turbulent  eddy  viscosity.  The  forcing,  F(x),  in  (6)  is  due  to  the 
cross-shore  gradients  of  radiation  stress,  setup  (setdown),  r\ ,  and  convective  acceleration  of 
the  depth-averaged  undertow  respectively.  Here,  F(x)  is  assumed  to  be  independent  of  depth 
based  on  empirical  evidence  from  laboratory  studies  (Nadaoka  and  Kondoh,  1982;  Stive  and 
Wind,  1982  and  1986). 

A  solution  for  the  vertical  distribution  of  the  mean  undertow  can  be  determined  by 
integrating  (6)  twice  over  depth,  to  give 


U(z)  =  F{x)  f-Z-dz  +  Cx(x)  [-4*-  *  C2(x)  (7) 

Jon  Jou  v  ' 


where  C,(x)  and  Cn{x)  are  spatially  varying  integration  constants.  Both,  the  Stive  and  Wind 
(1986)  and  Svendsen  et  al.  (1987)  solutions  are  used,  with  modifications  to  their  original 
developments  to  improve  the  derivation.  In  both  these  original  works,  a  different  solution  for 
the  undertow  forcing  (F)  is  used  as  they  assume  that  the  wave  contribution  to  stress  is  small 
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compared  with  the  Reynolds  stress  (puw  «  puw)  resulting  in  an  over  estimation  of  the 

forcing  by °  [w2  -vv2]  (Rivero  and  Arcilla,  1995).  Here,  (6)  will  be  used  instead. 

2  dx 

Although  different  formulations  for  the  eddy  viscosity  variation  with  depth  are 
investigated  (Appendix),  the  simplest  solution  for  a  depth-independent  eddy  viscosity  (u ) 
is  outlined  here.  For  this  case,  (7)  simplifies  to: 

U(z)  =  -1—F(x)z2  +  —Cx{x)z  +  C2(x)  (8) 

2pu  pu 

First,  following  Stive  and  Wind  (1986),  C^x)  is  determined  by  integrating  (6)  once 
over  depth  and  solving  for  the  shear  stress  at  the  trough  level,  and  C2(x)  by  applying 
conservation  of  mass  over  the  vertical  (see  Appendix  for  details)  to  give 


U(z)  =  UR  +  —F0c)gr+hz+  —  ]  +  —  [z  +  h  ■    -1  ]       (9) 


,  _ —  +  h  z  +  J  +  

pu  2  6  p(i 


where  xbx  is  the  mean  cross-shore  bed  shear  stress  that  is  calculated  using  a  quadratic 
formulation  (zbx  =  p  C,  ub  Ub,  where  Cf  =  0.01  is  a  constant  friction  factor  and  ub  is  the 
amplitude  of  the  near-bottom  wave-induced  velocity).  This  monochromatic  solution  is 

quadratic  in  z  with  coefficients  independent  of  wave  height  A2  =  — ,     A,  =  h,  and 

3/2 2  -  hf 

A0=  ,  which  can  be  extended  to  random  waves  by  ensemble  averaging  over  the 

6 

wave  distribution 


eo  7 

<U(z)>  =  UR  *  (A,z2  +A,z  +A0)—    f  F{x)p(H)dH  +  -Z*  (z  +  *  --  -i )  (10) 

pu    {  PM  2 
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In  the  cross-shore  direction,  xbx  is  small  compared  with  wave-breaking  dissipation  (Thornton 
and  Guza,  1983;  Svendsen,  1984;  Dally  and  Brown,  1995),  and  hence  the  vertical  structure 
of  the  undertow  is  mostly  determined  by  the  coefficients  A2,  Ajand  AQ.  The  amount  of 
curvature  in  the  profile  is  a  function  of  both  F  (calculated  from  wave  quantities)  and  the  eddy 
viscosity,  u .  Large  values  of  F  produce  more  vertical  shear  resulting  in  a  parabolic  profile, 
whereas  large  values  of  u  reduce  the  vertical  shear  producing  a  more  uniform  velocity 
profile  with  depth. 

The  influence  of  the  boundary  condition  choice  to  the  vertical  structure  of  the  mean 
undertow  is  investigated  by  comparing  the  results  obtained  by  (10)  with  the  Svendsen  et  al. 
(1987)  model  that  uses  a  no-slip  condition  at  the  bottom  to  replace  the  stress  at  the  trough 
level  as  the  second  boundary  condition.  This  no-slip  condition  is  obtained  by  coupling  the 
undertow  model  with  a  bottom  boundary  layer  model  (BBL,  hereafter).  Within  the  BBL,  the 
flow  is  a  combination  of  the  steady  streaming  induced  by  the  oscillatory  motion  and  the 
undertow  above  the  BBL,  which  results  in  a  mean  velocity  at  the  top  of  this  layer  ( Ub)  that 
is  obtained  by  requiring  continuity  in  velocity  and  shear  stress  between  these  two  regions 
to  give  (see  Svendsen  et  al.,  1987  for  details) 


<U(z)>  =  Ub  +  (A2z2  +A]Z  +A3)  J-    f  F(x)p{H)dH  +  ^^  (z  +  h)      (11) 

PM    {  PM 


h2 
where  A3  =  —  is  another  coefficient  independent  of  wave  height  and  xfa  is  related  to  the 


steady  streaming  in  the  BBL  ( xbs  =  --ppw 


\ 
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,  where  /  is  wave  frequency  and  uw 


is  the  eddy  viscosity  inside  the  BBL).  Their  model  further  assumes  that  uw  is  much  smaller 
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than  the  eddy  viscosity  in  the  middle  layer,  which  was  recently  verified  in  a  laboratory 
experiment  using  LDV  (Cox  and  Kobayashi,  1997).  Following  Putrevu  and  Svendsen  (1993), 

(J  ii  2 

the  eddy  viscosity  inside  the  BBL  is  estimated  by =  0.32  C,2  ( — )  . 

hyfgh  h 

Setup 

The  setup  gradient  is  a  dominant  driving  force  for  the  undertow.  The  setup  is 
calculated  by  depth-integrating  the  time-averaged  cross-shore  momentum  equation  (5)  from 
the  bottom  to  the  mean  water  level,  to  give 

dS„      dMr  _         dUR2  -         9= 

-^  +-T1  +  P(T|  +  *)  — -  =   -Pgm+h)?p-  (12) 

ox         ox  ox  ox 

where  S^  is  the  wave  radiation  stress  (Longuet-Higgins  and  Stewart,  1964),  Mr  =  cqr  is 
the  momentum  flux  associated  with  wave  rollers,  and  the  last  term  on  the  lhs  is  the 
convective  acceleration  of  the  depth-averaged  undertow.  In  the  derivation  of  (12),  it  is 
assumed  that  the  mean  bed  shear  stress  (Longuet-Higgins  and  Stewart,  1964)  and  the  cross- 
shore  wind  stress  (owing  to  small  observed  onshore  winds  in  the  data  described  later)  are 
negligible. 

DATA 
Field  measurements  were  acquired  as  part  of  the  DUCK94  experiment  {e.g.  Garcez 
Faria  et  al,  1997  [Chapter  1])  conducted  at  the  U.S.  Army  Corps  of  Engineers  Field 
Research  Facility  (FRF)  in  Duck,  North  Carolina.  The  data  selected  for  analysis  are  from  10- 
12  Oct.  when  strong  cross-shore  currents  (0.05-0.4  m/s)  caused  by  a  storm  with  predominant 
winds  and  waves  from  the  north  were  present. 
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The  vertical  structure  of  the  current  was  measured  with  a  vertical  stack  of  seven  two- 
component  Marsh-McBirney  electromagnetic  current  meters  (ems,  hereafter)  mounted  on  a 
mobile  sled  at  elevations  0.4,  0.7,  1.0,  1.5,  1.8,  2.2,  and  2.6  m  above  the  bed.  The  ems  were 
horizontally  displaced  at  least  1  m  from  the  sled,  and  oriented  such  that  the  vertical  stack  was 
in  the  up-drift  direction  of  the  longshore  current  to  minimize  interference  by  the  sled 
structure.  The  em  offsets  were  determined  in  situ  and  found  to  be  accurate  within  1  cm/s. 

The  sled  orientation  was  determined  using  a  digital  compass  mounted  on  the  sled 
with  accuracy  0(1  degree).  Measured  two-component  velocities  were  rotated  to  a  shore 
normal  right-handed  coordinate  system  by  using  compass  data  and  adding  at  each  sled 
position  any  deviation  of  the  bottom  contour  line  from  a  shore  parallel  direction.  A  precise 
determination  of  the  rotation  angle  is  important  to  avoid  contamination  of  cross-shore 
velocities  by  the  observed  strong  longshore  currents  0(1  m/s).  Velocity  errors  associated 
with  the  rotation  of  the  coordinate  system  were  comparable  to  the  determined  offsets,  and 
therefore  are  neglected. 

For  the  first  run  on  each  day,  the  sled  was  towed  by  the  Coastal  Research  Amphibious 
Buggy  (CRAB)  to  its  farthest  offshore  location  seaward  of  the  bar  (approximately  160  m 
from  the  shoreline).  A  forklift  on  the  beach  pulled  the  sled  shoreward  10  to  30  meters  for 
subsequent  measurement  runs  that  are  referred  to  in  the  text  by  sequential  numbers  within 
each  day.  Each  data  run  was  nominally  one  hour,  and  seven  to  eight  runs  were  made  across 
a  transect  during  each  day  spanning  the  high  tide  during  this  period. 

Waves  and  mean  water  level  were  measured  using  an  array  of  five  pressure  sensors 
mounted  on  the  sled.  Directional  wave  spectra  were  also  acquired  using  a  linear  array  of  10 
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pressure  sensors  in  8  m  depth.  Additionally,  a  13  element  cross-shore  array  of  pressure 
sensors  was  used  to  measure  wave  heights  spanning  the  width  of  the  surf  zone  (Elgar  et  al., 
1997).  The  fixed  array  was  located  approximately  25  m  to  the  north  of  the  sled  transect. 
These  data  were  continuously  sampled  at  2  Hz.  Video  observations  were  used  to  measure  the 
fraction  of  wave  breaking  along  the  same  transect  using  the  method  of  Lippmann  and 
Holman  (1991). 

Meteorological  information  of  wind  and  atmospheric  pressure  were  recorded 
simultaneously  at  the  seaward  end  of  the  600  m  long  FRF  pier  and  atop  the  FRF  building  in 
front  of  the  pier.  The  bathymetry  was  measured  daily  using  the  CRAB  and  it  was  found  that 
the  depth  contours  were  nearly  straight  and  parallel  for  the  three  days  under  consideration. 

RESULTS  AND  DISCUSSION 

Measured  and  modeled  undertow  flows  are  maximum  on  top  and  on  the  shoreward 
slope  of  the  bar,  with  increasing  magnitude  from  10  to  12  Oct.  as  wave  forcing  increased 
with  the  approach  of  a  storm  (Fig.  1).  The  vertical  structure  over  this  region  is  the  classic 
parabolic  shape  associated  with  strong  wave-breaking  turbulence  (Svendsen,  1984;  Stive  and 
Wind,  1986).  In  the  inner  trough,  the  return  flow  is  weak,  and  almost  no  vertical  structure 
is  seen.  At  the  seaward  slope  of  the  bar,  observed  profiles  on  10  Oct.  (Fig.  la)  are  nearly 
uniform  with  depth.  It  is  noted  there  is  significant  bottom  roughness  in  the  trough  of  the 
barred  profiles  related  to  the  development  of  lunate  and  long-crested  mega-ripples  as  a  result 
of  the  strong  longshore  currents  (Thornton  et  al.,  1997).  The  bar  migrated  off-shore  during 
this  period  O(20  m),  which  was  found  by  Gallagher  et  al.  (1997)  to  be  associated  with  the 
strong  undertow. 
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Model/data  comparisons  are  evaluated  by  calculating  absolute  and  relative  root- 
mean-square  (rms)  errors.  Absolute  error  has  dimensional  units  and  is  defined  by 
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where  P.  is  model  prediction,  Oi  is  the  observed  quantity,  and  N  is  the  number  of 
observations.  The  relative  or  percent  error  is  calculated  by 


e-i  =  100  [ 
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-£(- '-)   ]  (14) 


and  weights  the  difference  between  model  predictions  and  observations  against  a  measure 
of  expectation  that  is  represented  by  the  observation.  This  statistic  is  not  well-behaved  for 
small  Or  which  sometimes  is  the  case  for  observed  mean  cross-shore  velocities.  Thus 
measurements  smaller  than  the  in  situ  determined  offsets  (+/-  1  cm  s"1)  are  excluded  from 
(14). 

In  the  following,  the  accuracy  of  the  wave  transformation  model  is  described  first, 
which  is  then  used  to  predict  the  surface  mass  flux,  setup  and  vertical  profiles  of  the 
undertow  including  contributions  from  surface  rollers. 
Wave  Transformation 

The  rms  wave  height  is  approximated  by  H  =  y8o2  where  o2  is  the  variance 
calculated  from  the  surface  elevation  time  series.  Surface  elevation  was  calculated  by  Fourier 
transforming  a  one-hour  pressure  record,  applying  a  LWT  transfer  function  to  the  complex 
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Fourier  amplitudes  in  the  frequency  domain,  band-pass  filtering  by  zeroing  coefficients 
outside  the  range  of  interest  (0.05Hz  <  f  <  0.5  Hz),  and  inverse  transforming  to  obtain  the 
surface  elevation  time  series  (Thornton  and  Guza,  1982). 

The  sensitivity  of  the  models  for  surface  mass  flux,  setup,  and  undertow  to  errors 
associated  with  the  use  of  a  random  wave  transformation  model  (LBT)  is  investigated  by 
comparing  with  results  using  measured  wave  heights  interpolated  with  cubic  splines.  The 
difference  between  model  output  and  data  interpolation  methods  is  shown  in  Fig.  2  for  the 
first  run  of  1 1  Oct.,  which  has  the  largest  rms  relative  error  (7  percent). 

A  plot  of  measured  versus  modeled  Hrms  for  the  three  days  analyzed  (10-12  Oct.) 
shows  good  agreement  (Fig.  3).  The  mean  zrel  is  5  percent  for  all  runs,  with  largest  errors  at 
any  cross-shore  position  within  1 2  percent.  Best  fit  values  for  y  and  z  {  are  summarized  in 
Table  1 .  In  general  the  transformation  model  well  represents  measured  wave  heights. 
Surface  Mass  Flux 
Estimating  Surface  Mass  Flux  from  velocity  measurements 

The  surface  mass  flux  is  inferred  from  the  measured  return  flow,  which  is  sensitive 

to  the  depth  of  the  return  flow  layer.  For  field  experiments  as  well  as  for  modeling  efforts, 

the  height  of  the  return  flow  layer  is  assumed  to  be  coincident  with  the  trough  level,  ht .  Here, 

ht  is  calculated  using  both  linear  and  non-linear  wave  theories  and  compared  to  surface 

elevation  estimates.  Laitone's  second  order  approximation  (Laitone,  1961)  to  cnoidal  wave 

theory  (Korteweg  and  de  Vries,  1895)  is  used  as  this  theory  provides  an  analytical  solution 

to  h  ,  and  the  conditions  for  all  runs  being  examined  are  within  the  range  of  validity  of  this 

theory  determined  by  the  Ursell  number,  — >  0.6,  where  k  is  wavenumber  (Laitone, 

(kh)2 
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1962).  The  use  of  a  non  linear  wave  theory  generally  gives  higher  trough  levels  than  LWT, 
and  consequently  smaller  mean  return  flows  ( UR )  for  the  same  onshore  mass  transport  (Fig. 
4).  However,  the  effect  on  the  calculated  total  return  mass  flux  is  minimal  as  an  increase  in 
ht  is  counterbalanced  by  a  decrease  in  UR .  No  statistically  significant  difference  at  the  97 
percent  confidence  level  is  found  between  the  results  obtained  using  cnoidal  and  linear  wave 
theory  for  all  22  runs  (Fig.  5). 

A  problem  associated  with  practical  applications  of  cnoidal  wave  theory  is  related 
to  its  mathematical  complexity,  as  it  contains  both  the  Jacobian  elliptic  functions  and  the 
complete  elliptic  integrals  of  the  first  and  second  kinds.  An  attractive  alternative  is  the  use 
of  its  asymptotic  approximation,  or  hyperbolic  wave  theory  (Iwagaki,  1968),  which  is  valid 
when  the  complete  integral  of  the  second  kind  (K)  is  greater  than  3,  and  the  Ursell  number 
is  greater  than  0.6  (Svendsen,  1974),  which  is  the  case  for  the  data  examined  here.  The 
calculated  difference  between  observed  surface  mass  flux  given  by  LWT  and  Iwagaki' s 
approximation  is  zrd  =  7.9  percent,  indicating  that  the  hyperbolic  wave  theory  is  accurate 
enough  to  replace  the  more  complex  cnoidal  wave  theory.  Applying  Iwagaki 's  solution,  an 
expression  for  ht  similar  to  LWT  is  derived  that  is  a  function  of  only  local  depth  and  rms 
wave  height. 


1  H 

'       2  3  \ 


H      2       H* 
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The  observed  mean  trough  level  is  obtained  by  applying  a  zero-up-cross  method  to 
the  measured  surface  elevation  time  series  to  define  individual  waves  and  then  averaging  the 
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minima  between  each  up-cross  interval.  A  comparison  with  measured  mean  trough  level 
(Fig.  6),  shows  that  this  simple  expression  (15)  provides  a  better  estimate  for  the  mean 
trough  level  within  the  surf  zone  (rms  over-prediction  of  4  cm)  than  LWT  (rms  under- 
prediction  of  17  cm).  It  can  also  be  noticed  in  Fig.  6  that  discrepancies  between  predicted  and 
observed  mean  trough  level  increase  with  rms  wave  height  for  LWT,  while  (14)  provides  an 
estimate  of  ht  that  does  not  appear  sensitive  to  H . 
Surface  Mass  Flux  model 

The  inferred  surface  mass  flux  obtained  by  numerically  integrating  the  observed 
cross-shore  mean  currents  from  the  bottom  to  the  trough  level  given  by  (14)  are  compared 
to  model  estimates  (4).  qr  is  determined  from  the  calibrated  LT97  model  and  qw  is 
determined  first  using  linear  wave  theory  and  then  with  non  linear  theory  for  comparison. 
Predicted  qr  are  dependent  on  the  percentage  of  breaking  waves  in  the  LT97  model.  For  the 
22  profiles  examined,  best  fit  values  of  B  range  0.60-0.95  and  i|r  range  1.1-1.9  (Table  1),  and 
the  error  between  observed  and  modeled  percentage  of  breakers  is  trel=  5.6  percent.  An 
example  of  observed  and  predicted  percentage  of  breakers  is  shown  in  Fig.  7  for  the  third  run 
of  12  Oct.,  which  has  the  largest  rms  error  (9  percent). 

Inferred  and  predicted  surface  mass  flux  using  linear  wave  theory  are  compared  for 
the  entire  ensemble  in  Fig.  8.  Due  to  significant  changes  associated  with  wave-breaking 
characteristics,  the  data  have  been  divided  into  four  regions:  shoreward  slope  of  the  bar, 
seaward  slope  of  the  bar,  trough,  and  foreshore.  The  rms  relative  error  between  observations 
and  predictions  using  only  qw,  given  by  linear  theory  is  40  percent.  Including  the  mass  flux 
contribution  from  wave  rollers,  qr,  improves  the  overall  agreement  (rms  error  =  28  percent). 
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Largest  values  of  qr  occur  at  the  shoreward  slope  of  the  bar  (Fig.  8,  lower  panel),  where 
wave-breaking  is  most  intense.  In  this  region,  qr  is  on  average  72  percent  of  qw,  but  can  be 
as  large  as  144  percent  (fourth  run  of  12  Oct.).  This  is  in  accordance  with  earlier  model 
studies  (Svendsen,  1984;  Dally  and  Brown,  1995),  but  contrary  to  Masselink  and  Black 
(1995)  who  contend  that  the  roller  contribution  to  the  mass  transport  is  of  secondary 
importance  based  on  the  results  from  field  experiments  at  two  near  planar  beaches. 

The  impact  of  using  a  non  linear  wave  theory  to  calculate  qw  is  examined  next.  For 
this  investigation,  stream  function  theory  (SFT,  hereafter)  is  used  as  it  provides  a  higher 
order  solution  for  qw  than  cnoidal  wave  theory,  and  it  has  already  been  shown  to  give 
accurate  results  with  laboratory  data  (Dally  and  Brown,  1995).  Surface  mass  fluxes  were 
calculated  using  SFT  both  including  and  excluding  qr .  For  the  field  data  examined,  there  is 
little  difference  between  the  use  of  LWT  or  SFT,  with  LWT  giving  values  on  average  8 
percent  larger  than  SFT  (Fig.  9). 

Errors  in  surface  mass  flux  predictions  could  also  arise  from  the  misfit  in  estimating 
wave  heights  with  a  transformation  model.  Calculated  rms  relative  error  in  mass  flux 
between  using  LBT  model  predicted  H  and  applying  a  cubic  spline  to  the  measurements 
for  all  runs  is  9  percent,  with  a  maximum  relative  error  of  19.5  percent  (third  run  of  12  Oct.). 
These  results  suggest  that  the  surface  mass  flux  is  not  overly  sensitive  to  the  choice  of  the 
random  wave  transformation  model. 
Undertow  Jets 

The  onshore  transport  is  observed  to  extend  below  the  mean  through  level  during 
high  wave  conditions  and  where  a  high  percentage  of  wave-breaking  occurs,  such  as  over 
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of  the  bar  (Fig.  4b),  which  violates  a  basic  assumption  of  existing  surface  mass  flux  models. 
Inferred  downward  momentum  mixing  produced  from  breaking-injected  turbulence  is 
probably  responsible  for  this  anomalous  onshore  flow  in  the  upper  layer  penetrating  beneath 
the  trough  level  (Masselink  and  Black,  1995;  Garcez  Faria  et  al.  1997  [Chaper  1]).  To 
conserve  mass,  undertow  velocities  in  this  constricted  return  flow  layer  are  enhanced 
resulting  in  an  undertow  jet.  This  phenomenon  was  qualitatively  observed  by  Nadaoka  et  al. 
(1989)  in  a  laboratory  experiment  using  flow  visualization  techniques  and  a  fibre-optic  LDV 
system.  In  their  experiment,  obliquely  descending  eddies  behind  the  breaking  wave  crest 
were  observed  to  substantially  transport  vorticity  into  an  otherwise  nearly  irrotational 
velocity  field,  resulting  in  an  increased  mass  transport. 

On  the  shoreward  slope  of  the  bar  where  strong  wave-breaking  was  observed  (Fig. 

7),  the  model  considerably  under-predicts  mass  transport  (Fig.  8).  The  effect  of  the 

q  ,     -  [q     +  q  ] 
constricted  return  flow  layer  is  estimated  for  this  region  by  — — .  Enhanced 

"obs 

undertow  velocities  account  for  an  average  of  39  percent  of  the  calculated  surface  mass  flux, 
and  are  in  qualitative  agreement  with  Nadaoka  et  al.  (1989). 

These  observed  strong  jets  are  not  believed  to  be  associated  with  three-dimensional 
circulation  cells,  as  measured  bathymetry  was  essentially  uniform  alongshore  during  the 
period  being  analyzed  and  no  qualitative  evidence  of  stationary  rip  currents  was  observed  in 
the  video. 
Empirical  estimate  for  Surface  Mass  Flux 

It  was  found  empirically  that  the  inferred  surface  mass  flux  (Q  =  p  Uobsht)  is 
significantly  correlated  at  the  99.5  percent  confidence  level  (linear  correlation  of  0.89  and 
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e  ,=  27  percent)  with  the  mean  rate  of  wave-energy  dissipation  per  unit  area,  D  (Fig.  10) 


.1/3  r»2/3 


G  =  5.27 — -I  (16) 
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The  spatially  varying  D  is  estimated  from  the  energy  balance  equation  using  the  LT97  model. 

This  empirical  relation  provides  a  simple  alternative  to  estimate  the  surface  mass  flux  within 

the  surf  zone  that  easily  can  be  incorporated  into  nearshore  current  models  for  field 

applications  with  the  same  accuracy  of  existing  surface  mass  flux  models.  A  better  agreement 

is  obtained  for  the  region  of  the  shoreward  slope  of  the  bar  with  this  relation  (Fig.  8-10). 

Setup 

Setup  is  calculated  using  a  finite  centered  difference  method  to  numerically  solve 

(12),  with  the  condition  that  the  setup  is  assumed  zero  at  the  most  offshore  grid  point. 

Contributions  from  each  term  of  (12)  are  examined  for  the  fourth  run  of  12  Oct.  (Fig.  1 1), 

which  corresponds  to  the  most  energetic  period.  The  cross-shore  gradient  of  the  momentum 

flux  associated  with  wave  rollers  is  calculated  using  the  LT97  model,  the  radiation  stress 

term  is  calculated  using  linear  wave  theory  (Longuet-Higgins  and  Stewart,  1964),  and  the 

convective  acceleration  of  the  current  (  — — )  is  estimated  by  applying  cubic  splines  to 

dx 

measured  depth-integrated  cross-shore  mean  return  flow.  The  largest  contributions  are  due 
to  the  radiation  stress  and  roller  terms.  The  convective  acceleration  term,  although  generally 
an  order  of  magnitude  smaller  than  the  other  terms,  cannot  be  neglected,  as  it  can  be  of  the 
same  order  as  the  sum  of  the  larger  terms. 

The  effects  on  setup  by  contributions  from  the  roller  momentum  flux  and  using  non 
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linear  wave  theory  to  calculate  the  radiation  stress  are  examined  with  one  example  from 
each  day  (Fig.  12).  The  effect  of  the  roller  is  to  redistribute  momentum  laterally  into  the 
trough  of  the  bar,  shifting  the  point  where  the  setup  begins  onshore.  This  shift  has  a 
significant  impact  on  the  setup/setdown  profile  within  the  surf  zone.  The  use  of  a  non  linear 
wave  theory  (SFT  is  used  here)  to  estimate  the  wave  radiation  stress  does  not  appear  to  alter 
the  setup  profile  dramatically  (Fig.  12). 
Undertow 
Cross-shore  variation  of  the  eddy  viscosity 

An  implicit  assumption  of  the  closure  model  used  in  the  undertow  solution  is  that  the 
eddy  viscosity  coefficient  ( u )  is  proportional  to  turbulence  intensity.  Conceptually,  large 
cross-shore  variation  of  this  parameter  is  expected  throughout  the  surf  zone  over  a  natural 
barred  beach  associated  with  significant  changes  in  wave-breaking  generated  turbulence. 
Haines  and  Sallenger  (1994),  in  an  earlier  field  experiment  limited  to  3  ems  over  the  vertical 
at  this  same  beach,  used  a  different  model  for  the  vertical  variation  of  the  undertow  and 
found  that  best  fit  u  for  each  location  varied  by  more  than  an  order  of  magnitude  across  the 
surf  zone  (u  =  0.0055  -  0.075  m2  s"1 )  in  qualitative  agreement  with  theoretical  expectations. 
Smith  et  al.  (1992)  on  the  other  hand  applied  a  constant  u  across  the  surf  zone  (u  =  0.05  m2 
s"1 )  and  found  good  agreement  with  the  DELILAH  data. 

In  order  to  investigate  the  cross-shore  variation  of  the  eddy  viscosity,  best  fit  u  (Table 
2)  were  calculated  using  (10)  by  minimizing  both  the  zabs  and  zrd  at  each  cross-shore 
position  occupied  by  the  sled  (Fig.  1).  For  all  22  runs,  zabs  =  2.2  cm  s"1  (erel  =  19  percent) 
with  maximum  8^=  6.1  cm  s  '  (zrel  =  44  percent)  for  the  fourth  run  of  12  Oct.  Observed 
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mean  undertow  velocities  ( Uobs)  and  splined  rms  wave  heights  were  used  in  the  calculations 
to  isolate  effects  due  to  the  choice  of  u .  Using  modeled  rms  wave  heights  (LBT96),  zabs  = 
2.5  cms"1  {zrel  =  22  percent),  which  indicates  that  the  vertical  profile  of  the  undertow  is  not 
sensitive  to  errors  associated  with  the  use  of  a  random  wave  transformation  model.  Using  the 
surface  mass  flux  model  to  estimate  UR  results  in  rms  errors  of  5.1  cms"1  and  48  percent, 
which  suggests  that  largest  errors  in  the  vertical  profile  of  mean  undertow  are  related  to 
failure  of  existing  models  to  correctly  predict  the  surface  mass  flux  in  breaking  waves. 

Models  for  the  vertical  structure  of  the  undertow  that  use  an  eddy  viscosity  closure 
to  solve  for  the  turbulent  shear  stress  (Reynolds  stress),  independent  of  the  choice  of 
boundary  conditions  and  depth-dependence  of  the  eddy  viscosity,  result  in  a  general  solution 
of  the  form: 

U(z)  =  U     *  ^±  S(z)  (17) 

pu 

where  Uref  is  a  reference  velocity  that  could  be  either  the  depth-averaged  undertow  or  the 
velocity  at  the  top  of  the  bottom  boundary  layer.  A  nearly  uniform  vertical  profile  of 
undertow,  as  observed  in  the  inner  trough  and  seaward  slope  of  the  bar  (Fig.  1),  requires 
either  no  forcing  (F  =  0)  or  an  infinite  eddy  viscosity  to  be  correctly  reproduced.  An 
example  of  the  spatial  distribution  of  each  term  and  the  total  F(x)  for  the  first  run  of  10  Oct. 
is  shown  in  Fig.  13.  The  dominance  of  the  setup  gradient  term  is  evident,  and  hence  the 
model  predicts  offshore  directed  flow  throughout  the  entire  surf  zone,  which  agrees  with 
observations.  It  can  also  be  noticed  that  modeled  total  forcing,  although  small  offshore  of  the 
bar  and  within  the  trough  where  the  undertow  was  near  uniform,  was  never  nil  within  the 
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cross-shore  transect  occupied  by  the  sled,  which  was  also  the  case  for  all  the  other  runs 
during  the  period  being  analyzed.  Minimum  rms  errors  over  these  regions  are  obtained  for 
large  eddy  viscosity  values  (Table  2),  with  a  maximum  cutoff  level  arbitrarily  set  at 
u  =  0.5 .  These  unrealistically  large  values  of  u  are  associated  with  the  sensitivity  of  the 
eddy  viscosity  model  to  small  errors  in  the  dynamical  forcing,  which  is  calculated  as  the 
difference  between  two  large  numbers  (setup  and  radiation  stress  gradients). 

The  parameterization  for  the  eddy  viscosity  is  of  the  form  (j  ~  /  C,  where  /  and  C 

are  characteristic  length  and  velocity  scales  (Battjes,  1975).  Several  dimensionally  consistent 

H  2 
parameters   for    u    are   found  in   the   literature,    such   as       ( — )/   [Thornton,  1970]; 

h(-)     [Battjes,   1975];    hfgh    [Stive  and  Wind,    1986];    ^^-  {—)m [Haines  and 
P  /        P 

Sallenger,  1984].  An  attempt  was  made  to  relate  best  fit  uwith  these  parameters  after 
excluding  from  the  ensemble  runs  that  reached  the  cutoff  level  (Table  2).  No  statistically 
significant  correlation  at  the  95  percent  confidence  level  was  obtained  between  these 
parameters  and  best  fit  u .  On  the  bar  crest  and  shoreward  slope  of  the  bar  the  mean  undertow 
profile  assumes  a  parabolic  shape  and  best  fit  eddy  viscosity  approaches  a  constant  value 
( |i  ~  0.04 ).  Although  contrary  to  theoretical  expectations,  these  results  suggest  that  the  use 
of  a  constant  eddy  viscosity  across  the  surf  zone  provides  a  reasonable  approximation  for  the 
modeling  of  the  vertical  structure  of  the  undertow. 
Vertical  variation  of  the  eddy  viscosity 

Here,  the  impact  of  both  a  depth  varying  (jt  and  different  boundary  conditions  on  the 
vertical  structure  of  the  undertow  under  field  conditions  are  investigated.  It  is  expected  from 
measurements  (King  and  Kirby,  1994;  Cox  and  Kobayashi,  1997)  that  u,  should  increase 
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from  the  bottom  towards  the  surface  with  increasing  levels  of  wave-breaking  generated 
turbulence.  Several  mathematical  formulations  for  the  vertical  variations  of  u,  were 
investigated  and  the  three  solutions  that  give  best  agreement  with  data  are  described  in  the 
Appendix.  A  comparison  between  the  solutions  for  a  linear  (A6)  and  quadratic  (All)  u  that 
give  the  smallest  overall  rms  errors  for  the  entire  ensemble,  together  with  the  undertow 
solution  (10)  for  a  constant  eddy  viscosity  with  depth  are  shown  in  Fig.  14  for  the  two 
stations  within  each  day  that  have  largest  observed  vertical  structure  of  the  undertow.  No 
significant  improvement  in  the  rms  relative  or  absolute  errors  between  observations  and  the 
solutions  using  variable  \iz  with  depth  is  observed.  These  results  confirm  earlier  modeling 
findings  that  a  depth-dependent  eddy  viscosity  does  not  substantially  improve  the  description 
of  the  vertical  structure  of  the  mean  undertow  (Svendsen  and  Buhr  Hansen,  1988;  Nadaoka 
etal,  1989). 

The  influence  of  the  boundary  condition  choice  to  the  vertical  structure  of  the  mean 
undertow  in  the  middle  layer  is  investigated  next.  Comparison  between  predictions  by  (10) 
that  uses  conservation  of  mass  over  the  vertical  and  the  stress  at  the  trough  level  as  boundary 
conditions  and  (11)  that  uses  a  no-slip  condition  at  the  bottom  to  replace  the  stress  at  the 
trough  level  as  the  second  boundary  condition  are  shown  in  Fig.  1 5  for  the  same  stations  used 
in  Fig.  14.  Again,  no  significant  improvement  in  the  total  rms  errors  between  observations 
and  predictions  by  these  two  solutions  is  obtained,  although  there  are  noticeable  differences 
in  the  predicted  currents  by  each  model  as  a  function  of  depth.  Equation  (10)  shows  better 
overall  comparison  with  data,  while  (11)  better  represents  the  structure  of  the  flow  in  the 
lower  half  of  the  water  column.  Although  models  coupling  the  middle  layer  and  BBL  flows 
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(1 1)  are  expected  to  provide  a  more  realistic  description  of  the  undertow  structure  close  to 
the  bed  (Battjes  et  al.,  1990),  their  applicability  to  field  conditions  is  still  limited  by  the  lack 
of  data  to  constrain  their  free  parameters. 

CONCLUSIONS 

The  predicted  spatial  distribution  of  mean  cross-shore  currents  (undertow)  over  a 
barred  beach  are  compared  with  field  observations  during  the  DUCK94  experiment  to 
quantify  the  relative  importance  of  contributions  from  the  various  terms  in  the  cross-shore 
momentum  equation  and  to  identify  physical  mechanisms  not  yet  incorporated  in  existing 
nearshore  models. 

The  surface  mass  flux  model  based  on  conservation  of  mass  over  the  vertical  and 
including  contributions  associated  with  wave  rollers  is  compared  with  field  observations 
obtained  by  integrating  the  measured  return  flow  over  depth.  It  was  found  that  the  roller 
contributions  to  the  mass  transport  can  be  larger  than  the  contributions  from  the  organized 
wave  motion  when  high  energetic  waves  are  present.  An  increase  of  8  percent  in  the  surface 
mass  flux  prediction  is  obtained  using  linear  wave  theory  combined  with  contributions  from 
wave  rollers  compared  to  the  solution  given  by  non  linear  streamfunction  wave  theory. 

The  largest  discrepancies  between  model  predictions  and  observations  of  undertow 
velocities  occur  over  the  sand  bar,  associated  with  the  failure  of  existing  models  to  correctly 
predict  the  surface  mass  flux  under  breaking  waves.  It  is  postulated  that  the  downward 
mixing  produced  from  breaking-injected  turbulence  penetrating  beneath  the  trough  level 
results  in  a  constricted  return  flow  layer  and  enhanced  undertow  velocities  (undertow  jet), 
a  process  not  accounted  for  in  the  model.  Good  agreement  with  data  is  obtained  for  this 

93 


region  with  an  empirical  relation  that  correlates  the  surface  mass  flux  with  the  mean  rate  of 
wave-energy  dissipation  per  unit  area.  This  empirical  relation  provides  a  simple  alternative 
to  estimate  the  surface  mass  flux  within  the  surf  zone  that  can  be  easily  incorporated  into 
nearshore  current  models  for  field  applications  with  the  same  accuracy  of  existing  surface 
mass  flux  models. 

Setup  (setdown)  is  a  dominant  forcing  mechanism  for  the  undertow  and  is  calculated 
from  the  depth-integrated  and  time-averaged  cross-shore  momentum  balance.  It  is  found  that 
despite  the  dominance  of  the  setup  (setdown)  and  radiation  stress  gradients  in  this  balance, 
the  contribution  from  the  convective  acceleration  of  the  mean  current  is  significant  during 
energetic  wave  events.  It  is  also  shown  that  the  inclusion  of  contributions  from  wave  rollers 
results  in  an  onshore  shift  of  the  point  where  the  setup  begins,  which  has  a  significant  impact 
on  the  dynamical  balance  of  forces  within  the  surf  zone. 

The  vertical  structure  of  the  undertow  is  modeled  using  a  turbulent  eddy  viscosity 
closure,  and  it  is  found  that  a  depth-dependent  eddy  viscosity  does  not  substantially  improve 
the  agreement  with  data  under  field  conditions.  The  effects  of  using  different  boundary 
conditions  to  solve  for  the  vertical  structure  of  the  mean  undertow  is  investigated  by 
extending  the  monochromatic  formulations  of  Stive  and  Wind  (1986)  and  Svendsen  et  al. 
(1987)  to  random  waves  by  ensemble  averaging  over  the  wave  distribution  (Equations  10, 
and  11,  respectively).  No  significant  improvement  in  the  total  rms  errors  between 
observations  and  predictions  given  by  these  two  solutions  is  obtained.  Despite  the  better 
representation  of  the  undertow  structure  close  to  the  bed  obtained  by  using  (11),  which 
couples  the  middle  layer  and  bottom  boundary  layer  flows,  its  applicability  to  field 
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conditions  is  still  limited  by  the  lack  of  data  to  constrain  their  free  parameters. 

In  the  inner  trough  and  seaward  slope  of  the  bar,  the  measurements  show  almost  no 
vertical  structure  for  the  mean  undertow,  which  would  require  a  local  balance  between  setup 
gradient  and  wave  forcing  (F  =  0)  or  a  large  eddy  viscosity  (u  -  °°)  to  be  properly  modeled. 
Although  small,  model  predicted  forcing  was  never  nil  within  these  regions,  so  that 
unrealistically  large  values  of  U  are  required  to  model  the  observed  uniform  profile  of  the 
undertow.  These  unrealistically  large  values  of  u  are  associated  with  the  sensitivity  of  the 
eddy  viscosity  model  to  small  errors  in  the  dynamical  forcing,  that  is  calculated  as  the 
difference  between  two  large  numbers  (setup  and  radiation  stress  gradients). 

On  the  bar  crest  and  shoreward  slope  of  the  bar,  where  strong  wave-breaking  was 
observed,  the  mean  undertow  profile  assumes  a  parabolic  shape  and  best  fit  eddy  viscosity 
approaches  a  constant  value  (u  =  0.04).  These  results  suggest  that  the  use  of  a  constant 
eddy  viscosity  across  the  surf  zone  provides  a  reasonable  approximation  for  the  modeling  of 
the  vertical  structure  of  the  undertow. 

The  absence  of  both  setup  measurements  and  observations  of  currents  close  to  the  bed 
are  major  limitations  of  the  present  study.  Setup  is  an  integral  measure  of  the  dynamic 
response  of  the  mean  water  level  to  cross-shore  gradients  of  momentum  fluxes,  and  its 
gradient  is  an  important  driving  force  for  the  undertow  within  the  surf  zone.  The  lack  of 
current  measurements  close  to  the  bed  prevents  a  quantitative  evaluation  of  the  effect  of 
applying  different  boundary  conditions  to  solve  for  the  vertical  profile  of  the  undertow. 
APPENDIX 

In  this  appendix,  the  solution  for  the  case  of  an  eddy  viscosity  that  varies  linearly  with 
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depth,  as  suggested  by  Okayasu  et  al.  (1988)  is  derived  in  more  detail.  Following  Stive  and 
Wind  (1986),  the  stress  at  the  trough  level  and  conservation  of  mass  over  the  vertical  are 
used  as  boundary  conditions.  The  final  solutions  for  other  depth-dependent  formulations  for 
the  eddy  viscosities  investigated  in  this  study  are  also  included. 
1.  Linear:   \i7  =  a  +  $z! 

(z  +  h) 

In  this  solution,  a  normalized  vertical  coordinate  is  introduced  z '  = -,  and  the 

K 

primes  are  omitted  hereafter  for  simplicity  of  notation.  Assuming  the  undertow  forcing  to  be 
constant  over  depth  and  using  the  eddy  viscosity  closure  for  the  turbulent  shear  stress,  a 
solution  for  the  vertical  distribution  of  the  mean  undertow  is  obtained  by  integrating  (6)  over 
depth,  to  give 

pHzi|£>   =x(z)  =  F(x)z  +  c,  Al 

oz 

where Cj  is  an  integration  constant  that  is  solved  by  evaluating  (Al)  at  the  bottom  (z  =  0) 

C,  =  *b     ■  A2 

Integrating  Al  a  second  time  and  using  A2,  gives 

U(z)  =  ^  [z  -  %  ln(a  +  Pz)]  +  ^  ln(a  +  pz)  +  C2  A3 

PP  P  PP 

where  C,  is  an  integration  constant  that  is  solved  by  applying  conservation  of  mass  over  the 
vertical,  to  give 
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c^u'~-^[l'jc"]-^  c°  M 


where 


a.         tt  +  ftd. 

Co  =  T7  ln( ">  +  ]n<<a  +  K)  ~  1  A5 


substituting  A4  in  A3,  the  final  solution  for  the  vertical  profile  of  the  undertow  is  given  by 

Viz)  =  Ur  +  ^  [  (z  -  -*)  -  |  [ln(a  +  Pz)  +  C0]  ]  +  -^M  [ln(a  +  pz)  -  C0]    A6 
PP  2        p  pP 


2.  Quadratic:    (i    =  a  +  P(z) 


1/2 


F(jc)  H         a2  t(z)  2aJ',/2 

t/(z)  =  U+^  [MAz)  -  — —  -  —  MA  +  -^  [Af2(z)  +  — -f—  -  M3]        A7 
PP  15        p2  PP  P 


where 


2z3/2_  «z+ 2Afzi/2_  2a3 
3  P  p2  P2 


Ml(z)  =  4z3/2  -  ^z  +  i^-z1/2  -  ^L  ln(a  +  pz1/2)  A8 


M7(z)  =  2z  m  -  2aln(a  +  Pz 1/2)  A9 
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M3  =  -L[*Ld?/2  +  a$2dt  +  2a(cc2-p2)ln(a  +  pj//2)  -  a(21n(a)  +  1)]         A10 
P  "f 


3.  Parabolic:    u.  =  a  +  p  [z  -  (z)2] 

This  parabolic  form  for  the  vertical  variation  of  the  eddy  viscosity  was  suggested  by 
Roelvink  and  Reniers  (1994),  and  observed  by  Cox  and  Kobayashi  (1997)  with  LDV  in  a 
laboratory  experiment  on  regular  waves  spilling  on  a  rough,  plane  slope. 


*"(*>   r*T,^       /V3-,,     T(Z) 


tf.  W^  W 


4 


TO  =  ^ +   at  [  *VZ) "  ^r ]  +   «r  W(z)  + "  "f ]  ah 


where 


iV(z)  =  Ptanh-,(-fc£)--|ln[a  +  P(z-z2)]  A12 

ij;  2 


iV2(z)  =  2Ptanh-1(^£2P 


'*"""--    -"^H  A13 


2BJ-6  r      Urrf 

A^3  =  i|/Jr-(2a+p^)tanh"1( J — )  -2atanh"1(-^-)  -— ln[a+P(rff-rf2)]  +  N4     A14 


,   2P</-P  ,    6      iir      a+P(</-d,2) 

JV4  =  (2Prf-P)tanh-'(        f      )  -ptanh^-)  -~ln[ —      '        ]  A15 

ijj  i|j      2  a 


98 


\}i=Vp2+4aP  A16 

In  the  above  listed  solutions,  the  eddy  viscosity  (uj  is  a  function  of  two  non- 
dimensional  parameters  ( a  and  P).  The  value  of  the  eddy  viscosity  at  the  bottom  (  u  _  0)  is 
given  by  a ,  which  in  this  study  is  assumed  to  be  represented  by  the  molecular  viscosity,  due 
to  the  lack  of  reliable  measurements  of  currents  close  to  the  bed  under  field  conditions  to 
better  constraint  this  value.  Although,  the  solution  within  the  middle  layer  does  not 
significantly  change  if  the  assumed  constant  eddy  viscosity  within  the  bottom  boundary  layer 
(Putrevu  and  Svendsen,  1993)  is  used  for  a. 
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Figure  1.  Measured  (  *  )  and  rms  best  fit  model  predicted  (Equation  10)  vertical  profiles  of 
mean  cross-shore  return  currents  (heavy  line)  superposed  on  bottom  profile  with  meantrough 
level  indicated  by  (o)  for  the  entire  period  being  examined. 
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Figure  2.  LBT  "worst  case"  model  prediction  of  H  versus  cross-shore  distance  (dash- 
dotted  line)  compared  with  cubic  spline  interpolation  (solid  line)  of  observations  (  *  ),  first 
run  of  1 1  Oct.  (upper  panel).  The  bottom  profile  is  shown  in  the  lower  panel. 
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Figure  3.  Predicted  (LBT  model)  versus  observed  Hma.  The  solid  line  represents  perfect 
agreement  and  the  dashed  lines  the  ±10  percent  error  bounds. 
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Figure  4.  Observed  vertical  profiles  of  mean  cross-shore  currents  (  *  )  for  the  eighth  run  of 
1 1  Oct.  (a)  and  for  the  fourth  run  of  12  Oct.  (b).  Horizontal  heavy  lines  represent  mean 
trough  level  (ht)  and  vertical  heavy  lines  the  depth  averaged  undertow  ( £/.)  calculated  using 
linear  (dashed)  and  cnoidal  (solid)  wave  theories. 


108 
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Figure  5.  Comparison  between  inferred  surface  mass  flux  ( p  Uobs  ht)  from  mean  velocity 
measurements  using  linear  and  cnoidal  wave  theories  to  estimate  the  mean  trough  level,  hr 
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Figure  8.  Inferred  (p  Uobsht)  versus  linear  wave  theory  predicted  surface  mass  flux  (qw, 
only)  for  the  entire  ensemble  of  22  runs  (upper  panel).  The  lower  panel  shows  the  effect  of 
including  contributions  from  wave  rollers  ( qw  +  qr). 
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Figure  9.  Comparison  between  predicted  surface  mass  flux  (qw  +  qr)  given  by  linear  and 
streamfunction  (Dean,  1974)  wave  theories.  The  dashed  line  represents  perfect  agreement, 
and  the  solid  line  represents  a  linear  regression  with  a  slope  of  0.92. 
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Figure  10.  Inferred  versus  empirically  estimated  (Equation  16)  surface  mass  flux.  The  line 
represents  perfect  agreement. 
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Figure  1 1.  Terms  of  the  momentum  balance  equation  (12)  versus  cross-shore  distance  for  the 
fourth  run  of  12  Oct.  (upper  panel).  The  bottom  profile  is  shown  in  the  lower  panel. 
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Figure  12.  Setup  calculated  using  linear  wave  theory  (dashed  line),  streamfunction  theory 
with  roller  (dash-dotted  line)  and  linear  wave  theory  with  roller  (solid  line)  versus  cross- 
shore  distance  for  the  fourth  run  of  each  day. 
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Figure  13.  Undertow  dynamical  forcing  terms  (6)  versus  cross-shore  distance  for  the  first 
run  of  10  Oct.  (upper  panel).  The  bottom  profile  is  shown  in  the  lower  panel.  The  vertical 
dashed  line  indicates  the  cross-shore  position  occupied  by  the  sled. 
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Figure  14.  Comparison  between  measured  undertow  (*)  and  predictions  using  different 
formulations  for  the  vertical  variation  of  the  eddy  viscosity  (constant  =  solid  line;  linear  = 
dash-dotted  line;  and  parabolic  =  dashed  line)  for  the  third  and  fourth  runs  of  Oct.  10  (a,  and 
b),  Oct.  1 1  (c,  and  d),  and  Oct.  12  (e,  and  f). 
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Figure  15.  Same  as  Figure  14,  but  comparing  predictions  given  by  Equation  10  (solid  line) 
and  Equation  1 1  (dashed  line). 
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Table  1.  Best  fit  model  parameters 


Model 

Wave  Transformation 

Roller 

(LBT  86) 

(LT  97) 

Day 

Run 

a(°) 

7 

error  (%) 

B 

¥ 

error  (%) 

10 

1 

10 

0.33 

2.3 

0.70 

1.3 

2.0 

10 

2 

10 

0.32 

4.9 

0.60 

1.6 

3.8 

10 

3 

10 

0.32 

4.5 

0.60 

1.2 

3.5 

10 

4 

10 

0.32 

4.9 

0.70 

1.1 

5.2 

10 

5 

10 

0.32 

4.1 

0.65 

1.7 

4.2 

10 

6 

10 

0.31 

4.1 

0.75 

1.7 

4.7 

10 

7 

10 

0.31 

5.0 

0.75 

1.9 

4.3 

11 

1 

10 

0.32 

7.0 

0.85 

1.4 

7.5 

11 

2 

10 

0.33 

4.7 

0.65 

1.7 

6.5 

11 

3 

10 

0.33 

4.3 

0.70 

1.1 

7.3 

11 

4 

10 

0.33 

3.9 

0.65 

1.2 

5.7 

11 

5 

10 

0.32 

4.8 

0.60 

1.7 

3.8 

11 

6 

10 

0.32 

4.1 

0.85 

1.2 

7.3 

11 

7 

10 

0.32 

5.5 

0.70 

1.9 

3.4 

11 

8 

10 

0.31 

6.5 

0.85 

1.8 

5.0 

12 

1 

10 

0.34 

5.9 

0.90 

1.3 

3.8 

12 

2 

10 

0.33 

5.1 

0.80 

1.6 

7.8 

12 

3 

10 

0.34 

5.2 

0.80 

1.3 

9.0 

12 

4 

10  . 

0.34 

4.1 

0.80 

1.1 

8.2 

12 

5 

10 

0.34 

4.6 

0.75 

1.3 

7.8 

12 

6 

10 

0.35 

5.8 

0.65 

1.7 

5.3 

12 

7 

10 
mean 

0.34 

5.8 
4.9 

0.95 

1.4   • 

7.0 
5.6 
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Table  2.  Wave  conditions  and  undertow  measurements 


depth- 

rms  best  fit 

rras  wave 

peak 

mean 

tide 

averaged 

eddy 

Day 

Run 

height 

frequency 

direction 

undertow 

viscosity 

Ho 

fP 

Oo 

U0bs 

u 

(m) 

(Hz) 

0 

(m) 

(ms"1) 

(mV1) 

10 

1 

1.70 

0.171 

44 

0.40 

0.06 

0.500 

10 

2 

1.81 

0.162 

38 

0.79 

0.08 

0.500 

10 

3 

1.81 

0.162 

38 

0.98 

0.16 

0.058 

10 

4 

1.68 

0.152 

36 

0.77 

0.14 

0.017 

10 

5 

1.68 

0.152 

36 

0.44 

0.05 

0.500 

10 

6 

1.44 

0.162 

24 

0.04 

0.05 

0.500 

10 

7 

1.44 

0.162 

24 

-0.12 

0.05 

0.500 

11 

1 

2.11 

0.142 

18 

0.22 

0.10 

0.240 

11 

2 

1.88 

0.142 

16 

0.64 

0.12 

0.250 

11 

3 

1.75 

0.142 

17 

0.86 

0.19 

0.033 

11 

4 

1.70 

0.142 

18 

0.92 

0.19 

0.049 

11 

5 

1.70 

0.142 

18 

0.85 

0.05 

0.063 

11 

6 

1.70 

0.142 

18 

0.62 

0.08 

0.500 

11 

7 

1.60 

0.142 

16 

0.32 

0.07 

0.500 

11 

8 

1.60 

0.142 

16 

0.03 

0.15 

0.140 

12 

1 

1.91 

0.162 

18 

0.05 

0.13 

0.130 

12 

2 

2.29 

0.142 

12 

0.34 

0.20 

0.240 

12 

3 

2.29 

0.142 

12 

0.71 

0.30 

0.032 

12 

4 

2.32 

0.142 

10 

0.83 

0.30 

0.036 

12 

5 

2.32 

0.142 

10 

0.84 

0.20 

0.035 

12 

6 

2.35 

0.142 

10 

0.69 

0.07 

0.020 

12 

7 

2.35 

0.142 

10 

0.42 

0.10 

0.500 

H0 ,  fp  ,  and  a© were  measures  in  8m  depth  (FRF  linear  array) 
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IV.  A  QUASI-3D  MODEL  FOR  LONGSHORE  CURRENTS 


(This  chapter  consists  of  work  to  be  submitted 
to  the  Journal  of  Geophysical  Research) 
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A  QUASI-3D  MODEL  FOR  LONGSHORE  CURRENTS 
A.F.  Garcez  Faria,  E.B.  Thornton,  T.C.  Lippmann1,  T.P.  Stanton,  R.T.Guza1,  and 

S.  Elgar2 
Naval  Postgraduate  School,  Monterey,  California  93940 

ABSTRACT 
The  spatial  distribution  of  mean  longshore  currents  over  a  barred  beach  is  examined 
by  comparing  predictions  given  by  a  quasi  three-dimension  model  presented  here  with  data 
acquired  during  the  nearshore  dynamics  experiment  DUCK94  at  Duck,  North  Carolina.  The 
model  includes  forcing  due  to  breaking  waves  described  using  the  roller  concept  (Lippmann 
and  Thornton,  1997),  alongshore  wind  stress,  cross-shore  advection  of  mean  momentum  of 
the  alongshore  current,  and  a  full  non  linear  bottom  shear  stress  with  a  variable  bed  shear 
stress  coefficient,  Cf,  constrained  by  observations  (Garcez  Faria  et  ai,  1997a  [Chapter  1]). 
Contributions  from  the  alongshore  wind  stress  are  mostly  evident  offshore  and  over  the  inner 
trough  of  the  sand  bar  due  to  the  relative  increase  in  the  wind  force  to  wave  force  ratio  as 
wave  forcing  decreases  over  these  regions.  The  advection  of  the  momentum  of  the  longshore 
current  by  mean  cross-shore  currents  is  shown  to  improve  the  agreement  with  observations 
within  the  surf  zone,  0[  10  percent].  The  use  of  a  non  linear  bed  shear  stress  formulation  with 
a  variable  C,  is  shown  to  improve  model/data  comparison,  O[20  percent],  compared  to  the 
use  of  a  constant  C,.  The  largest  overall  improvement  with  observations  is  obtained  by 
incorporating  the  roller  contribution,  O[50  percent]. 


Scripps  Institution  of  Oceanography,  La  Jolla,  California  92903 
Washington  State  University,  Pullman,  WA  99164-2752 
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INTRODUCTION 

The  longshore  current  maximum  observed  in  the  trough  of  a  barred  beach  during  the 
nearshore  dynamics  experiments  DELILAH  and  DUCK94  at  Duck,  North  Carolina,  is  not 
predicted  by  present  theory.  The  simplest  longshore  current  models  balance  cross-shore 
changes  in  the  alongshore  wave  momentum  (radiation  stress)  with  the  alongshore  bottom 
shear  stress.  Waves  break  over  the  bar  reform  in  the  trough  and  again  break  on  the  foreshore. 
Wave  breaking  results  in  changes  in  the  radiation  stress  predicting  two  jets,  one  over  the  bar 
and  the  other  at  the  foreshore,  which  does  not  agree  with  the  observed  current  maximum  in 
the  trough.  An  example  of  the  measured  and  modeled  wave  height  and  longshore  current 
distributions  using  DUCK94  data  are  shown  in  Fig.  1 .  The  predictions  suggest  that  a  transfer 
of  momentum  is  required  to  account  for  the  current  deficit  in  the  trough. 

A  number  of  mechanisms  have  been  proposed  to  mix  momentum  laterally  into  the 
trough  region  to  drive  the  longshore  current.  Traditional  turbulent  mixing,  usually 
parameterized  using  classical  eddy  viscosity  concepts  associated  with  the  shear  of  the 
longshore  current  (e.g.,  Bowen,  1969;  Longuett-Higgins,  1970;  Thornton,  1970),  would 
require  up  gradient  mixing  of  the  longshore  current  which  is  not  feasible.  Battjes  (1975) 
formulated  turbulent  mixing  induced  by  breaking  waves,  where  the  scale  of  the  turbulent 
mixing  is  the  same  order  of  the  wave  height  or  the  local  depth,  which  are  much  too  short 
length  scales  to  explain  the  observations.  Smith  et  al  (1993)  described  wave  breaking  as 
rollers  that  propagate  with  the  wave  at  the  phase  speed;  they  applied  a  turbulent  kinetic 
energy  equation  and  argued  that  turbulence  was  diffused  downward  into  the  water  column 
generating  an  additional  alongshore  thrust.  They  applied  their  formulation  to  the  DELILAH 
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data  and  partially  explained  the  momentum  deficit  in  the  trough. 

Changes  in  the  bottom  shear  stress  due  to  turbulence  being  injected  from  the  surface 
by  breaking  waves  and  modifying  the  vertical  profile  of  the  longshore  current,  which  in  turn 
modify  the  magnitude  of  the  longshore  current,  were  investigated  by  Church  and  Thornton 
(1993).  They  allow  a  spatially  variable  bed  shear  stress  coefficient  dependent  on  the  breaking 
wave  induced  near  bottom  turbulence  levels.  The  model  predicted  cross-shore  profiles  of  the 
longshore  current  improved  agreement  with  observation  compared  with  treatments  using 
constant  bed  shear  stress  values,  but  did  not  completely  account  for  the  momentum  deficit 
in  the  trough. 

Instabilities  of  the  longshore  current  have  been  identified  as  a  mechanism  for 
turbulent  mixing  of  the  longshore  current  originally  suggested  by  Bowen  and  Homan,  (1989). 
Dodd  and  Thornton  ( 1 990)  showed  that  if  shear  instabilities  exist,  there  is  an  accompanying 
cross-shore  mixing  of  momentum.  Putrevu  and  Svendsen  (1992)  carried  out  a  numerical 
study  of  shear  instabilities  over  various  topography  and  using  an  order  of  magnitude  analysis 
concluded  that  even  a  weak  shear  in  the  longshore  current  might  be  capable  of  producing 
significant  mixing.  Significant  shear  instabilities  of  the  longshore  current  have  been  observed 
in  the  field  (Oltman-Shay  et  al,  1989).  Dodd  et  al.  (1997)  calibrated  the  amplitude  of  the 
shear  instabilities  using  field  measurements,  and  then  calculated  the  Reynolds'  stress 
associated  with  the  instabilities.  They  found  that  the  mixing  predicted  due  to  shear 
instabilities  to  be  in  qualitative  agreement  with  that  required  for  modeled  longshore  current 
profiles  to  agree  with  observed  profiles. 

Putrevu  et  al.  (1995),  using  an  order  of  magnitude  analysis,  suggested  that  alongshore 
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pressure  gradients  induced  by  alongshore  variation  in  bottom  topography  could  contribute 
substantially  to  the  forcing  for  longshore  currents.  Reniers  et  al.  (1997)  examined  this  effect 
using  DELILAH  data  and  found  that  even  a  relatively  small  alongshore  pressure  gradient 
acting  in  the  same  direction  as  the  wave  forcing  over  the  trough  of  the  bar,  can  have  a 
significant  effect  on  the  longshore  current  distribution. 

The  advection  of  the  momentum  of  the  longshore  current  by  mean  cross-shore 
currents  as  a  source  of  momentum  mixing  was  suggested  by  Putrevu  and  Svendsen  (1993). 
The  longshore  current  is  strongest  toward  the  surface  and  decreases  to  zero  at  the  bottom. 
The  cross-shore  mean  current  has  an  onshore  transport  in  the  wave  crest/trough  region  and 
an  offshore  transport  below  (undertow).  Svendsen  and  Putrevu  (1994)  showed  that  the  net 
interaction  could  induce  significant  mixing  of  the  momentum  of  the  mean  currents. 

The  objective  of  this  paper  is  to  formulate  a  simplified  quasi-three-dimension  model 
to  describe  longshore  currents  assuming  that  the  bottom  profile  has  straight  and  parallel 
contours.  Included  in  the  formulation  are  forcing  due  to  breaking  waves  described  by  a  roller 
(Lippmann  and  Thornton,  1997)  and  alongshore  wind  shear  stress,  cross-shore  advection  of 
mean  momentum  of  the  alongshore  current  by  the  shear  of  the  mean  cross-shore  current  as 
suggested  by  Svendsen  and  Putrevu  (1994)  and  a  full  non  linear  bottom  shear  stress.  The 
observations  acquired  during  DUCK94  will  be  used  to  test  the  model  and  to  quantify  the 
relative  importance  of  contributions  from  the  various  terms  of  the  alongshore  momentum 
equation  within  the  surf  zone. 

MODEL  DEVELOPMENT 

A  three-dimensional  model  of  the  nearshore  circulation  is  derived  assuming 
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stationary  wave  conditions,  straight  and  parallel  bottom  contours,  and  random  waves  that  are 
narrow-banded  in  both  frequency  and  direction.  A  right-handed-coordinate  system  is  used 
with  x  positive  onshore,  y  positive  to  the  south  and  z  positive  upward  from  the  sea  surface. 
In  the  following,  the  energy  flux  balance  is  considered  first  to  describe  wave  transformation, 
which  includes  contributions  by  breaking  waves  described  with  the  roller  concept.  The 
random  wave  distributions  are  described,  which  are  used  to  ensemble  average  the  various 
equations.  The  conservation  of  mass  is  described  next,  which  provides  an  integral  condition 
for  the  solution  of  the  vertical  profile  of  the  cross-shore  velocity.  The  conservation  of  cross- 
shore  and  alongshore  momentum  fluxes  are  considered  next.  The  depth  integrated  cross- 
shore  momentum  equation  describes  wave  set-up/down,  which  is  the  primary  forcing  for 
cross-shore  currents  (undertow).  The  alongshore  momentum  equation  contains  the  lateral 
transfer  of  mean  momentum  term  that  couples  cross-shore  and  longshore  currents,  which 
requires  specifying  the  vertical  profile  of  U(x,z)  by  solving  the  cross-shore  momentum 
equation.  The  depth  integrated  alongshore  momentum  equation  is  then  solved  to  find  the 
cross-shore  distribution  of  the  mean  longshore  current,  V(x). 
Conservation  of  Energy  Flux 

Wave  heights  across  the  surf  zone  are  determined  using  the  energy  flux  balance 
equation 

4~[E„C  cos(a)]  +  i-[£  Ccos(a)]  =  D  (1) 

ox  s  ox 

where  the  energy  is  partitioned  into  wave,  Ew,  and  roller,^,  contributions  (Svendsen,  1984), 
C  and  C  are  the  wave  group  velocity  and  phase  speed  described  by  linear  theory,  a  is  the 
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mean  incident  wave  angle  with  respect  to  the  shore  normal  calculated  using  Snell's  Law,  and 
D  represents  dissipation  of  roller  energy  (Nairn  et  al,  1990).  Following  Deigaard  (1993), 
the  roller  dissipation  is  modeled  as  D  =  C  xs,  where  xs  is  the  mean  surface  shear  stress  at 
the  wave/roller  interface. 

Both  cross-shore  distributions  of  wave  height  and  the  roller  dissipation,  D,  are 
calculated  using  the  Lippmann  and  Thornton  (1997)  model.  This  model  assumes  that  random 
wave  heights  are  Rayleigh  distributed  everywhere,  even  when  breaking  (Thornton  and  Guza, 
1983)  and  was  shown  to  give  good  results  for  the  DUCK94  data  set  (Garcez  Faria  et  al., 
1997b  [Chapter  2]).  Therefore,  wave  properties  in  the  various  equations  are  ensemble 
averaged  by  integrating  through  the  Rayleigh  distribution,  p(H),  and  using  a  weighting 
function,  W(H),  (Whitford,  1988)  to  calculate  the  percentage  of  waves  that  are  breaking, 
pb(H)  =  W(H)P(H). 
Conservation  of  mass  flux 

The  conservation  of  mass  for  straight  and  parallel  contours  is  given  by 


f-p[U(z)  +  u(z)+u(z)]dz  =  0  (2) 

-h 

where  the  boundary  condition  of  no  flow  through  the  beach  has  been  utilized  and  the 
velocities  have  been  partitioned  into  mean,  wave  and  turbulent  contributions.  In  an  Eulerian 
reference  frame,  there  is  an  onshore  mass  transport  by  waves  and  an  additional  transport  by 
the  rollers  of  the  breaking  waves  (Svendsen,  1984)  limited  to  an  upper  region  between  the 
crest  and  trough.  The  ensemble  and  time-averaged  onshore  mass  transport  in  the  upper  region 
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including  contributions  from  wave  rollers  is  balanced  by  a  mean  return  flow  below  the  trough 
(undertow) 


Ti, 


<q>  =  -<  f  pU{z)dz>  =   ~pUrht  (3) 


-h 


where  x\t  is  trough  elevation,   Ur  is  the  depth-averaged  return  flow,  ht  is  the  depth  below 
the  trough,  and  <  > indicates  ensemble  averaging  (see  Garcez  Faria  et  al.,  1997b  [Chapter 
2],  for  details). 
Conservation  of  momentum  flux 

The  depth  integrated  and  time  averaged  cross-shore  momentum  equation  is  given  by 


dS        dM  dU2  -         dn 

-JEL   +^   +  p(T1+/2)_^_   =    -pg(T1+/l)^n  (4) 

ox         ox  ox  ox 


where  r)  is  the  mean  wave  setup  (down),  S  is  the  wave  radiation  stress  associated  with  the 
momentum  flux  due  to  the  wave  motion,  the  second  term  on  the  lhs  is  the  cross-shore 
gradient  of  the  momentum  flux  associated  with  wave  rollers,  and  the  last  term  on  the  lhs  is 
the  convective  acceleration  of  the  depth-averaged  undertow. 

The  depth  integrated  and  time  averaged  alongshore  momentum  equation  is  given  by 


F    + 


dL       d 


^  \ 


[fpuvdz*  jp(a+u)(y+v+v)dz  ]  =  -tJ  +  t"  (5) 


dx        dx 

-h  n, 


where/7 ,  is  the  ensemble  averaged  alongshore  wave  forcing,  which  is  modeled  here  using 
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a  two-layer  concept  where  waves  are  breaking  (Lippmann  el  ai,  1996).  The  lower  layer 
contains  the  organized  wave  motion  and  the  upper  layer  describes  wave  breaking  as  an 
elevated  body  of  turbulent  fluid  riding  on  the  front  face  of  the  wave  (wave  roller),  consistent 
with  spilling  breakers.  Following  Deigaard  (1993),  the  exchange  of  momentum  across  the 
wave/roller  interface  is  described  by  the  upward  and  downward  volume  fluxes  between  these 
two  layers.  The  downward  flux  produces  a  shear  stress  at  the  surface,!  ,  that  is  associated 
to  the  total  wave  forcing,  F  ,  by  conservation  of  momentum.  Therefore,  following  Stive  and 
De  Vriend  (1994),  F  can  be  written  in  terms  of  the  average  shear  stress  at  wave/roller 
interface  and  local  mean  incident  wave  angle 

dS  dS 

y  dx  dx 


where  S      is  the  alongshore  wave  radiation  stress  and  S      is  the  roller  contribution  to  the 

xy.w  °  xy,r 

radiation  stress.  The  surface  shear  stress  is  calculated  from  the  roller  dissipation,  D,  using 
Lippmann  and  Thornton  (1997)  model,  which  is  calibrated  with  video  records  of  the  cross- 
shore  variation  of  wave  breaking  (see  Garcez  Faria  et  ai,  1997b  [Chapter  2],  for  details). 

The  second  term  on  the  Ihs  of  (5)  is  the  lateral  mixing  associated  with  the  integrated 
turbulent  Reynolds  stress  that  is  parameterized  using  an  eddy  viscosity  model  (Bowen,  1969; 
Thornton,  1970) 


;■ -  dV 

Sxy    =      J    P»vdz      =    -pVtx(h+T])—  (7) 

-h 
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where   v    is  a  kinematic  eddy  viscosity  that  is  estimated  following  Battjes  (1975)  by 

-      D  m 
v    =  N  (h  +  r\)  ( — )     ,  where  N  is  an  adjustable  coefficient  of  0(  1 ),  the  turbulent  length 

P 

D  m 
scale  is  order  of  the  local  depth,  and  the  turbulent  velocity  scale  is  ( — )    . 

P 
The  third  and  fourth  terms  on  the  Ihs  of  (5)  represent  the  lateral  mixing  associated 

with  the  interaction  between  cross-shore  and  alongshore  currents  (Putrevu  and  Svendsen, 

1993),  which    require  specifying  the  vertical  profiles  of  both  longshore  and  cross-shore 

currents  and  will  be  addressed  later. 

Following  Thornton  and  Guza  (1986),  the  time  and  ensemble  averaged  non-linear 

alongshore  bottom  shear  stress,  xb ,  is  modeled  as 


<t*>  =    f  [-  fpCfJv2  +  u2  +  2Vusin(a)   (V  +  usin(a))  dt  ]  p(H)dH        (8) 

0  T 


where  Cf  is  a  bed  shear  stress  coefficient. 

The  time  averaged  alongshore  wind  stress,  xw,  is  calculated  using  the  drag  coefficient 
method 

t;  =  PaCdWWy  (9) 

where  p  is  the  atmospheric  density,  W  and  W ,  are  the  total  and  alongshore  components  of 
wind  speed  and  Cd  is  an  atmospheric  drag  coefficient.  For  Cd,  the  value  obtained  by  the 
WANDI  Group  (1988)  from  observations  at  intermediate-depth  water,  is  enhanced  by  33 
percent  to  take  into  account  the  increased  surface  roughness  associated  with  breaking  waves 
(Whitford  and  Thornton,  1993,  1996). 
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Vertical  profile  of  longshore  velocities 

Garcez  Faria  et  al.  (1997a  [Chapter  1])  using  data  from  DUCK94  showed  that  the 
vertical  structure  of  mean  longshore  currents  on  a  barred  beach  is  well  described  by  a 
logarithmic  profile 


v.         7  +h 
V(z)  =  —  ln(^-^)  (10) 

K  Z„ 


where  k  is  the  von  Karman  constant  (0.4),  v,  is  the  alongshore  shear  stress  velocity,  and  z 
is  an  apparent  roughness  due  to  the  increase  in  effective  roughness  owing  to  the  presence  of 
waves  (Grant  and  Madsen,   1979).  The  vertical  profile  of  the  longshore  current  is 
parameterized  by  v.  and  za ,  and  the  mean  bottom  shear  stress,  ib ,  is  related  to  the  alongshore 
shear  stress  velocity  by  t,  =  pv„2. 

To  compare  model  predicted  longshore  currents  with  measurements,  the  mean 
longshore  current  found  by  integrating  (10)  over  depth  is  defined  as  the  reference  velocity 
given  by 


1      r    v-  ,    ,z+h,  ,         v.  ri    ,  h  ,      Zc 


Vm(x)=±     f  ^\n(±-Z)dz  =  —  [ln(-^)  +  -f  -1]  (11) 

h   -h{Za  K  Za  K  Za  h 

Vertical  profile  of  cross-shore  velocities 

Garcez  Faria  et  al.  (1997b  [Chapter  2])  extended  the  monochromatic  wave 
formulation  for  undertow  by  Stive  and  Wind  (1986)  using  the  random  wave  formulation  of 
Thornton  and  Guza  (1983).  The  vertical  structure  of  the  ensemble  averaged  undertow  is 

133 


obtained  by  depth-integrating  the  cross-shore  momentum  twice  and  solving  for  the 
integration  constants  by  applying  conservation  of  mass  over  the  vertical  (2)  and  the  stress  at 
the  trough  level  (by  integrating  the  vertical  momentum  equation  from  the  bottom  to  crest 
elevation),  to  give 

i        -2  3/z2-n2        "  <T*>  Tl 

<U(z)>  =  Ur  +  —  [^-+hz  + — -2-]   f  Fxp(H)dH  +  — x—[z+h--±]  (12) 

P(J     2  6  o  PM 

where  u  is  a  depth-independent  vertical  eddy  viscosity  coefficient  (Garcez  Faria  et  al., 
1997b  [  Chapter  2]),xb  is  the  mean  cross-shore  bed  stress,  and  F  is  the  undertow  forcing, 


given  by 


„        dpU2        1  dp  rT3    ~i  dr\  .._, 

Fx  =  — —  + -[uz-w*]  +  pg — L  (13) 

dx  2  dx  dx 


The  cross-shore  bottom  stress,  t  ,  has  already  been  shown  to  be  small  (Longuet-Higgins, 
and  Stewart,  1964),  and  thus  the  last  term  on  the  rhs  of  (12)  is  neglected  here. 
Momentum  transfer  by  mean  currents 

Momentum  mixing  is  the  result  of  the  vertical  shear  between  the  cross-shore  and 
longshore  mean  currents  as  described  by  the  third  and  fourth  terms  on  the  lhs  of  (5).  In  the 
water  column  beneath  the  trough,  the  alongshore  momentum  of  the  mean  longshore  current, 
p  V,  is  advected  offshore  by  the  undertow,  giving 


f  pUVdz  =  p<U>— /r[ln(  — )-l]  +  -^-h?<F>  (14) 

{  K  *a  9kVz 
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which  can  be  rewritten  as  a  function  of  the  reference  velocity,  Vm  (11),  neglecting  the  term 

—  «  1 ,  to  give 
h 

pUVdz  =  p<U>htVm  *  —  hf<Fs>  (15, 

J  y  KU 

-h  ^ 

In  the  surface  region  between  the  crest  and  trough  of  the  wave,  the  alongshore 
momentum  of  the  mean  longshore  current,  p  V,  is  advected  shoreward  due  to  the  mean  mass 
transport  velocity  of  the  waves  and  rollers  of  breaking  waves,  which  is  described  by  the 
fourth  term  on  the  lhs  of  (5).  This  term  can  be  simplified  assuming  that  waves  and  turbulence 
are  uncorrelated,  small  incident  wave  angles  and  that  the  radiation  stress  contributions  ( p  u  v) 
in  the  crest-trough  region  are  of  higher  order,  giving 


[p(u+u)(V  +  v  +  v)dz  ]  =  pV(ri=0)  f(u+u)dz  =  <q  >—  ln(— )  (16) 

J  J  K  Z 

Tl,  n,  a 


The  cross-shore  distribution  of  the  mean  longshore  current,  Vm(jc) ,  is  calculated  by 
solving  (5)  using  the  Newton-Raphson  iterative  method,  applying  a  centered  finite 
differencing  scheme  to  rewrite  the  partial  derivatives  with  the  boundary  condition  of  no  flow 
through  the  beach. 

FIELD  DATA 

Field  data  to  test  the  model  were  acquired  during  the  DUCK94  nearshore  experiment 
at  the  U.S.  Army  Corps  of  Engineers  Field  Research  Facility  (FRF),  Duck,  North  Carolina, 
in  October  1994.  For  details  of  the  experiment  see  Thornton  et  al.  (1997),  Elgar  et  al,  1997, 
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and  Garcez  Faria  et  al.  (1997a  [Chapter  1],  and  1997b  [Chapter  2]).  The  data  presented  here 
are  from  10-12  Oct.  when  strong  currents  caused  by  a  storm  with  predominant  winds  and 
waves  from  north  were  present  and  the  vertical  structure  of  both  longshore  and  cross-shore 
currents  have  been  studied  by  Garcez  Faria  et  al.  (1997a  [Chapter  1],  and  1997b  [Chapter 
2]).  During  this  period,  observations  spanned  the  entire  surf  zone  and  conditions  approximate 
the  assumptions  of  steady  state  and  straight  and  parallel  contours. 

Cross-shore  transects  of  wave  heights  and  currents  were  measured  using  a  cross-shore 
array  of  two-component  Marsh-McBimey  electromagnetic  current  meters  (ems,  hereafter) 
and  pressure  sensors  at  each  of  13  locations  (See  Elgar  et  al,  1997,  for  details).  The  vertical 
structure  of  the  current  was  measured  with  a  vertical  stack  of  seven  ems  mounted  on  a 
mobile  sled,  which  was  deployed  approximately  25  m  to  the  south  of  the  fixed  cross-shore 
array.  Additionally,  directional  wave  spectra  were  acquired  using  a  linear  array  of  10  pressure 
sensors  in  8  m  depth. 

For  the  first  run  on  each  day,  the  sled  was  towed  by  the  Coastal  Research  Amphibious 
Buggy  (CRAB)  to  its  farthest  offshore  location  seaward  of  the  bar  (approximately  160  m 
from  the  shoreline).  A  forklift  on  the  beach  pulled  the  sled  shoreward  10  to  30  meters  for 
subsequent  measurement  runs  that  are  referred  to  in  the  text  by  sequential  numbers  within 
each  day.  Each  data  run  was  nominally  one  hour,  and  seven  to  eight  runs  were  made  across 
a  transect  during  each  day  spanning  the  high  tide  during  this  period. 

For  the  cross-shore  array,  longshore  currents  were  measured  near  mid-depth,  but  at 
arbitrary  elevations,  z  •  To  make  comparisons  with  the  model,  measured  longshore 
currents,  V      ,  are  corrected  to  correspond  to  the  reference  mean  longshore  current,  V  , 
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using 


ln(A)  +  f£_i 

V                                              z         h 
V     = V      (z  =  z      )  = V      (z=z      )  cn\ 

Z. 


"a 


where  Vcor  is  the  corrected  velocity.  The  elevations  of  the  ems  relative  to  the  mean  water 
depth,  z  ,  were  adjusted  using  bathymetric  profiles  measured  daily  by  the  CRAB  and  the 
measured  mean  tidal  elevation.  The  apparent  roughness,  za,  at  each  sensor  location  within 
the  cross-shore  transect  occupied  by  the  sled  was  obtained  by  interpolating  the  values 
obtained  at  the  sled  locations  found  by  fitting  a  logarithmic  profile  to  the  measured  vertical 
distribution  of  longshore  currents  based  on  a  linear-regression  least-squares  method  (Garcez 
Faria  et  ai,  1997a  [Chapter  1]).  For  sensors  located  off-shore  of  the  sled  transect,  a  constant 
value  equal  to  the  measured  z  at  the  most  off-shore  position  occupied  by  the  sled  was  used. 
An  example  of  the  measured  and  corrected  mean  longshore  velocities  during  the  first  run  of 
1 1  Oct.  are  shown  in  the  upper  panel  of  Fig.  2.  The  root  mean  square  (rms  ,  hereafter)  value 
of  the  velocity  corrections  for  all  22  runs  considered  is  7.1  cm  s"1  which  results  on  a  rms 
reduction  of  16  percent,  and  maximum  correction  of  16.4  cm  s"1  amounting  to  a  27.5  percent 
reduction  (Fig.  2,  lower  panel).  Data  analysis  showed  that  the  velocity  correction  is  not  a 
strong  function  of  z  ,  but  rather  mostly  determined  by  the  ems'  elevation  and  local  depth. 
If  instead  of  an  interpolated  value,  a  constant  za  equal  to  the  mean  of  all  measured  values 
within  the  days  being  analyzed  were  used,  the  rms  difference  between  the  velocity 
corrections  given  by  these  two  approaches  is  2.1  cm  s"1  (5  percent). 
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RESULTS  AND  DISCUSSION 

The  relative  importance  of  the  various  terms  in  the  depth  integrated  and  time 
averaged  alongshore  momentum  balance  (5)  is  investigated  by  comparing  model  predictions 
with  field  measurements.  Contributions  from  each  term  in  (5)  are  addressed  separately,  first 
by  using  a  constant  bed  shear  stress  coefficient  ( CX  and  then  by  specifying  cross-shore 
variation  of  C,.  Model/data  comparison  are  quantified  using  a  relative  rms  error  defined  by 


8  =  100  [ 


E(V       -  V     ) 
v     meas  mod ' 


\       Ev, 


2 
meas 


2 

]  (18) 


In  model  comparisons  with  data,  the  vertical  structure  of  mean  currents  is  examined 
at  the  single  sled  location  during  each  run,  while  the  cross-shore  distribution  of  mean 
longshore  velocity  is  examined  using  the  array  data  combined  with  the  single  sled  data  point. 
The  runs  used  in  the  following  figures  to  illustrate  contributions  from  different  terms  in  the 
alongshore  momentum  balance  are  arbitrarily  chosen  to  cover  as  much  of  the  data  as  possible 
Alongshore  wave  forcing 

The  importance  of  including  contributions  from  wave  rollers  is  examined  in  a  simple 
balance  between  wave  radiation  stress  and  bottom  shear  stress  using  the  Lippmann  and 
Thornton  (1997)  model  that  includes  rollers.  The  basis  for  comparisons  is  the  Thornton  and 
Guza  (1986)  model  that  does  not  include  rollers,  which  implies  that  dissipation  of  wave 
energy  by  wave  breaking  occurs  locally;  this  model  has  good  success  in  describing  wave 
driven  longshore  currents  on  near  planar  beaches.  A  comparison  of  the  two  model 
predictions  with  observations  of  mean  longshore  currents  for  the  third  run  of  each  day  being 
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examined  is  shown  in  Fig.  3.  Incorporating  the  roller  contribution  into  the  alongshore  wave 
forcing  significantly  changes  the  cross-shore  distribution  of  the  mean  longshore  current, 
resulting  in  a  better  agreement  with  observations  (Tables  1  and  2). 

The  effects  of  including  wave  rollers  is  to  displace  the  forcing  of  the  longshore 
current  shoreward.  The  magnitude  of  this  onshore  displacement  is  controlled  by  the  steepness 
of  the  front  face  of  the  breaker  and  wave  phase  speed  (Nairn  et  al.,  1990;  Deigaard,  1993; 
Thornton  and  Lippmann,  1997).  For  highly  asymmetric  waves,  the  advection  distances  are 
smaller  and  the  longshore  current  maximum  will  be  located  closer  to  the  sand  bar  crest  than 
for  more  symmetric  waves.  Advection  distances  also  increase  with  increasing  wave  phase 
speed,  which  is  proportional  to  water  depth  in  shallow  water,  and  hence  more  current  is 
forced  in  the  trough  at  higher  tidal  elevations.  This  tidal  modulation  of  the  longshore  current 
is  consistent  with  observations  during  the  DELILAH  experiment  (Thornton  and  Kim,  1993). 

This  good  agreement  with  data  suggests  that  the  alongshore  wave  forcing  within  the 
surf  zone  can  be  modeled  using  the  average  shear  stress  at  the  wave/roller  interface,  x  ,  as 
long  as  model  predicted  shear  stress  is  constrained  by  wave  breaking  observations. 

It  can  also  be  noticed  in  Fig.  3  that  good  agreement  is  obtained  between  the  sled  data 
and  the  depth-corrected  data  from  the  fixed  array,  despite  the  approximate  25  meters 
separation  in  the  alongshore  direction  between  the  sensors.  These  results  support  the 
assumption  that  there  was  no  alongshore  variation  of  the  mean  longshore  current  in  the 
vicinity  of  the  measurements  during  the  period  being  analyzed. 
Wind  forcing 

During  the  three-day  period  studied  here  the  mean  alongshore  wind  speed  was  9. 1 
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m  s"1  (range  7.7  -  12.6  m  s"1)  resulting  in  an  average  wind  force  to  wave  force  ratio  of  0.21 
(0.10  -  0.50)  [Chaper  1].  Despite  the  wave  force  dominance  characterizing  wave-driven 
longshore  currents,  wind  force  contribution  is  not  negligible.  The  effect  of  including  the  time 
averaged  wind  stress  contribution,  x  * ,  to  the  alongshore  momentum  balance  is  shown  in  the 
upper  panel  of  Fig.  4  for  the  time  period  with  largest  wind  force  to  wave  force  ratio  (seventh 
run  of  10  Oct.).  Better  agreement  with  observations  offshore  and  within  the  trough  is 
obtained,  which  is  related  to  the  relative  increase  in  the  wind  force  to  wave  force  ratio  as 
wave  forcing  is  decreased  over  these  regions. 

xw  is  a  direct  function  of  the  wind  drag  coefficient,  Cd,  that  is  yet  poorly  constrained 
by  measurements  within  the  surf  zone.  The  constant  value  for  this  parameter  adopted  here 
is  questionable,  as  the  increased  surface  roughness  over  the  bar  associated  with  wave 
breaking  is  theoretically  expected  to  result  in  larger  values  for  Cd  over  this  region.  The 
sensitivity  of  the  model  to  uncertainties  in  Cd  is  shown  in  the  lower  panel  of  Fig.  4  by 
varying  this  value  within  the  range  of  values  obtained  by  Hsu  (1970)  for  the  upper  swash 
zone  of  a  Florida  beach  ( Cd  =  1 .0  -  7.0  x  10"3).  These  results  indicate  that  uncertainties  in  Cd 
could  lead  to  underestimation  of  the  wind  stress  contribution  to  the  alongshore  momentum 
balance.  Hence,  under  certain  conditions  xwcan  be  an  O(l)  term  even  for  wave  dominated 
conditions  inside  the  surf  zone,  especially  within  the  trough  of  a  barred  beach. 
Turbulent  Reynolds  stress  mixing 

Little  difference  is  found  over  most  of  the  cross-shore  transect  covered  by  the 
measurements  between  model  predictions  with  and  without  the  inclusion  of  the  turbulent 
Reynolds  stress  mixing,  S    ,  parameterized  using  an  eddy  viscosity  model(Fig.  5,  upper 
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panel).  The  large  differences  close  to  the  shore  are  not  of  a  physical  nature,  but  rather 
limitations  of  the  dissipation  model  as  the  depth  approaches  zero  at  the  shoreline  (see 
Lippmann  and  Thornton,  1997  for  details). 

The  gradual  change  in  the  predicted  energy  dissipation  given  by  the  Lippmann  and 
Thornton  (1997)  model  (Fig.  5,  lower  panel)  results  in  an  already  smooth  longshore  current 
distribution  that  does  not  require  this  ad  hoc  mixing  to  smooth  current  gradients.  Another 
possible  reason  that  this  turbulent  mixing  term  has  little  effect  on  the  longshore  current 
profile  is  the  small  turbulent  mixing  scale  by  breaking  waves  used,  0(h)  or  0(H  ),  which 
is  much  less  than  the  length  scale  of  the  width  of  the  surf  zone  during  the  days  considered 
0(50-100  h).  Other  turbulent  mixing  scales  order  of  the  surf  zone,  such  as  shear  instabilities, 
may  contribute  to  momentum  mixing  (Dodd  et  ai,  1997),  but  are  beyond  the  scope  of  this 
paper. 

Following  these  results,  and  to  allow  an  independent  evaluation  of  the  importance  of 
the  momentum  mixing  by  mean  currents  to  the  alongshore  momentum  balance,  the  turbulent 
Reynolds  stress  mixing  is  neglected  for  the  remainder  of  the  paper,  by  setting  N  to  zero. 
Momentum  mixing  by  mean  currents 

The  contribution  from  the  momentum  mixing  by  mean  currents  to  the  alongshore 
momentum  balance  is  illustrated  by  comparing  model  predictions  including  and  neglecting 
this  term  with  observations  of  mean  longshore  currents  for  the  second  run  of  each  day  being 
examined  (Fig.  6).  The  inclusion  of  this  type  of  mixing  results  in  a  shoreward  displacement 
of  the  peak  of  the  predicted  mean  longshore  current  improving  the  agreement  with 
observations  within  the  surf  zone  for  either  a  constant  or  a  variable  C,  by  O[10  percent] 
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(Table  2).  Outside  the  surf  zone,  where  both  the  magnitude  and  vertical  structure  of  the 
cross-shore  current  decrease,  this  mixing  effect  is  not  as  important,  although  it  still  leads  to 
significantly  smaller  errors  in  the  case  of  a  constant  C,  (Table  1).  It  can  also  be  noticed  in 
Fig.  6  and  Table  2  that  better  agreement  with  data  is  obtained  towards  the  end  of  the  three- 
day-period.  It  is  postulated  that  this  improvement  in  model/data  comparison  is  due  to 
enhanced  non  linear  mean  current  interactions  associated  with  the  observed  strong  cross- 
shore  currents  caused  by  increased  wave  heights  and  smaller  angles  of  incidence  with  respect 
to  beach  normal  connected  with  the  arrival  of  a  storm  (Garcez  Faria  et  ai,  1997b  [Chapter 

2]). 

The  solution  for  the  vertical  structure  of  the  undertow  used  in  describing  the 
momentum  mixing  has  as  a  free  parameter  the  vertical  eddy  viscosity  coefficient,  u ,  that  was 
kept  constant  (u  =  0.04)  in  both  the  cross-shore  and  over  the  vertical,  which  was  shown  to 
be  a  reasonable  approximation  under  field  conditions  (Garcez  Faria  et  ai,  1997b  [Chapter 
2]).  The  robustness  of  this  approach  is  shown  by  changing  the  values  of  uby  an  order  of 
magnitude  (Fig.  7).  The  amount  of  mixing  is  directly  proportional  to  u  and  the  model  is 
stable  for  increasing  values  of  u  .  However,  reducing  u  by  a  factor  larger  than  5  leads  to 
numerical  instability  as  the  vertical  profile  of  the  undertow  (12)  is  an  inverse  function  of  u  . 
These  results  indicate  that  the  dispersion  mixing  is  enhanced  by  an  increase  in  the  vertical 
gradient  of  the  undertow  (decrease  in  \i ),  which  is  consistent  with  theoretical  expectations, 
and  that  the  model  is  numerically  stable  for  a  wide  range  of  values  of  the  free  parameter  u  . 
Bed  shear  stress 

The  alongshore  forcing  by  winds  and  waves,  including  the  effects  of  wave-breaking 
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modeled  with  the  roller  concept,  plus  the  momentum  mixing  associated  with  non  linear  mean 
current  interactions  is  balanced  by  the  bed  shear  stress,  modeled  with  a  non  linear  quadratic 
formulation  (Thornton  and  Guza,  1986).  Best  fit  values  for  the  constant  bed  shear  stress 
coefficient,  C„  listed  in  Table  1  are  found  by  iteratively  minimizing  the  error,  s ,  between 
measurements  and  the  complete  model  predictions  (waves,  rollers  and  UV  mixing).  Only 
ems  located  shoreward  of  the  mean  breaker  line  (at  approximately  150  m  offshore)  are 
included  in  the  error  estimate  in  order  to  avoid  biasing  the  rms  fit  by  large  errors  from  ems 
outside  the  surf  zone.  Nevertheless,  the  values  of  e  listed  in  Table  1  are  calculated  using  all 
ems  in  order  to  compare  the  effects  of  incorporating  a  variable  Cf.  Due  to  the  non-linear 
formulation  used  in  the  model,  the  values  of  e  for  partial  solutions  of  the  model  (waves  and 
rollers)  are  obtained  by  re-running  the  model  until  the  Cf  value  that  gives  a  rms  best  fit  is 
found  for  each  case.  The  average  values  of  efor  the  ems  only  within  the  surf  zone  are  listed 
in  Table  2. 

Although  a  constant  C,  value  provides  a  reasonably  accurate  and  robust  numerical 
solution  for  the  cross-shore  distribution  of  mean  longshore  currents  within  the  surf  zone,  it 
does  not  provide  insight  into  the  physics  of  the  problem,  nor  does  it  describe  the  circulation 
outside  of  the  surf  zone  well  (Figs.  3-7).  Feddersen  etal.  (1997),  using  measurements  during 
the  entire  DUK94  experiment  balanced  cross-shore  integrated  wind  and  wave  forcing  with 
a  quadratic  bottom  stress  formulation  and  found  that  best  fit  C,  values  within  the  surf  zone 
are  three  times  larger  than  the  values  outside  of  the  surf  zone.  Garcez  Faria  et  al.  (1997a 
[Chapter  1])  equated  bed  shear  stress  calculated  from  logarithmic  velocity  profiles  to  a 
quadratic  bottom  shear  stress  formulation  and  found  that  the  associated  C,  values  varied  by 
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more  than  an  order  of  magnitude  across  the  surf  zone  (0.0006-0.012)  and  were  positively 
correlated  with  rms  bottom  roughness  estimated  from  sonar  altimeter  measurements.  Their 
measurements  show  that  smaller  values  of  Cf  occurred  over  the  bar  associated  with  small 
bottom  roughness  owing  to  the  planning  action  of  the  intense  wave-breaking  in  this  region. 
Larger  values  of  Cf  were  measured  in  the  trough,  where  the  roughness  increased  due  to  the 
presence  of  mega  ripples  aligned  with  the  strong  observed  longshore  currents.  They  found 
the  following  empirical  relation  between  C,  and  apparent  roughness  length  scale  of  the  bed, 
k    =  30  z  ,  normalized  by  water  depth. 

r   =  0.011  (— )"*  (19) 

1  h 

Here,  a  first  attempt  is  made  to  estimate  the  effects  of  a  cross-shore  varying  C,  on  the 
spatial  distribution  of  mean  longshore  currents  by  comparing  model  predictions  using  this 
simple  relation  with  measurements  in  a  diagnostic  approach.  In  order  to  use  (19),  measured  z 
values  for  the  entire  three-day-period  were  linearly  regressed  as  a  function  of  cross-shore 
position  with  small  values  over  the  bar  and  increasing  shoreward.  A  constant  za  equal  to  the 
furthest  offshore  regressed  value  is  used  for  cross-shore  positions  offshore  of  the  transect 
occupied  by  the  sled,  where  no  measurements  were  available  (Fig.  8,  upper  panel).  The 
points  with  large  estimated  errors  ( 8  >  50  percent)  in  either  za  or  Cf  were  not  included  in 
the  regression.  These  regressed  values  are  then  used  in  (19)  to  estimate  the  average  cross- 
shore  variation  of  Cf  within  the  period  being  analyzed  (Fig.  8,  lower  panel).  Despite  the  large 
data  scatter  resulting  in  low  correlation  between  measured  and  regressed  za  (r  =  0.53),  a 
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statistically  significant  correlation  at  the  95  percent  confidence  level  (  r  =  0.96)  between 
predicted  and  observed  C,  is  obtained.  This  is  clearly  not  the  case  if  an  average  value  for  C, 
(dashed  line)  is  adopted. 

The  simplistic  assumption  that  the  cross-shore  distribution  of  the  bottom  roughness 
length  scale  (ka)  remains  unchanged  during  the  span  of  these  three  days  which  have 
significant  changes  in  wave  and  wind  forcing  is  certainly  not  realistic.  Therefore,  it  was 
expected  that  the  magnitude  of  the  estimated  cross-shore  distribution  of  C,  would  have  to 
be  adjusted  for  each  run  by  multiplying  the  results  given  by  (19)  by  a  factor  of  O(l). 
Nevertheless,  no  change  was  required  for  all  the  runs  on  10  and  1 1  Oct.,  which  suggests  that 
the  cross-shore  variation  of  C,  is  not  overly  sensitive  to  changes  in  ka,  and  is  mostly 
controlled  by  bathymetrical  changes  associated  with  tidal  variation.  On  12  Oct.,  all  the  runs 
had  to  be  increased  by  a  constant  factor  of  1.7  to  match  observations.  This  is  not  surprising, 
as  the  storm  waves  on  this  day  increased  bottom  roughness  within  the  trough  (Thornton  et 
ai,  1997)  and  enhanced  the  non-linear  interactions  between  waves  and  currents,  which  are 
the  two  major  factors  defining  the  magnitude  of  ka  (Grant  and  Madsen,  1979).  Increased  za 
and  C,  values  were  measured  on  12  Oct.  (Fig.  8),  although  these  anomalously  large  values 
were  excluded  from  the  regression  due  to  their  large  uncertainties. 

A  comparison  between  model  predictions  using  the  best  fit  constant  C,  and  the  cross- 
shore  variable  C,  given  by  (19)  is  shown  in  Fig.  9  for  the  first  run  of  each  of  the  days 
examined.  These  results,  together  with  smaller  calculated  8  (Table  1)  show  that  despite  the 
crude  approximation  adopted  here,  a  variable  Cf  significantly  improves  the  agreement  with 
observations  outside  the  surf  zone,  O[20  percent].  Even  inside  the  surf  zone,  where  the 
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constant  C,  value  was  adjusted  to  give  the  best  rms  fit  to  observations,  a  similar  or  even 
better  agreement  is  obtained  (Table  2).  The  overall  improvement  of  model/data  comparison 
obtained  by  applying  the  quasi-three-dimensional  model  with  a  cross-shore  varying  C, 
indicates  that  the  bottom  roughness  length  scale  is  an  important  factor  controlling  the  bed 
shear  stress.  It  also  suggests  that  the  empirical  relation  obtained  by  Garcez  Faria  et  al. 
(1997a  [Chapter  1]),  (19),  provides  a  reasonable  estimate  for  the  cross-shore  variation  of  C, 
under  field  conditions,  although  the  prognostic  use  of  this  relation  is  still  limited  due  to  our 
present  inability  to  predict  the  bottom  roughness  length  scale. 

SUMMARY  AND  CONCLUSIONS 

The  predicted  cross-shore  distribution  of  mean  alongshore  currents  over  a  barred 
beach  given  by  the  quasi-three-dimensional  model  formulated  here  are  compared  with  field 
observations  acquired  during  the  DUCK94  experiment  in  an  attempt  to  improve  our  physical 
understanding  of  the  hydrodynamic  processes  governing  longshore  currents  in  the  surf  zone. 
The  model  includes  forcing  due  to  breaking  waves  described  using  the  roller  concept 
(Lippmann  and  Thornton,  1997),  alongshore  wind  stress,  cross-shore  advection  of  mean 
momentum  of  the  alongshore  current  (Svendsen  and  Putrevu,  1994),  and  a  full  non  linear 
bottom  shear  stress  with  a  variable  bed  shear  stress  coefficient,  C,,  constrained  by 
observations  (Garcez  Faria  et  al.,  1997a  [Chapter  1}). 

Incorporating  the  roller  contribution  results  in  a  shoreward  displacement  of  the 
forcing  of  the  longshore  current  controlled  by  the  wave  asymmetry  and  phase  speed.  Using 
the  shear  stress  at  the  wave/roller  interface,  calibrated  with  observations  of  the  spatial 
variation  of  wave  breaking  (Lippmann  and  Thornton,  1997),  to  force  the  longshore  current 
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significantly  shifts  the  longshore  current  maximum  towards  the  trough  of  the  barred  beach, 
resulting  in  largest  improvement  with  observations,  O[50  percent],  for  the  various 
mechanisms  examined. 

Improved  agreement  with  observations  offshore  and  over  the  inner  trough  of  the  sand 
bar  is  obtained  by  including  the  time  averaged  wind  stress  contribution,  which  is  modeled 
using  the  drag  coefficient  method.  This  better  agreement  is  related  to  the  relative  increase  in 
the  wind  force  to  wave  force  ratio  as  wave  forcing  decreases  over  these  regions. 
Measurements  of  wind  stress  within  the  surf  zone  are  needed  to  better  constrain  the 
magnitude  of  the  wind  drag  coefficient,^,  and  to  investigate  the  expected  cross-shore 
variation  of  this  parameter  with  increased  surface  roughness  associated  with  breaking  waves. 

Two  mixing  mechanisms  for  the  mean  longshore  current  were  investigated.  First,  the 
turbulent  Reynolds  stress  mixing  is  modeled  using  an  eddy  viscosity  formulation  (Battjes, 
1975),  and  it  is  shown  that  it  does  not  significantly  improve  the  agreement  with  observations. 
The  failure  of  this  mechanism  is  possibly  related  to  the  scale  used  to  parameterize  this 
turbulent  mixing  by  breaking  waves,  O(h)  or  0(H  ),  which  is  much  less  than  the  length 
scale  of  the  width  of  the  surf  zone  during  the  days  considered  0(50-100  h).  The  second 
mechanism,  the  momentum  mixing  by  non  linear  mean  current  interactions,  is  modeled  by 
combining  the  solutions  for  the  vertical  profiles  of  longshore  and  cross-shore  currents 
(Garcez  Faria  et  ai,  1997a  [Chapter  1]  and  1997b  [Chapter  2]).  This  mechanism  is  shown 
to  displace  the  peak  of  the  predicted  longshore  current  shorewards  resulting  in  a  better 
agreement  with  observations  within  the  surf  zone  O[10  percent]. 

The  use  of  a  non  linear  bed  shear  stress  formulation  with  a  cross-shore  varying  bed 
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shear  stress  coefficient,  Cf,  constrained  by  observations  significantly  improves  the  agreement 
with  observations,  O[20  percent],  compared  with  using  a  constant  C,  value  adjusted  to  give 
the  best  rms  fit  to  observations  inside  the  surf  zone.  These  results  indicate  that  incorporating 
the  effects  of  a  changing  bottom  roughness  length  scale  through  a  variable  Cf  is  an  important 
step  for  understanding  nearshore  hydrodynamics.  However,  a  predictive  relation  for  the 
bottom  roughness  length  scale  as  a  function  of  mean  currents  in  the  presence  of  waves  is 
required  to  accurately  specify  Cf  values. 

The  overall  good  agreement  obtained  between  model  predictions  and  observations 
indicates  that  alongshore  pressure  gradients  were  weak  during  the  period  being  analyzed  (10- 
12  Oct.),  and  hence  the  circulation  was  essentially  two-dimensional.  Therefore,  despite  the 
lack  of  setup  measurements,  the  DUCK94  data  provide  a  detailed  data  set  to  verify  existing 
2D  and  quasi-3D  models  for  longshore  currents  under  field  conditions. 
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Figure  1.  Cross-shore  distributions  of  predicted  (Lippmann,  Brookins,  and  Thornton,  1996 
model)  H  wave  heights  (upper  panel)  and  mean  longshore  currents  (Thornton  and  Guza, 
1986  model)  Vm  (lower  panel).  The  open  circles  are  the  measured  values  on  10  October  1994 
during  DUCK94  experiment.  The  bottom  profile  (not  to  scale)  is  shown  in  the  upper  panel. 
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Figure  2.  Cross-shore  distribution  of  measured  (+)  and  depth-corrected  (o)mean  longshore 
currents,  Vmm,  for  the  first  run  of  11  Oct.  (upper  panel).  Cross-shore  distribution  of  the 
calculated  velocity-correction  using  (17)  applied  to  the  measured  longshore  velocities  in  the 
fixed  cross-shore  array  for  10-12  Oct  (lower  panel).  The  bottom  profile  (not  to  scale)  is 
shown  in  the  upper  panel. 
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Figure  3.  Comparison  between  measured  mean  longshore  velocity  at  the  sled  (o)  and  fixed 
array  (+)  and  predictions  forcing  the  model  with  waves  only  (solid  line)  and  waves  and 
rollers  (heavy  line)  for  the  third  run  of  10  Oct.  (a),  1 1  Oct.  (b),  and  12  Oct.  (c).  Respective 
bottom  profiles  (not  to  scale)  are  shown  in  all  panels. 
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Figure  4.  Comparison  between  measured  mean  longshore  velocity  at  the  sled  (o)  and  fixed 
array  (+)  and  model  predictions  neglecting  (solid  line)  and  including  (heavy  line) 
contributions  from  the  alongshore  wind  stress  for  the  seventh  run  of  10  Oct.  (upper  panel). 
Model  sensitivity  to  the  wind  drag  coefficient,  Cd  (lower  panel).  The  bottom  profile  (not  to 
scale)  is  shown  in  the  upper  panel. 
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Figure  5.  Comparison  between  measured  mean  longshore  velocity  at  the  sled  (o)  and  fixed 
array  (+)  and  model  predictions  neglecting  (solid  line)  and  including  (heavy  line)  the 
Reynolds  stress  mixing  for  the  fifth  run  of  1 1  Oct.  (upper  panel).  Cross-shore  distribution  of 
predicted  dissipation  given  by  the  Lippmann  and  Thornton  (1997)  model  (lower  panel).  The 
bottom  profile  (not  to  scale)  is  shown  in  the  upper  panel. 
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Figure  6.  Comparison  between  measured  mean  longshore  velocity  at  the  sled  (o)  and  fixed 
array  (+)  and  model  predictions  neglecting  (solid  line)  and  including  (heavy  line)  the 
momentum  mixing  by  mean  currents  for  the  second  run  of  10  Oct.  (a),  1 1  Oct.  (b),  and  12 
Oct.  (c).  Respective  bottom  profiles  (not  to  scale)  are  shown  in  all  panels. 
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Figure  7.  Model  sensitivity  to  the  vertical  eddy  viscosity  coefficient,  u ,  compared  with 
measured  mean  longshore  velocity  at  the  sled  (o)  and  fixed  array  (+)  for  the  first  run  of  1 1 
Oct.  superposed  on  the  bottom  profile  (not  to  scale). 
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Figure  8.  Apparent  roughness,  z  ,  versus  cross-shore  distance  (upper  panel).  The  dashed  line 
represents  the  regressed  values  used  in  (19).  Bed  shear  stress  coefficient,  C,,  versus  cross- 
shore  distance  (lower  panel).  The  solid  line  represents  predictions  given  by  (19)  using 
regressed  za  values  shown  in  the  upper  panel  and  the  dashed  line  is  the  average  Cf. 
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Figure  9.  Comparison  between  measured  mean  longshore  velocity  at  the  sled  (o)  and  fixed 
array  (+)  and  model  predictions  using  a  constant  bed  shear  stress  coefficient,  Cf,  (solid  line) 
and  a  variable  C,  (heavy  line)  for  the  first  run  of  10  Oct.  (a),  1 1  Oct.  (b),  and  12  Oct.  (c). 
Respective  bottom  profiles  (not  to  scale)  are  shown  in  all  panels. 
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Table  1.  Best  fit  model  parameters  and  percent  errors  for  all  current  meters 


Constant  Cf 

Variable  Cf 

Day 

Run 

Percent  Error 

Percent  Error 

UV 

Average 

UV 

cf 

Waves 

Rollers 

Mixing 

cf 

Waves 

Rollers 

Mixing 

10 

1 

0.0050 

63 

50 

51 

0.0031 

54 

43 

36 

10 

2 

0.0065 

65 

69 

50 

0.0030 

65 

44 

43 

10 

3 

0.0050 

60 

43 

45 

0.0029 

95 

37 

37 

10 

4 

0.0065 

70 

46 

45 

0.0030 

65 

38 

37 

10 

5 

0.0060 

276 

52 

50 

0.0030 

68 

47 

46 

10 

6 

0.0060 

72 

50 

50 

0.0032 

54 

37 

38 

10 

7 

0.0055 

85 

55 

48 

0.0032 

51 

36 

38 

11 

1 

0.0045 

98 

62 

50 

0.0031 

57 

23 

23 

11 

2 

0.0055 

67 

41 

40 

0.0030 

52 

32 

33 

11 

3 

0.0050 

68 

44 

44 

0.0030 

54 

33 

33 

11 

4 

0.0060 

74 

47 

46 

0.0030 

61 

37 

37 

11 

5 

0.0070 

78 

54 

48 

0.0030 

62 

43 

44 

11 

6 

0.0065 

90 

61 

45 

0.0030 

62 

32 

33 

11 

7 

0.0060 

65 

58 

54 

0.0032 

54 

37 

41 

11 

8 

0.0090 

83 

64 

59 

0.0031 

63 

35 

35 

12 

1 

0.0090 

74 

53 

41 

0.0030 

62 

33 

23 

12 

2 

0.0085 

74 

45 

20 

0.0053 

53 

30 

28 

12 

3 

0.0090 

66 

54 

35 

0.0052 

56 

37 

29 

12 

4 

0.0065 

64 

48 

35 

0.0050 

56 

29 

27 

12 

5 

0.0070 

66 

44 

40 

0.0050 

60 

36 

36 

12 

6 

0.0090 

72 

62 

46 

0.0050 

62 

48 

45 

12 

7 

0.0100 

68 

139 

98 

0.0052 

67 

137 

105 

Mean 

0.0068 

82 

58 

47 

0.0036 

61 

41 

39 
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Table  2.  Average  percent  errors  for  current  meters  located  within  the  surf  zone 


Constant  Cf 

Variable  Cf 

Day 

Percent  Error 

Percent  Error 

Waves 

Rollers 

uv 

Mixing 

Waves 

Rollers 

UV 
Mixing 

10 

110 

47 

37 

73 

42 

38 

11 

89 

52 

34 

65 

40 

34 

12 

78 

52 

29 

65 

39 

34 

Mean 

92 

50 

33 

68 

40 

35 
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